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FURTHER RESULTS ON ‘THE ANALYSIS 
OF THE VARIATION OF LATITUDE 


A, M. Walker and Andrew Young 
(Received 1957 January 1) 


Summary 


In this paper we give further results of our analysis of the observations of 
the variation of latitude. We find that the best results are those derived from 
the unsmoothed series of values given at monthly intervals by the International 
Service. 

There is evidence of inhomogeneity in the data because of which it is 
difficult to obtain a reliable estimate of the damping factor. We discuss the 
bias introduced into the estimates of the period and damping by the inhomo- 
geneity of the data and by other causes. 





1. In our previous paper (1) we formulated a statistical model of the variation 
of latitude and gave a detailed analysis of its properties. We obtained a 
computational procedure which we applied to a 35-year series of observations 
given by the International Service. ‘These particular observations were 
examined because they were unsmoothed and we hoped that the effect of 
smoothing might be seen when the results were compared with those obtained 
by Jeffreys (2) who used smoothed observations. However, the two sets of 
observations differed in another important respect—the z-term had not been 
taken into account in the set which we analysed but had been in the set analysed 
by Jeffreys. 

In this paper we give much more detailed results and examine some of the 
points raised in the previous paper. 

2. Through the good offices of Professor T. Nicolini of Naples we have 
obtained much more extensive sets of observations, one set of which forms 
series of unsmoothed results covering the interval from 1g00 to 1955 at monthly 
intervals, the other forming series of smoothed results for the interval from 
18g0 to 1955 at 1/10 yearly intervals. In both sets the z-term has been removed. 
These sets are collated from several sources and, as some of the results are as yet 
unpublished in easily accessible form, they are given in ‘Tables Ia and Ib. 
The first series (Table Ia) are exhibited graphically in Fig. 1 (a,c). 

Observations of the variation of latitude were made at a number of 
uncoordinated observatories prior to the inception of the International Service. 
The results were collated by Albrecht and form the basis for the values given 
in ‘Table 14 for the interval 1890-0 to 18g9g9°8. ‘The international programme 
started at the end of 1899. Until 1922-7, the work was under the direction 
firstly of Albrecht and then Wanach and Mahnkopf at Potsdam. From 1922-7 
to 1934°9 the work was directed by Kimura at Mizusawa; he was succeeded 
in 1935 by Professor Carnera and in 1948 by Professor Cecchini of Turin, the 
present director. ‘lhe methods of reduction used by the different directors 
have not been the same and according to Melchior (3, p. 38) the results require 
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corrections; the published results can, he states, be brought into agreement by 
applying the following corrections* 


w =X, +0"-048 = x, + 0":040 = x, + 0004, (1) 
Vw=VE+0"'072 =, + 0":060 = y,, — 0-002, 


where (x, y) are the coordinates of the polet as published and w, k, c, ce refer to 
Wanach, Kimura, Carnera and Cecchini respectively. As the values of x and y 
seldom exceed +0"-300, these corrections seem to be large. 

3. For the international programme, observations are made at several 
observatories on the parallel of latitude 39°8’N. Ideally there should be a 
large number of stations distributed symmetrically on the parallel. However, 
at most, six stations have been in use, and for long periods only three have 
functioned. The stations with details of the intervals during which they 
functioned are given in Table II. 

If Ad; is the value of the variation of latitude of the station, 7, measured 
from its mean latitude, then (x,y) are given by 


x cosA;+y sin A;= Ad; — (2) 


where A; is the longitude of the station and z is Kimura’s z-term which arises 
from a number of causes which are not all perfectly understood. These causes 
include errors in the adopted values of the proper motions and declinations of 
the groups of stars observed, and also perhaps local effects such as refraction 
and oscillations of the vertical. 

When more than three stations were included, equations (2) were solved 
by least squares. For example, from 1go1-7 to 1906-0 the values of x were 
found from the formula 


x= —0°437A¢, + 0'139A¢g, + 0°426A¢, + 0-101 Ag, + 0°:042A¢g, — 0°272Ady, 


whereas, when only three stations were used from 1922-7 to 1935-0, the formula 
became 

x= —0°396A¢, + 0°591Ad; — 0'195Adgg. 
There are evidently quite considerable changes in the relative weights of the 
contributions of different stations. 

Because of precession and proper motion, no individual star can be used in 
the observing programme indefinitely and the programme has been changed on 
a number of occasions. From 1899°9 to 1905-9 each of the twelve groups of 
stars included eight pairs, of which two were “ refraction-pairs” of larger 
declination than the others. At 1906-0 these were replaced by normal pairs, 
and, in effect, the number of stars used was thus increased from 144 to 192, of 
which 132 were in the original programme. Further changes were made at 
1912°0, 1922°7 and 1935-0. At each change the number of new stars introduced 
was not great. In fact, in the programme used since 1955-0, there are still ten 
of the original pairs of stars used since 1899 (see Cecchini, 4). From con- 
versations which one of us (A. Y.) has had with Professors Cecchini and 

* If these suggested corrections are valid they should agree with the means of c, c’ given by 
Jeffreys (2, p. 146.) The latter however average to 0”-o1 at most over the intervals mentioned. 

+ In (1) (x, y) are the coordinates of the pole referred to a left-handed set of axes. We follow 
Jeffreys, and use a right-handed set with /=x and m=—y. We use (x, y) to denote the free com- 
ponent of the motion after the systematic, annual component (/s, m,) has been deducted from (/, m)— 
cf. z, equation (6). 
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Taste la Results of Variations of Latitude : (Mi 
l 
Jan. Feb. Mar. Apr. May June July Aug. Sep. Oct. Nov. Dec. 
1899 7 
1900 69 49 40 . ~?s - G6 -—f -4e - PRP -& - & - SB 
I 8 44 53 68 102 94 143 9! 47 -10 -—7! — 99 
2 —102 — 76 — 48 — 15 71 128 203 216 167 53 — 34 —I119 
3 —162 —199 —136 —103 — 46 37 95 165 226 175 16 — 3 
4 —104 —153 —172 —163 —161 —102 — 14 67 133 183 153 117 
1905 30 — 30 —130 —122 —146 —134 — 90 — 27 61 115 144 132 
6 44 12 — 20 —42 —094 —117 —110 — 9 — 36 — 9 20 40 
7 66 yp | 52 67 66 10 I — II — §7 —147 — 322 — 77 
8 ee 35 99 153 205 223 232 177 79 —- 7 -—209 —246 | 
9 —276 —223 —193 — 85 43 189 280 347 309 163 — 45 —166 
1910 —205 —253 —307 —277. —216 — 78 96 244 324 322 230 95 
1 — 47 —-™%45 —197. —237 —207 —172 —126 25 144 249 314 263 
2 148 73 9 —72 %—114 —154 —116 — 96 — 64 — 16 82 116 
3 116 143 138 136 136 123 108 8 — 26 —106 — 83 — 72 
4 — 8 — 36 74 144 196 182 198 168 92 8 —115 —158 
1915 —196 —158 —150 — 34 24 107 205 270 289 183 71 — 86 
6 —144 —168 —191 —186 —120 — 2 80 191 289 310 222 116 
7 —- 5 — 38 —-112 —172 —126 — 70 — 30 59 160 174 172 160 
8 88 -—- 9 — 28 — 43 — 66 — 76 — 74 12 58 82 99 100 
9 112 68 52 92 80 86 70 50 —- 9 — 58 — 68 — 42 
1920 21 59 28 126 152 198 212 203 170 82 22 — 72 
I — 76 —102 — 62 12 65 176 214 244 236 146 72 22 
2 — 82 -—-104 — 90 — 38 2 92 172 254 210 167 95 — 26 
3 — 89 —125 —178 -—-127 — 90 —I1 79 176 207 198 160 72 
4 —- ¢§ -—-G@ -—t ~8 =F = B - g 84 94 109 94 71 
1925 oe -4 -8 - FP - 8M - MB KC 54 58 76 -  -)\ 4 
6 16 — 75 — 99 — 95 — 95 — 61 — 8 16 54 52 16 28 
7 sae 7a = % 39 21 24 36 59 37 19 . ae. 
8 —- 75 —-7o — 68 —22 — 13 — 12 6 56 —-11 —59 — 53 — 73 
9 — 58 —-116 — 44 33 46 101 110 117 75 — 2 ~~ 2 — 97 
1930 —146 —148 — 88 — 50 5 79 154 173 127 60 II —121 
I —15t —148 —-172 — 95 — 4! 70 165 193 206 72 25 — 56 
2 —168 —238 —186 —I111 —119 — 69 15 137 187 159 100 36 
3 — 85 -—120 ~—138 —162 -—-117 — 91 — 28 68 73 74 55 27 
4 — 26 — 72 -—130 —167 -—123 -—- 74 — 11 37 28 51 20 37 
1935 — 2 9 11 8 II 23 50 42 10 — 22 — 46 — 60 
6 —67 —62 —44 — 8 45 92 102 89 50 ae 
7 —129 -—-133 -—- 8 — 12 33 89 128 139 75 370 =~ 233 — 
8 —134 —125 —102 — 66 — 22 37 129 176 172 100 36 — 76 
9 —im —155 112 -@ -—- %4 36 100 172 179 163 105 22 
1940 — 42 -—92 —122 —133 — & 8 65 99 121 120 58 22 
I 55 15 3. = 7 16 19 57 67 79 93 23 
2 eis = 2 12 18 3 —-24 —12 13 ° 22 
3 43 72 107 98 147 132 98 120 38 — (23 7 — 67 
4 ca seme 10 76 191 278 274 279 I5t — 19 — 55 —152 
1945 “353 19s “197s = OS 92 171 234 279 305 239 100 74 
6 “1g —176 193 ~—195 =& — 39 114 166 270 283 244 208 
7 73° -— 33 ~+~106 —183 -—194 —129 — 20 73 148 226 254 253 
8 245 198 164 79 — 18 —75 — 58 — 74 —112 — 48 — 14 120 
9 124 234 240 238 205 204 m5 — 5 —101 —144 —178 —200 
1950 —— ~~ 76 169 266 352 370 307 1332 — 18 —185 —281 
I —324 —328 -199 — 55 139 259 397 454 377 259 95 —112 
2 —308 —369 -420 —247 —158 45 216 364 428 415 289 66 
3 — 80 —199 —336 —353 —351 —192 — 42 120 280 336 298 126 
4 — 12 — 93 —186 -—222 -—269 —144 — 80 27 34 7° 88 4 
Sources: 1899 Dec.—1905 Dec. Albrecht and Wanach : 
1906 Jan. -1911 Dec. Wanach I 
1912 Jan. —-1922 July Wanach and Mahnkopf I 
1922 Aug.—1934 Dec. Kimura I 
1935 Jan. -1940 Dec. Carnera y ri 


1941 Jan. —1948 Dec. Carnera 


1949 Jan. -1954 Dec. Cecchini 
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-158 
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116 
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(Monthly interval) : 


Jan. 


<a 
— 61 
— 30 
17 
116 
147 
40 
— 94 
—211 
eden 
190 
280 
156 
— 70 
—122 
I4 
231 
298 
190 
— 53 
— 188 


25 
51 
50 
18 
84 
125 
89 
10 
- 29 
67 
95 
117 
25 
— 
—134 
—103 
19 
158 
172 
— 30 
—259 
— 267 
—118 
141 
225 
189 


Feb. 


33 
— 26 
— 94 
= 2 

66 

150 

51 
— 47 
— 208 
—179 

74 

235 

210 
— 44 
—116 

124 

253 

118 
—112 

106 

109 

153 


a 
— 59 
4 

67 
71 
25 

— 
— 61 
40 
7° 

8 

— 20 
— 144 
53 
86 
132 
54 
—228 
— 347 
—241 
171 
224 


unit 0” oor. 


Mar. 


35 

— 147 
—126 
7 
122 
55 
= 23 
— 206 
a7 
— 46 
130 
186 
= we 
— Say 
—198 
= so 
188 
61 

— 82 
—176 
— 80 
12 

— 10 
85 
46 
59 
—124 
— 109 
— 91 
3 

= <3 
17 

— 74 
~ 
Te 
37 
10 
29 
— 68 
—tS3 
—148 
4 

g2 
63 
is? 
— 339 
— 366 
—210 
69 
240 


Apr. 


54 
— 48 
— 169 
mt 
— 56 
102 
31 


= 


— 167 
—298 
“557 

32 

189 

54 
—270 
—108 

84 

24 
— 23 
—I114 
~~ a 
—174 
— 22 

27 
— 6 

31 
— 8 
—115 
—129 
—229 
— 180 
— oF 
“= 
— 2 
— 81 
— 11 
—119 
— 2 
— 2 
_ 31 
— 45 
—I41 
—219 
—118 
— 6 

123 
— 69 
—274 
— 460 
— 348 
— $3 

141 


May 


—440 
—210 
87 
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Resultate .... 
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June 


80 
33 
—184 
—234 
—189 
— 21 
— 22 
56 

= 2 
—290 
— 320 
— 197 
54 
114 
ae 
—224 
—276 
— 98 
— 94 
4 

= 40 
— 109 
—190 
— 187 
—110 
— II 
— 64 
— 
— 89 
—152 
—238 
—201 
— 188 
—130 
— 21 
~~ 2 
— 100 
— 145 
—126 
—140 
=< on 
— eo 
ie 
amen go 
—236 
— 267 
—161 
62 
67 

— 100 
—=—3e7 
— 468 
— 328 
— 62 


m 
Jul. 


86 
73 
—I122 
—238 
—117 
77 
40 
—210 
— 349 
—3lr 


66 
— 12 
—210 
— 282 
—220 
— 139 

14 
— 21 
= 33 
— 140 
—147 
— 102 


— Bo 
— 40 


—113 
—198 
—229 
— 193 
—177 
— 21 
— 18 
— os 
— 96 
Sar 
—109 
—112 

42 

22 
—I119 
—284 
—190 

40 

86 
— 2 
—201 
— 445 
— 400 
= 37s 


Aug. 


63 
143 
— 43 
—197 
—193 
—178 
—101 
79 
132 
— 61 
—314 
— 346 
—124 
102 
79 
— 92 
— 184 
—201 
—126 
57 
62 
— 61 
—158 
—114 
10 

— 88 
6 

6 

61 
—134 
— 188 
—173 
—141 
20 

4 

- a 
— 
— 56 
— 136 
— 
32 
go 
—257 
—229 
— “4 
138 
117 
— 12 
— 352 
— 391 
—259 


Band III, pp. 227-228. 
Band V, pp. 184-185. 
Band VI, pp. 219-221. 
Vol. VIII, pp. 170-172. 
Definitive results communicated to Dublin meeting of LA.U. 1955. 


19 
130 
37 
— 91 
— 142 
—192 
—129 
33 
170 
95 
— 163 
—311 
— 134 
109 
136 


—17!1 
—122 
14 
143 
54 
109 

- a 
~ ee 
64 


64 
17 
99 
47 
—150 


— 96 
—100 


113 
125 
55 
— 61 
— 162 
—131 
— 26 
167 
218 
68 
—190 
— 166 
45 
210 
204 
82 
— 42 
— 66 
— 10 
156 
150 
204 


-— 12 
43 
14 
5° 
34 
104 
89 
102 

~~ ae 


74 
74 
81 


37 

35 

— 68 
—10gG 
180 
183 
— 
—176 
—161 
35 
229 
238 


—260 


123 


Nov. 


— 58 
82 
115 
149 

6 


—130 
—113 
120 
284 
192 
21 
—170 
— 32 
129 
268 
222 
24 
= 
— 29 
154 
107 
185 
83 


83 
— 43 
89 


137 
92 
114 
105 
14 
17 
19 
81 
go 
130 
85 


— 102 
134 
217 
25 

—158 

—247 

—114 
200 
25 
119 

=e 

—271 


2, p. 18 and no. 6, p. 44. 


1955- 


Dec. 


— 61 
— 26 
40 
71 
146 
79 
ae 
—147 
—177 


266 
286 


—117 
— 89 


196 
251 
98 
— 68 
—107 
77 
70 
233 
III 
82 
81 
44 
62 
— 6 
100 
101 
126 
131 
71 
oe 


33 

92 
127 
140 
80 

=— 
— 137 
78 
187 
87 

— 103 
—249 
—216 
29 
246 
229 
41 
—201 


$¢ 
8s 


wu pwn 


190 


191 


wi &WNR OO OND 


1g!I5 


192 


= 
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| 124 TaBLe 1b Results of Variation of Latitude 
' l 
‘°° 1 2 3 4 5 6 or | 8 ‘9 
1890 —271 —258 —188 — 73 55 154 216 202 90 — 55 
I —200 —260 -—260 —210 — 70 65 184 257 254 186 
; 2 61 — 72 —189 —249 —228 — go 61 I 179 158 
3 104 “4 =-C -% -— -I= = <6 9 94 144 
i 4 126 84 s -F -Bo = ee lUSrR UO UUme 
i 1895 — 33 ;-7 -8 8 43 71 Ss -8g -— 
| 6 —115 —III — 8 — 23 53 139 181 162 gl — 8 
| 7; -—F - -— -37F - 7 51 153 204 225 148 
H 8 53 — 65 —172 —202 —174 —105 27 158 199 151 
: 9 gI 26 — 38 — 70 — 8 —104 — 80 — 17 63 12 
: 1900 39 60 ww - &§ -2 -2 ~g2 -G@ -@& -g 
I —- 1 26 58 83 106 122 88 30 — 25 — 80 
i 2 —102 — 7 — 51 38 134 205 200 139 44 — 56 
3 —144 —-179 —142 — 68 29 114 188 209 171 86 
: 4 — 43 —146 -—170 —162 — 94 I 87 151 181 148 
1905 92 — 11 -—I12t —144 —131 — 82 I 84 122 148 
6 97 24 —-26 — 71 —110 —109 — 68 — 33 — 2 17 
i 7 51 64 57 57 %@ - 4-29 —8& 135 —115 
: 8 — 63 15 112 183 226 215 150 33 — 104 —212 
9 ———e ~645 —175 — 6 174 295 338 277 122 — 56 
1910 ~—185 —254 —295 —239 — 73 121 275 323 300 199 
I 36 -—117. —208 —227. —167 — 79 59 175 251 297 
3 246 102 — 22 —109 ~—140 —123 — 91 — 44 20 5 
3 121 136 141 137 123 85 8 — 60 — 99 — 100 
4 — 80 — 20 81 171 197 192 151 7o — 24 —112 
1915 —174 —181 — 109 ° 115 223 284 271 164 28 
6 —102 -—167 -—103 —152 — 43 95 231 306 301 210 
7 85 — 30 —130 —151 — 80 ° 89 159 182 174 
8 131 26 — 33 —61 — 76 — 51 26 66 87 101 
9 108 " 63 89 87 69 390 — 27 — 6 — 63 
| 1920 — 2! 28 80 138 198 213 197 147 72 — 10 
] I — 76 —103 — 55 42 153 238 250 210 140 61 
i 2 — 28 —100 — 80 — 8 92 191 272 217 167 75 
i :_ “2 -. a ee hl 19 134 199 199 148 
i 4 46 — 36 — 94 -102 —90 — 55 44 94 107 go 
3 1925 3S ~- 4¢ ~~ — 6 — 27 — 22 24 60 I 
' 6 8 -— 32 —90 —97 —87 — 43 5 45 pd a 
' 7 15 32 27 21 25 27 49 42 24 6 
: 8 ae io. i we a 34 6 —49 — 62 
i 9 — 67 — 88 — 61 14 64 104 114 86 — 3 — 65 
' 1930 —323 — 147 —102 — 1 17 104 165 139 63 — 6 
: I —130 —150 ~—162 —9g95 — 23 104 180 202 77 II 
' 2 —- & -200 -197 ~—132 — 9% — 39 85 174 160 89 
3 -— Fs -— -% -S -—87 - © 27 72 74 5° 
4 9 — 49 —-116 —167 -116 — 54 14 36 40 33 
1935 wo a 4 12 8 II 52 20 19 —20 — 46 
6 — 63 = =< ee eee 49 98 98 67 12 — 37 
7 “ee Se — se. — 40 102 140 94 35 — 25 
8 —11z2. —133 —-110 — 70 — 18 62 160 180 108 28 
9 — 8§ -110 -116 —55 — 10 54 140 179 166 98 
1940 2 — 68 —-114 —134 — 76 28 84 115 122 
I 70° 38 $.--—-] = @ 21 40 66 82 oe 
2 23 7 — 4 1 9 9 —- 7 — 8 2 14 
3 30 56 87 118 145 129 104 63 7 — -— 
4 — 66 — 26 4 60 210 268 272 182 — 10 — 8&4 
1945 —144 39 — 162 — 60 92 C 
ee ee i a a a 
2 170 1° -% —183 —199 — 8 32 330 : 6 356 
8 251 224 170 99 — 28 —70 — 45 som ge a 
9 a 178 244 248 220 160 4 = 74 —I15§2 an 
1950 — wn js 
1. = 2 SS eS 
: a <a met _ 302 445 408 260 58 
35 392 298 135 104 302 420 420 254 
3 io +) —-317 .—32 —38E. <1 S 6 8 
4 100 © -100 —180 <2 be 4 245 33 253 
40 140 — 40 ° 40 7° 


Sources: 1890-0-1899°8 Albrecht 
1899°9-1905'9 Albrecht and Wanach 
1906-0-1911'9 Wanach 
1912°0-1922'6 Wanach and Mahnkopf 
1922°7—1934'9 Kimura 
1935°0-1940'0 Carnera 
1940°0-1948-9 Carnera 


1949°0-1954°9 Cecchini 














(o'1 year interval) (unit 0”:001) 


m 
Te) “rl 2 3 *4 5 6 "7 8 
83-110 —246 —279 —221 —126 22 144 232 
240 100 — 50 —190 —240 -—166 — 34 95 157 
134 71 19 — 4 — 50 —156 —199 —160 — 32 
107 12! III 87 43 — 11 — 71 —I1! —112 
— 2 = 12 43 95 123 118 70 10 
— 87 -117 — 87 — 28 37 106 101 110 go 
— 29 —102 —139 —153 —140 — 95 II 110 222 
117 22 — 83 —143 —181 —171 —139 — 73 —- 7 
147 107 55 — 10 — 87 ~—146 —176 —140 — 65 
99 150 140 98 25 — 45 —-11I9 —136 —118 
(90 15 40 57 7 76 5° Ss =-* 
— 47 —. §2 — 54 — 28 27 100 141 130 110 
II — 77 —149 — 186 —180 —110 — 22 61 118 
51 — 22 —-110 —IgI —240 —228 -—165 — 51 76 
141 86 —- 5 — 97 — 183 —21!I —175 —123 — 55 
101 144 125 54 — 38 —127 —182 —189 —145 
10 50 45 19 —- 20 — 7 —113) —127. =—130 
—119 -— 60 — 16 15 53 75 67 22 — 4! 
—192 —215 —198 —124 — 17 67 14! 169 160 
— 15 —I52 —261 —312 —290 -—174 — 8 127 228 
237 “o = -—a -—I oe 6 OC = 81 
288 239 122 — 36 —205 —32I —342 —286 —158 
112 199 200 138 55 — 31 — 109 — 160 —I171 
—100 — 48 10 57 g! 109 III gI 24 
— 99 —116 —I114 —I101 — 66 —- !I 80 170 210 
23 — 92 —-205 -—282 -—260 -—175 — 63 70 211 
235 129 I —132 —260 —273 —144 — 10 112 
283 271 154 ao - 99 ~. —7 -—t35 = — & 
123 132 55 —- 20 — 93 —-132 —129 —100 — 70 
—60 —57 — 38 — 3 22 30 30 20 — 9 
—124 —159 —154 —119 — 55 20 go 149 165 
49 — 29 — 81 —110 — 99 — 59 3 80 130 
39 — 40 —117. —185 —181 —128 — 60 98 187 
199 128 34 — 71 —172 —13724 —I5I — 74 4 
117 118 18 — 36 — 87 —107 —106 — 69 — 24 
82 130 100 30 o —-26 — 8 37 60 
109 142 64 j-s -8 = +2 14 
28 38 44 32 — 22 — 40 — 14 42 65 
64 29 — 6 — 42 — 72 — 45 — 18 16 33 
— 10 — 24 — 62 —112 —130 — 65 24 81 113 
$83 —-28 -—107 —143 —178 —139 — 79 19 83 
73 — 22 — 99 %—-207 —228 —216 —161 — 76 go 
112 43 -— 69 —168 -—190 -—212 -—-207 —159 — 29 
129 96 1% — 94 —148 -—190 —182 —118 — 41 
75 76 yy -S =F -s. =. eS. . =o 
8 20 23 I — 16 — 22 — 20 — 12 6 
3 -s —- PP - se — FR — ee I 44 72 
20 — 40 — 8 —110 —I113 — 8 — 18 34 71 
88 29 — 60 —I1I5 —149 — 133 — 71 — 2 73 
119 4 +$-14 — 61 —107 —130 — 99 — 16 33 
139 93 50 — 25 —100 —138 — 84 — 23 31 
— 16 13 12 — 6 — 49 — 93 —124 —I121 — 81 
—-m —-6--—> 1 — 16 —27 — 28 — 34 — 76 —103 
— 6 =< - 7 — $5 —.2 5 37 56 33 
— 58 —125 —156 —144 —116 — 84 60 150 180 
62 — 12 —135 —-211 —237 —228 — 65 38 78 
176 129 3 —104 —215 —276 —276 -—200 — 72 
125 160 97 — I —100 —169 —214 ~-218 —182 
— 82 15 63 113 96 58 10 — 66 —160 
—283 —252 —181 — 74 21 83 128 131 35 
—237. —324 —348 -—280 -—186 — 72 74 176 230 
— 18 —182 —336 -—460 -—406  —286 -—100 102 238 
225 we =I —33 “oe -—e ee lh Ke eS ee 
237 204 se -@ *-S.. + -. -  - 
95 190 200 140 40 —100 -—210 —285 —290 


Astron. Nach. 3333 and 3633. 
Resultate des Internationalen Breitendienstes Band III, pp. 223-224. 


Resultate.... Band V, pp. 180-182. 
Ergebnisse . . . Band VI, p. 222. 
Results...... Band VIII, p. 232. 


Definitive results reported to Dublin meeting of the LA.U. 1955. 
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(These results are not definitive.) 4 
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Nicolini, it appears that the observers are satisfied that the overlap in the 
programmes has been sufficient to minimize any inhomogeneity of the results 
because of these changes. 

Taste If 


Stations used in International Programme 


i Station Longitude (A) Dates of functioning 

1. Mizusawa sax” & E. 1899 Dec. to date. 

2.  Tschardjui (i) 63° 29’ E. 1899 Dec. to 1909 July.* 
Tschardjui (ii) 63° 35’ E. 1909 July to 1919 May. 
Kitab (iii) 66° 53’ E. 1930 Nov. to date.t 

4. Carloforte 8° 19’ E. 1899 Dec. to 1943. 

1946 to date. 
4. Gaithersburg 97° 12' W. 1899 Dec. to 1914 Dec. 
1932 Apr. to 1953 Aug.t 
5. Cincinnati 84° 25’ W. 1899 Dec. to 1915 Dec. 
6 Ukiah 123° 13’ W. 1899 Dec. to date. 


* The site of the observatory was moved in 1909 July. 
+ Kimura did not include the results of these stations when observations recommenced 
there, but Carnera did so after 1935. 


However, the star catalogue used from 1899 to 1935 was that of Cohn in 
which some of the declinations are in error. Since 1935, the more accurate 
Boss’s catalogue has been used but, even in this, errors large enough to affect the 
accuracy of the results are known to exist. 

The reliability of the recorded results might be affected by the fact that 
equal weight cannot really be given to each station. Because of their favourable 
climates, Carloforte and Ukiah usually contribute many more observations than 
the other stations. Between 1900 and 1924 from two to three thousand pairs 
of stars were observed at Carloforte every year, but in 1925, 1926 and 1927 the 
numbers of pairs observed were respectively only 1248, 928 and 1358. At that 
particular time only three stations, Mizusawa, Carloforte and Ukiah, were in 
operation, so the results for this period might well be less reliable than those 
for other periods. ‘There have been other intervals when there have been 
relatively few observations from particular stations. 

There is, finally, the “closing error” of the observations (see Melchior, 3, 
p. 17). This has always been systematically negative and positive respectively 
at stations in the Northern and Southern Hemispheres but, following the 
modifications in the observing procedure made by Kimura in 1922-7, the error 
was reduced in magnitude. 

We see, therefore, that the series of observations really cover a fairly large 
number of sub-intervals of time, each of which differs from the others in one or 
more of the number of stations in action, the star-programme used, the method 
of observing used and the method of reduction followed. 

We might expect that the sub-intervals are sufficient in number for our 
statistical theory to remain approximately true despite this inhomogeneity, 
systematic irregularities arising in a sub-interval due to the organization of the 
programme persisting only for a short time and changing sufficiently frequently 
for us to neglect the resulting bias in the estimates of the harmonic components 
comprising the assumed systematic part of the variation (see 1, equation (18)). 
In (1, p. 449), however, we suggested that the adequacy of our representation 
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of the systematic variation might be tested by dividing the series of observations 
into several consecutive sub-series and obtaining estimates from these which 
could be compared with estimates for the whole series. This no longer seems 
reasonable. If, on the one hand, we choose several short sub-series it seems 
probable that they would not include enough of the sub-intervals described 
above for our theory to be applicable. On the other hand, the inclusion of 
enough sub-intervals would give too few sub-series for useful comparisons to 
be made. 

4. The computational procedure follows that described in our previous 
paper. There we discussed, at some length, the factors governing the choice 
of the lags which ought to be adopted in evaluating the covariances of the various 
series and argued (pp. 450-451) that it seems best in general to use the smallest 
possible lags. We have, however, made a number of analyses using higher lags; 
for these we need formulae which are slight generalizations of those given before. 
In (1, p. 450) equation (29) generalizes to 

cot ry*h=p(r)/q(r) (3) 
where p(r) and q(r) are as given in equation (30), r being the lag, 27/y the period 
and A the interval between observations. ‘The damping factor, x, is estimated 
from 





e2k* k(r-1) — tpl) }? + ig)? (4) 
(P(r)? + ta(r)}* 
The equation 
tP(rs)}* + tq(r2)}* 
(P(r)? + igra)s 
may also be used, but detailed discussion of this is given later when we examine 
the question of bias (Sections 15 to 18).* 

5. We have analysed the series and sub-series which are listed in Table III. 
They include unsmoothed and smoothed series in which the z-term has been 
taken into account and unsmoothed series in which the s-term has been ignored. 
The series have been given code numbers for ease of reference. The first digit 
of the code denotes the type of series (e.g. unsmoothed series in which the 3-term 
has been taken into account are type-5 series), and the second digits distinguish 
the series according to the interval they cover. 

As one of the major purposes of the work is to determine the period and 
damping factor of the free motion of the pole, values of these constants are given 
in the table, but it should be noted that the values given are those which have 
been obtained by using the smallest lags. 

For the sake of completeness we have added the series which we considered 
in our previous paper, together with those considered by Jeffreys (2). 

6. First we considered in more detail the data which we originally examined 
in (1). In the reduction of these observations, the z-term was ignored, and as 
it seems certain that this term is physically real we did not think it necessary 
to try to locate data for the years after 1934 (if indeed such data are to be found). 
The thirty-five year interval is covered by series 31. Series 32 cover the 
twenty-two complete years, 1g00 Jan. to 1921 Dec. when the programme was 
directed at Potsdam and series 33 cover the twelve complete years 1923 Jan. 

* An asterisked quantity is the estimated value of the unasterisked one, except that /* and m* 
are estimates of /—/s (=x) and m—m,(=vy). In the tables asterisks are dropped because all the 
tabular values are estimates. 


oK* oo 
e2’* hiri—r2) — 


(4) 
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Tasze III 
Description of data and principal results 
ai Free period Damping factor 

Code Length Description (27/y) in years (x) in years~ 
(ole) 1890°0-1955°0 Smoothed; o-1 year 1°304+ 0°024 0°506+ 0'095 
oI 1900°0-1934'9 interval ; z-term taken 1'287+ 0°028 0°361+0°107 
02 1900‘0-1920°9 into account. Data of 1'249+ 0°025 0'198+ 0°100 

Table Ib. 
15 1892°0-1932°9 Smoothed; o-1 year 1'223+0°019 0058+ 0°015 
16 1908°3-1921°5 interval; z-term taken 1'202+ 0016 0:066 + 0°007 

into account. Data 

reduced by taking non- 

overlapping means of 

three consecutive values. 

Jeffreys’s data. 
31 1900 Jan. -1934 Dec. Unsmoothed; monthly 1'267 + 0'039 0°322+ 0°103 

2 1900 Jan.-1921 Dec. interval; z-term not 1'248+ 0'034 0°261+0°113 

33 1923 Jan. -1934 Dec. taken into account. 1°355+0°099 0°661 + 0°340 
34. 1900 Jan. -1905 Dec. 1°241+ 0°054 0°366+ 0221 
50 1899 Dec.-1954 Dec. Unsmoothed; monthly 1:287 + 0°026 0°441 + 0°095 
51 1900 Jan. —1934 Dec. interval; z-term taken 1'267+ 0°041 0°356+0°109 
52 1900 Jan. —1920 Dec. into account. Data of 1°238+ 0°033 o'239+ O11! 

Table Ia. 


to 1934 Dec., during which Kimura was director. The change in the method 
of making the observations was made in 1922 July, so this year has been omitted. 
Series 33 should be homogeneous in every respect, but in 32, although the 
methods of reduction and observation remained unchanged, there were changes 
in the star-programmes and numbers of stations in operation. We have also 
examined series 34 which cover the first six years of the work of the International 
Service and should be homogeneous. ‘The main results are given in Table IV. 
Series 33 appear to give results which are not in accord with the other three 
series of this group although it is unlikely that the difference would be highly 
significant. The physical explanation of this apparent anomaly is possibly to 
be found in the fact that the interval covered was short and seriously affected 
by the paucity of observations to which we have referred earlier. Series 34 give 
a value of the free period in reasonable agreement with the remaining two but 
the value of x* obtained from these series appears to be subject to considerable 
uncertainty. In this case, the fact that A* exceeds unity by more than three times 
its estimated standard error may well be attributed to the inapplicability of the 
asymptotic theory (developed in Sections 4 and 5 of 1) to such short series. 
(It will be recalled that A is the parameter that measures the effect of observational 


" error.) 


7. Next we examined the smoothed series in which account was taken of 
the z-term. We have analysed the longest available series, 00, covering the 
sixty-five years from 1890-0 to 1954-9 and also the series or and o2 which 
correspond to the series 31 and 32. The results for the type-o series are given 
in Table V. 

In this group of results there is a striking decrease of both the estimated 
values of the free period and the damping factor as the series shorten. ‘The 


results are not particularly consistent with each other, especially in the case of 
values of x*. 
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The most serious feature of these results is that A*>1. ‘Taking A=A* has 
little effect on the estimate of var«*h, but makes a great deal of difference to those 
of varA* and vary*h, the latter being in fact useless since it is then negative. 
The calculation of estimated variances is particularly suspect for this group 
of smoothed series in view of the values of the ratio of A*—1 to o/{(varA*), -1}: 
These are 12-6 for series 00, 10-4 for series 01, and 10°3 for series 02, all of which 
are so highly significant as to make the results unacceptable. ‘The discrepancy 
must of course be due to the behaviour of the auto-covariances and cross- 
covariances with lags 0, 1, 2 since only these are used in calculating the ratio. 

Since the three series overlap, the estimates of yh, « and o are related to each 
other; because of this it is very troublesome to carry out even approximate 
large sample tests for the significance of the differences between the estimates. 
The extent of the correlation could be worked out but the calculations would 
be so tedious that we do not think they would be justified, in view of the 
behaviour of A*. 


TABLE V 
Results for series of type-o 


(a) Harmonic components of systematic annual motion (unit 0”-001) 


Series oo Series o1 Series 02 

l m l m l m 
c +17°81 —28°10 + 19'00 — 20°92 27°40 ~— 22°33 
ay — 36°51 +75°25 — 39°76 -+-'79°87 41°12 + 66-92 
A, — 89-30 — 28-97 —77°36 — 30°69 — 64°80 — 31°24 
a2 — 090 — 0°36 + 1°05 — 2°93 - 9°26 — 2°43 
hy + 1°32 — 0°70 + 1°68 0°25 + 2°19 + 0:06 
a; + o61 +- 0°36 + 2°62 — 90:96 t+ 1°80 — 1°03 
hb, — 1°93 + 078 — 1°50 t+ 1°42 — 0°08 — ol! 
a, — 1°58 — 0°49 — IIs — 0°04 0°00 — 1°35 
by + 0°06 + 0°55 — o'21 +- 0°53 + 0°22 - o-61 
as — 0°69 — 0°03 + oO'or — 0°23 + 0°45 — 0°59 

(5) Free Motion (calculated with least possible values of lag) 

yh 27° 37'+ 31" 27° 58'+ 36’ 28° 50'+ 34! 
period 1°304+ 0°24 years 1'287+ 0°028 years 1'249+ 0°025 years 
K 0°506+ 0°095 year=! 0°361+0°107 year~! 0198+ 0'100 year™! 
A 1037 + 0°003 1°029+ 0'003 1°020+ 0'002 
Go 37°3 28°6 25°0 
Co eee 


N.B. The standard errors were calculated using A=1. 

8. Finally, we examined the unsmoothed series in which account was taken 
of the s-term. The longest suitable series, 50, cover the fifty-five years 1900 Jan. 
to 1954 Dec. Series 51 and 52 correspond to series 31 and 32 respectively. 
The results derived from this group of type-5 series are given in Table VI. 
Here we see again that the estimates of the length of the free period and of 
the damping factor decrease with the length of the interval. This group of 
results shows much more consistency than do either the groups of type-o or 
type-3 series and there is no trouble over the values of A* which are all significantly 
less than unity. 

Again the overlapping of the intervals makes comparisons difficult but the 
estimates of the standard errors are now acceptable and so do give a rough idea of 
the consistency of the estimates. We can use the inequality o(X — Y) <o(X) +.o(Y) 
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(for any random variables X, Y) to obtain an under-estimate of the significance 
of a difference—though admittedly this is very much of an under-estimate if 
there is a fair degree of positive correlation between the two variables. Thus 
we would probably be justified in concluding that the difference between the A*’s 
for series 50 and 51 was significant since this is 0-0182 compared with an estimated 
standard error of at most 0:0022 + 00038 = 0:0060. 


Tasie VI 
Results for series of type-5 


(a) Harmonic components of systematic annual motion (unit 0”:oo1) 


Series 50 Series 51 Series 52 
est. est. est. 
l m stand. / m stand. / m stand. 
error error error 
¢ + 18°80 —31°72 + 3°17. +12°92 —22°02 + 3:20 +19°00 —23:30 + 3°86 
_ . +69: = . 1 . oe . ak. . 
ay 64°45 +69°51 + 10°98 55°33 +74°99 4 11.46 48°42 +65°67 + 15°75 
by, —71'4l —45°59 -69°59 —45°89 — 59°68 —37°05 
a — 2a — £96 . ~ 2°55 — 6-6o — 2°51 — 7°95 
. + 2a + 2-5 + 3°00 
by + 5°66 — 0°39 ~ 99 - 305 * *° | 837 — 545 > 3 
as — o'16 + 1°38 1 1°84 2°36 aid. + 2°91 + 2°29 
= + 3 1°96 
bs — 259+ 279> "44 — gets 393 ~ — 119 + 4387 *Y 
ay — 0-78 + 1°47 . oz + O65 | — 0°65 + 089 . 
b, + 1°62 + 1°§2 adi, 125 + 1°36 ~ | 53 + o81 + 4°28 — vines 
as. + o-21 + 0°97 + 508 0°58 + o-4!1 soe 1°85 + 1°85 +e 
b, + 1°59 + 0°63 - aa + 597 + 2°11 + 2°94 
a6 — o71 — 0°85 + 0°70 — 0°72 — 1:63 + 1°00 + O19 — 2°65 + 1°03 
(b) Free Motion (calculated with least possible values of lag) 
yh 23° 18'+ 28’ 23° 41' +46’ 24°14' +39 
period 1°287 +0-026 years 1°267 +0°041 years 1'238 +0°033 years 
K 0°441 +0°095 year! 0°356 +0°109 year"! 0'239 +o-'111 year~! 
A 0°9975 + 070022 0°9793 + 0°0038 0°9942 + 0'0024 
te 32°4 25°8 25°2 
a’ 6:0 15°6 9'8 


, 


g. The components of the systematic annual motion shown in Tables IV, V 
and VI vary according to type of series and interval, there being noticeable 
differences in both the amplitudes and phases of the several components. When 
the tests of significance used in (1) were applied, we found that on the whole 
only the annual and semi-annual components are highly significant. When 
we have actually evaluated the estimated standard errors of the components, 
they appear to be relatively large. ‘The differences of corresponding harmonic 
components of type-5 series are given in Table VII, together with the estimated 
standard errors of the pairs of components. ‘The overlapping of the three series 
has again the effect of making any precise assessment of the significance of the 
difference very troublesome, but in view of the magnitudes of the estimated 
standard errors no startling discrepancy is indicated, although the differences 
50-51 and 50-52 for the constant terms in the m components are suggestive. 

10. There may be systematic trends in the observational data caused by 
genuine secular trends in the systematic terms in F, and G;, (cf. 1, equation (3)). 
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We have examined the series of type-5 for such trends by fitting curves of the 
type 
1, = Cy +¢,t+acos vt+bsin vt (27/v=1 year) 

to each interval covered by only one director or only one star-programme. The 
results give values of cy markedly different from those of ¢ recorded in Table VI 
and the values of c, are so large that the removal of systematic trends of the 
kind assumed here yields series for the free motion which are visibly irregular, 
with large discontinuities at each change in star-programme. If the trends are 
real, we cannot explain them satisfactorily. However, the frequent changes in 
programme make it difficult to determine trends accurately, and the large 
discontinuities found at each change in star-programme are likely to be due 
more to accidental rather than true physical causes. 

11. We have examined the residual motion {e’*, 7'*} which includes the 
effects of both observational errors and the random disturbances which maintain 
the free motion. We have restricted our attention to those of type-5 series. 
The residuals have been estimated in the way described in (1, p. 458) for each 
of the series 50, 51 and 52, those for 50 being displayed in Fig. 1 (0, d). 

If the series were heterogeneous and corrections such as those given by (1) 
should be applied, then the numerical harmonic analysis would leave in /*, m* 
errors which remain constant through each homogeneous sub-interval, but 
which change between such sub-intervals. As a result {e’*} and {n’*} would 
also be affected by similar errors. We have therefore determined the means of 
fe’*} and {n’*} over various sub-intervals during which star-programme and 


Tasie VIII 


Mean residuals : series of type-5 (unit 0”:001) 


Series 50 
Star Programme Interval é 7’ 
I 1900-05 2°89+ 2°94 6-89+ 2-90 
II 1906-11 — 6:00+ 3°22 5°39+ 3°03 
Ill 1911-22 July ' 8:83+ 2-60 — 7:00+ 3°38 
I-III 1900-22 July 3°3It 1°71 — 0°02+ 1°97 
IV 1922 July—34 2 67+ 2°60 9°89+ 2°67 
Vv 1935-48 3°73 + 2°30 — §°24+ 2°24 
VI 1949-54 — 24°15 + 5°44 - 986+ 4°42 
I-VI 1900-54 0°28+ 1°28 — o19+1°28 
Series 51 
Star Programme é’ 7 
I — 0°06+ 2°89 3°58+ 2°93 
II - $33 23%9 221+ 2°94 
Ill 590+ 2°59 — 10°43 + 3°35 
I-III 0°33+ 1°70 — 3°354 1°95 
IV — 5742-71 6-17+ 2°65 
I-IV — 0°35+1°46 0°03 + 1°59 
Series 52 
Star Programme € 7 
I — 026+ 2°86 6°35 + 2°93 
II — 9°57+ 3°07 508+ 2°88 
III 5°44+ 2°80 — 660+ 3°41 
I-III — 048+ 1°73 0°44+ 1°91 
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direction remained uniform. The results are given in Table VIII. The 
uncertainties given in the table are the standard errors calculated in the normal 
way but, because of correlations introduced by the processes by which the 
estimated annual and forced motions were removed, and possibly also correlations 
between the true residuals ¢’,, and e’,,,, and 7’,, and 7’,,,,, they are not unbiased 
estimates of the true standard errors. The bias is very troublesome to calculate 
but may be expected in all cases to be positive, so that they are useful as estimates 
of upper limits to the uncertainties. The results suggest that the data do suffer 
from inhomogeneity. 

12. When the results are grouped according to the intervals covered by the 
series it is immediately evident that the results for the group 02, 32 and 52 are 
very consistent, and those for the group o1, 31 and 51 are slightly less so. ‘The 
agreement for the longer series 00 and 50 appears to be less good. 

The evidence of the statistical analysis so far presented does not give any 
conclusive reasons for preferring or rejecting any particular set of results, although 
these for type-o series are suspected to be unacceptable because A*>1. When 
the results are considered in conjunction with the known facts about the initial 
data it seems reasonable to draw certain inferences. 

The data might be inhomogeneous, and if the corrections given in (1) are 
applicable in full the long series might be expected to yield unreliable results. 
The results for series oo and 50 seem to support this, and indeed, the results 
derived from series 50 appear to be less untypical compared with the main body 
of results than do those for series 00 which differ from 50 by the inclusion of 
the results of the uncoordinated observations made in the decade preceding the 
beginning of the international programme. 

It might also be expected that the series covering the first twenty-one years 
of the international programme would give better results than those covering 
the first thirty-five years. We have already recorded that the differences in the 
constant terms of the harmonic analyses of the m components between the 
series 50 and 51 and between the series 50 and 52 are suggestive, and although 
these differences are not what they should be if they were due to the non- 
application of corrections of the order given by (1) they still support the opinion 
that the data which we have analysed are inhomogeneous to some extent. 

It also appears that the variability of the results depends more on the interval 
covered by the series than on the type of series. The fact that A*>1 in the case 
of type-o series does suggest that smoothing has had an effect, but the differences 
between the results of type-3 and type-5 series suggest that the inclusion or 
exclusion of the z-term has had little influence. On observational grounds, 
however, the s-term must be accounted for and we conclude that the most reliable 
results are those derived from type-5 series. 

The results of the analysis of the short series 33 and 34 give periods and 
damping factors which seem unsatisfactory, but their estimated standard errors 

_are large enough to account for the discrepancies. The low precision of the 
estimates and the probable inaccuracy of the asymptotic theory for short series 
are likely to be general and we have not made any other analyses of such short 
series. This has led to the disadvantage, already noted, that the intervals are 
overlapping, but that has not proved to be too serious a defect. 

13. Jeffreys’s results are given in Table III as series 15 and 16. He used 
essentially the same data as those of our type-o series. The periods which he 
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obtained are in fair agreement with those for our series 01 and 02 but the damping 
factors are markedly different. His value of A for series 15 was less than unity. 
However his method of computing « and A was different from ours, being 
equivalent to solving (in our notation) 


N2e-2khr — {p(r)}? + {a(r)}*: 
{P(0)}? 
by least squares using three values of the lag r. He advocated the use of fairly 
large values of r, whereas we have argued that the smallest possible lags should 
give greatest accuracy (1, p. 450). Jeffreys also reduced the volume of data by 
taking non-overlapping means of three consecutive observations so that his 
lags 0, 1 and 2 correspond roughly to our o, 3 and 6. 

In order to make a direct comparison with Jeffreys’s work we have examined 
two series, 05 and 06, for the intervals 1892-0 to 1932-9 and 1908-3 to 1921°5 
which were the intervals of Jeffreys’s series 15 and 16 in Table III. The data 
differ only in that we have not grouped them. We have calculated covariances 
only for every third lag and have used equation (4) in the modified form 











ee *h(r—-1) — {p(3)}?? ia {9(3)}? (5) 
tP(3r)}? + {9(37) 3” ; 
and (44a) with r, = 3r and r,=3r—3, so that 
eoxen — (P(37—3)}* + t9(37— 3)". (5.4) 


(P(37)}* + 19(37)}° 

The results given in Table LX show that for small lags, the estimates of the 
period agree favourably with those given by Jeffreys so it seems that his grouping 
of the data has had no serious effect. ‘The damping factors, except those for 
series 06 computed by formula (5), seem to be irregular, but for both period and 
damping factor it can be said that as the lag increases there is initially a marked 
decrease in the estimates followed by a more or less marked oscillation. 


TasLe IX 


Analysis of series 05 and 06 for various lags 


Period in years Damping factor in year~! 
3r Series 05 Series 06 Series 05 Series 06 
Formula Formula Formula Formula 
(5) (5a) (5) (5a) 

3 1°230 1'216 

6 1°193 1‘206 +0169 +o*169 +0°064 +0°064 

9 1°178 1*202 +0°037 —o0'096 +0°028 —o'008 
12 I*190 1°201 ~O°005 —o'088 +0024 +0°017 
15 1*209 1‘201 +0036 +Or159 +0022 +0°015 
18 1°189 1°203 +-0:060 +0°159 +0°022 +0021 
21 1°184 1°202 +0°035 —0*0go +o°022 +0°022 
24 1°187 5*200 +0021 —0'063 +0021 +0°012 


14. As we believe that the data of the type-5 series are superior for the purposes 

of analysis, we have carried out more calculations for the series 50, 51 and 52. 

The results are given in Table X. (Formula (4a) was used with 1, =7, 72=1r—1.) 
1Q 
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The estimated period is clearly seen to be a damped oscillatory function of 
the lag. For the damping, the estimates follow a similar pattern but there is an 
appreciable difference between the results according to which formula is used. 


TABLE X 
Analysis of series 50, 51 and 52 for various lags. 
Period in years Damping factor in year~' 
Series Series Series Series 50 Series 51 Series 52 

Ir 50 51 52 Form. 4 Form. 4a Form. 4 Form. 4a Form. 4 Form. 4a 
I 1'288 = 1268 1°237 

2 1°265 1°253 1'227. 0°443 +0°443 0°352 +0°352 0°231 +0°231 
3 1'241 1°239 1°219 0°457 +0°47I1 0°379 +0°406 O°241 +0°251 
4 1°220 1°226 I‘212 0°436 +0°394 0°389 +0°409 0'234. +0°220 
5 1°205 1'215 1'207. 0400 +0°291 0°387 +0°380 0:°223 +0°190 
6 I‘I194 1°203 1'203 0°348 +0°142 0°350 +0°200 o'199 +0°105 
9 1185 1°193 I*199 0°315 +0°I50 O°315 +0°I40 0°174 +0°045 
8 1178 =. 1186S sa"196. 0270 —0°003--0'254 —0'107. = 0138 —0°074 
9 1172s '181 1'194. 0'221 —o'rr8 0°200 —0'184 0-105 —0'127 
10 1°169 1179s: 193. -0"166. —0'276 = O51 —0'239 §=— 00077: —0'152 
II 1169 ~=.:11'178 I°192 O17 —0°318 =—0"109 — 0263-0051 —0'180 
12 I‘r7I 1180 1°193 «©0083. —0'257 90077 —0°243 «= (0034 — 0138 
13 1177 1183 1°194 0066 —0'128 0°066 —0'057 0°033 +0°024 
14 1°182 1°185 I'I95 0°070 +0'118 0069 +0°097 0°039 +0°114 
15 1°185 1187 = r'19g5 = 0080 +0°216 §=0°075 +0°165 07045 +0°126 
16 1°187 1189 = I'IQS5) «= ©092_- + 0'254 «= «081 ~+0°158 = 050s +. 0108 


15. Professor Jeffreys has suggested to us that the use of higher lags may be 
advisable because of the possibility that the disturbances « and 7 are not 
uncorrelated as we assume in our model. There are good physical reasons for 
believing that correlations between the disturbances do exist; e.g. the seasonal 
climatic changes which give rise, at least in part, to the systematic terms F(t) 
and G,(t) (1, equation (4)) are not entirely smooth, and irregularities in them 
tend to persist over intervals of several weeks, so that there may well be appreciable 
positive correlation between values of € or 7 at times separated by intervals up 
to 2 or 3 months. 

Such correlation may cause our estimates of the period and of the damping 
factor to be biased. We have investigated the bias on the assumption that 
COV (€), €/.¢n) =COV(M Hira) =r SAY, While cov(€,, 7,,,) remains equal to zero 
for all r. It is then easily shown that 





b(ry*h) = tan“ ae } —ryh (6) 
and for the estimate «* given by equation (4a) 
*{A(r,) +1 B(re 
b{(r, —172)e*h} = log as Wee (7) 


where 
Alp) +iBir)= _& ertarovineny, .., 


8,t=0 


and we use the notation 5(6*) for the bias in an estimate 6* of a quantity 0. 
In particular, when g,=gou" (r>o), where |u]<1, i.e. the auto-correlation 
function of the series {¢,}, {nm} is that of a first-order autoregressive process, 
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(6) and (7) give 








«}\ 7 2eKhsin yh achyt\ fe 
b(ry*h) = tpi apeueph a =a sin {(r + 1)yh—} (8) 
and 
h 
Bi(r, —1)K*h} = =a [a coe (71 + 1)yh — $} —p* cos {(r, + 1)yh—4}) (9) 
where 


$=atan-t(— Ler. ) 


1—pcosyh 

and it is assumed that r, 7,, r. are sufficiently small for e”**, e"*”, ev” to be 
taken approximately equal to unity. ‘Thus as 7 increases the bias in y* tends 
to zero like 1/r after an initial damped oscillation. When «* is calculated by 
formula (4) the bias is given by (g) with r,=r and rs=1. This behaves in the 
same way as the bias in y*. When «* is calculated by formula (4a), however, 
the bias is given by (9) with r,=r and r,=r—1 and oscillates with an amplitude 
that decreases like uw’. ‘The behaviour of the biases should in general not be 
very different from this since gp, g,, g2, ... may be expected to be a decreasing 
sequence of positive numbers. In fact, if g, can be put equal to zero for r>g, 
it follows from (6) and (7) that for r>g, the biases of y* and «* calculated from 
equation (4) are proportional to 1/r, while the bias in «* calculated from 
equation (4a) with r;=1r, rg=r—1 is negligibly small. ‘This should be the case 
for q=3 or 4 since correlations between the disturbances should not extend over 
intervals exceeding a few months. 

We conclude that correlation of this type would suffice to explain the 
behaviour of the estimates of T=27/y and « in Table X for small values of r. 
However, it cannot reasonably be considered responsible for the persistence 
of oscillations for the larger values of r (those of «* in the last three columns 
being particularly marked), since that would require the auto-correlations g,/g, 
to remain appreciable for these values of r. 

16. The imperfect removal of systematic variations from the residual series 
used for the calculation of p(r) and g(r) will also produce bias in our estimates. 
Suppose that {é,, ¢,} represent systematic variations which have not been allowed 
for in determining the free motion. Then p(r) and g(r) are calculated from 
{X,, Y,} instead of {x,, y,} where 

X,= + & Y;=¥+ % 
If this is so, it can easily be shown that the consequential biases in p(r) and q(r) 
are 


b( p(r)) = cs 2 (Evesent Soerrn)s (10) 


learn — Sternde)s (11) 





N being the total number of observations. 

The results of Section g suggest that the annual term is different in different 
sub-intervals so that the process which we adopted of removing a uniform annual 
term throughout the whole interval must lead to the presence of terms &, and € 
which take the form 

£,= a) + a, cos(vt +¢), 
f,=b) +6, cos (vt +), 
10* 
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where dy, 4, 5, 5;, ¢ and % are constants and 27/v=1 year; there may also be 
small harmonic terms. In the true situation, the constants change from one 
sub-interval to another, but for simplicity we will first consider them unchanged 
over the whole interval. We find that the biases in p(r) and q(r) lead to 
—e’*" sin ryh 
2V,(0) 
provided that A is nearly equal to unity. (V,,(0) is the expected value of }p(o) 
in the absence of bias—see (1, equation (14)). If the damping is estimated by 
using equation (4), the bias is given by 
b{2(r —1)x*h} = {cos yh[a,? + by? + $(a,? + 6,”) cos vh] 

— eh cos ryh[ay? + by? + 3(a,? + b,?) cos vrh}}/Vy,(0) (13) 

and if it is obtained by using equation (4a) with r,=r,+1=r, it is given by 


(r—1)xh 
b{2x*h} = ‘ ) {(a? + by?) [cos (r — 1)yh —cos ryh] 
1 


b(ry*h) = (a,2 + by? + 4(a,? + b,?) cos vrh) (12) 


V 

+ }(a,? +,*)[cos (r —1)yh cos o(r — 1)h—cosryhcos vrh]}. (13) 
In each of these formulae we have omitted ‘“cross-product” terms involving 
a,b, sin(¢—) which appear in the detailed working. Their inclusion would, 
at worst, double the estimated bias, but should have a much smaller effect. 

If the period 2z/y is denoted by 7, then 

B(T*) _ _ (ry*h) 
a, —_" 
and the bias in 7, regarded as a function of r, is an oscillating one with an 
amplitude which varies with e™"/r, which decreases as r increases, at least 
initially. 

For the damping, formula (4) leads to a bias in «* which consists partly of 
a term decreasing with 1/(r—1) and partly of an oscillation with an amplitude 
which decreases (initially) with e”-!*"/(r—1). Formula (4a) leads to an 
oscillatory bias whose amplitude varies with e*-*", 

The formulae given are not valid for very large lags, when the second order 
terms, which have been neglected, may become more important than the terms 
given. 

When this analysis is extended to deal with €, and ¢, whose amplitudes change 
from sub-interval to sub-interval, the results are similar but, because of the 
increasing overlap of sub-intervals as r increases, the biases tend to decrease 
more and for sufficiently long series with sufficiently many sub-intervals the 
biases could become negligible—thus giving a justification for a remark made 
to this effect in Section 3 above. 

17. Values of ay, bo, a, and 6, consistent with the corrections given by (1) 
yield values of b(p(r)) which are comparable with p(r) itself, and our results do 
not indicate such serious effects. However, the results of harmonic analysis 
given in Tables IV, V and VI show that the inhomogeneity could give values 
of the constants of order of magnitude 10 units (the unit being o”-oo1). Putting 
a) =b)=a,=b,=10 units, the largest b(p(r)) is about 50 units and hence the 
largest bias in ry*h is about 0-02 in the range of values of r up tor=16. The 
results correspond to a bias in the period which is approximately 0-06/r years, 
which is in rough agreement with the behaviour of our estimates of T. Again, 
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using formula (4), the bias in «* given by the above analysis is at most about 
0°36/(r—1) which also compares with the behaviour of the estimates of x. 

18. It can also be shown that other kinds of systematic inhomogeneity in 
the data, not necessarily directly connected with the amplitudes of the annual 
terms, lead to biases similar to that arising from the sinusoidal terms (&, ¢)). 
For example, systematic errors in individual star-places could give rise to a 
systematic contribution to the observational error which would repeat annually 
until the particular stars were replaced at a subsequent change in the observing 
programme; such an error would contribute to the bias in much the same way 
as considered above. 

We may conclude that the behaviour of our estimates is explained by biases 
introduced by inhomogeneity and auto-correlation of the disturbing function. 

19. We believe that the sub-intervals during which the data can be assumed 
to be homogeneous are too short ever to give a reasonable prospect of eliminating 
the systematic motion adequately ; this is really a problem which is better tackled 
in the design of the observing programme rather than in the subsequent analysis 
of the results. Despite the complexities introduced by the inhomogeneity and 
the auto-correlation of the disturbances, however, it seems practicable to make 
a reasonable estimate of the period. First we have to detect a lag beyond which 
it is unlikely that auto-correlation of the disturbances is effective. Secondly, 
considering the form of b(ry*h) in equation (12), we must judge where the 
oscillations have decreased to an extent such that a mean value is well determined. 
Inspection of the results in Table X suggests that both these requirements are 
met by taking r>6. The maximum variation in the estimated period thereafter 
is less than about 10 days for series 50 and 51 and 4 days for series 52. The means 
of the estimated values of T= 27/y for r=1 to 16 and r=6 to 16 are given below. 


Series 50 Series 51 Series 52 
r=1 to 16 1-198 years 1-202 years 1*201 years 
r=6 to 16 1-179 years 1-186 years I'193 years 


On observational grounds, series 52 appear more likely to be homogeneous 
than the others. Series 50 contain an interval for which the published results 
are not yet definitive, and both series 50 and 51 cover intervals before and 
after 1922-7 when a major change in the method of reducing the observations 
was made. ‘This view is supported by the fact that the estimates of the period 
made from series 52 show markedly less variation than the others. We therefore 
consider the best estimate of the period to be 

I*193 years = 435°8 mean solar days. 

Unfortunately there does not seem to be any acceptable way of estimating «. 
As far as freedom from bias is concerned, the most reliable estimate would be 
obtained from formula (4a), choosing r, to be sufficiently large to eliminate 
bias due to the auto-correlation between the disturbances {e,y}, and r,—7r, 
sufficiently large to eliminate the bias due to ignoring part of the systematic 
variation, On the other hand, to keep the standard error of «* small, r, and rz 
should be small, so that it becomes a question of balancing conflicting require- 
ments. To make a theoretical estimate of the best combination of values of r, 
and rz is, in our opinion, impossible in view of all the unknown factors. 

The results in Table X show that an oscillation persists in the estimates of « 
which is such as to make it very difficult (even bearing in mind the probable 
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form of the bias) to decide which value to adopt. For series 52, formula (4) 
gives small estimated values, but it is hardly justifiable to say more than that, 
in the absence of bias, there is some evidence that the relaxation time lies in 
the range 10 to 30 years. To attach a measure of probability to this statement 
seems to be impossible. 

20. None of the results will be acceptable if the model should be proved to 
be completely invalid, and it seems fair to point out that Melchior (5) and 
Gutenberg (6) have both expressed doubts over the question of the damping of 
the free motion. The case for the model was, however, very well stated by 
Jeffreys in (2); our model is a generalization of the one used by Jeffreys and is 
designed to take further account of observational error. As we showed in (1), 
it contains Jeffreys’s model as a special case, so that we have confidence in the 
model at least from this aspect. 

Jeffreys, however, isolated what we have called series 16 for special attention 
because they covered an apparently disturbance-free interval. We do not wish 
to add to the comments we have made on this isolation (1, p. 458) so far as the 
numerical analysis is concerned, but it should be understood that, in establishing 
the validity of our model, as well as that of Jeffreys, we used the assumption 
that there is a fairly large number of disturbances in the interval between 
successive observations (apart from the limiting case of completely disturbance- 
free motion, when its validity can be seen by putting f,(¢) and g,(t) and consequently 
e and 7 all equal to zero in (1, equations (4) and (13)). The model will certainly 
not be valid if there are, as Jeffreys thought, quiescent periods and highly 
disturbed periods, except perhaps if the observations were extended over a 
very large number of years indeed. ‘This is important in the light of current 
geophysical theories of the Earth’s core (see e.g. Revelle and Munk, 7). If the 
compensations envisaged in such theories at the interface between the crust and 
core occur only at infrequent intervals—and such seems to be believed—then 
doubt must arise over the validity of our model. 

21. We conclude: 

(a) The harmonic analysis of the data is unlikely to do more than reveal 
qualitatively that inhomogeneity exists in the data because changes in the 
programme have occurred too frequently to permit the statistical theory to apply 
over the short sub-intervals involved. 

(6) The unsmoothed data at monthly intervals (series of type-5) appear to 
be best, and of these the short series 52 are to be preferred. 

(c) Analysis of type-5 series using high lags indicates that the estimates are 
biased. Qualitatively, the bias is consistent with the hypotheses that the data 
are inhomogeneous and that auto-correlation exists in the series of disturbances. 

(2) The period of the free motion may be taken to be 


1°193 years = 435°8 mean solar days. 


(e) The bias is fatal for the accurate estimation of the damping factor for 
which, however, there is some evidence suggesting chat the relaxation time is 
between about 10 and 30 years. 
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THE THEORY OF NUTATION AND THE VARIATION OF 
LATITUDE 
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Summary 


The theory of the bodily tide and of the various nutations is developed. 
Elasticity of the shell and fluidity of the core are taken into account. ‘The 
model used for the shell is Takeuchi’s Model 2, based on one of Bullen’s. The 
core is replaced by a homogeneous incompressible fluid, with an additional 
particle at the centre chosen to make the mass and moment of inertia of the 
core correct. Effects of the ocean are neglected, but can be allowed for at a 
later stage. 

The period of the free nutation is found to be 392 days (which would be 
increased by the ocean). The bodily tide numbers for semidiurnal and 
long-period tides are : h=0-58, k=0-29, 10-082. They have other values 
for diurnal tides, the greatest observable differences from the statical values 
being for the lunar tide O. The correcting factor for the 19-yearly nutation 
is 0°9964. 





1. The treatment of periodic changes of the Earth’s rotation given in most 
books regards the Earth as a rigid body. It has been known since Newcomb that 
elasticity lengthens the free period; by itself its effects on the precession and 
nutations are negligible. On the other hand actual fluidity of the core shortens 
the free period, and in the simplified case of a rigid shell and a homogeneous 
incompressible core it produces a reduction of the amplitude of the 19-yearly 
nutation that is too great to agree with the facts. 

Allowance for elasticity of the shell has been found to reduce the effect of 
fluidity of the core on the nutation. The first treatment (H. Jeffreys, 1949) used 
a Wiechert—Herglotz model, with the shell and core homogeneous and incom- 
pressible. With this no more than a rough agreement was to be expected. 
Recent work on the velocities of elastic waves and on the density distribution has 
made a detailed discussion by numerical integration possible, and such a study 
has been carried out by H. Takeuchi(1950). His solutions are based on two of 
K. E. Bullen’s models (1936, 1940) but the models used actually differ in some 
minor respects from Bullen’s. Three simultaneous differential equations of the 
second order are solved by a Taylor series method, and the calculated values for 
the bodily tide numbers 4, k, / are in good agreement with observation. 

Takeuchi used a statical theory for both shell and core. It is known that 
this is valid for the shell for disturbances with the periods in question, and for 
the core for the semidiurnal, fortnightly, and semiannual tides. For distur- 
bances tending to alter the axis of rotation, and therefore for all diurnal tides and 
the variation of latitude, a hydrodynamical theory of the core is necessary. His 
work, however, greatly assists such a treatment because it gives the general statical 
theory of the shell in terms of six adjustable constants; thus the only change 
needed for the shell is in the expression of the boundary conditions. 

In a first attempt to use his solution the core was still treated as homogeneous 
and incompressible, but this was found to be unsatisfactory (H. Jeffreys, 1951b). 
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In the present paper, accordingly, this assumption is modified. ‘Ihe actual core 
is too complicated for detailed treatment, and simplified models are used. 

We use a Lagrangian specification as before. The axis of z is in a fixed 
direction, while those of x and y rotate about it with a constant angular velocity 
w. If x; are the coordinates of a particle in its mean position, €; in its actual 
position, where we shall exchange xj, x», x3, with x, y, z, and €,, &, 3, with , n, f 
as convenient, the kinetic energy will be 


bf ff pol(E—on)? + (9 + 8)? + Phar, (1) 
where py is the undisturbed density and the integral is through the undisturbed 
position. ‘There is actually a second-order variation of the angular velocity 
about 2, but its reaction on the displacements in x and y will be of the third order. 

The possibility of adapting the statical solution for the shell depends on the 
fact that if the displacements arose from a pure rotation £7 — &n would be of the 
form ¢ (£% —7y/) and the terms in w and w? could be represented by an addition 
to the potential function. ‘This requires care, however, since € and y are not 
small and second-order terms in €, 7 need attention. We suppose the particle at 
x; displaced to x; + u,’, and then to be carried to €; by a rotation, so that 

€,=%,(1— $27) +u} + U(x3 +5) — b/mx,, 

£,= x_(1 — $m?) + uy + m(x3 +5) — b/mx,, (2) 

&, =X (1 — $1? — $m?) + uy — U(x, + u}) — m(x2 + U5). 
The w* terms are equivalent to an addition $w?(é? + 7») to the potential function, 
and can therefore be included in the work function W. ‘To develop this we first 
consider a displacement u;. The gravitational potential consists of the un- 
disturbed part Up, a small external disturbance U,, and a part U, due to displace- 
ments of the body. Then in a small variation of u; the work done by gravity is 





9 
sw, = || J po ag (Uo + Us + Us)udr, (3) 
where the derivative must be evaluated at x;+u;. To the second order 
bie aU, 2U, aU, au, 
6W,= {lf Po (= _ Ox;OX, " Ox; : OX; a, (4) 


all derivatives being now taken at x;. 
Since u; 0U,/0x; is quadratic in the displacements this gives 





, . aU, . , aU, dU, , ,. dU, 
~— uj z— + hue, -— + uj 5— + hu, 5— ) dr. 
Wy IJ po( ' OX; * 2MiUe OX}, * Mi OX}, ™ = OX; r (5) 
The stress is 
— Podix + ANB in + 2HEjKs (6) 
where py is the initial pressure at x;, 
ou Ou; Ou, 
A= > m : 2e ? k (7) 





a* = + =. 
m Ox, Ox; 
The relative increase of volume of an element is, to the second order, 
ou; > CU, OU, Olly Os 


OX, 


O(x; , = 
II ¢ (x; + u;)/Ox, lI OX, OX OX; OX, 


) 


Ou; | , Ou; Ou, , Ou; OU, 
~ as; * 8a, Oey bay Ba ™ 
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Thus py does work in the nae equal to 
Ou; Ouy, Cu, CU; 
W,, = | Po (= +3 — -1—- — dr 
af ’ "Ox; OX; ~ OX; OX, 


ine ‘ ~ the Ft) 48 


[foam E(nB)-m(nB}e — 


The surface integral is zero. ‘The part over the outer surface vanishes because 
p)=o. Over the inner boundary fy is constant and is continuous. ‘To the first 
order the total volume of matter carried through dS is JjudS. ‘The second order 
terms might differ according as they are evaluated just inside or just outside the 
core. Consider the displacement imposed gradually ; then u; is the displacement 
of the particle that starts at x;, while the particle that finishes at x; starts, to the 
first order, at x,’ such that 


A 


F ou; 
x= xj + U(x;')= x +, +(x hz, = 


and its displacement is u;—u,0u,/dx,. Thus the mean displacement of the 
particles that cross dS is u;— }u,0u,/ox,. With an interchange of suffixes this 
gives the first and third terms of the surface integral. ‘The second term gives 
the increase of volume due to expansion of the matter transferred. Hence the 
integral is zero, since, apart from po, it represents the difference of the calculated 
volumes within the core boundary according as the calculation is based on core 
or shell values of the displacements, and we are supposing that the core and shell 
continue to fit everywhere. 


Since the initial stress is hydrostatic, 


op aU, 
Os, OO, (11) 


and the first order terms in W, and W, cancel. The second order terms give, 
in all, 
2U, aU, OU, , OU, du, |, dU, du, 
W= {ff po (du ke Tx Ox, a “5 +13 1 Se Set + bt Se) 


-3 ff | (y+ ae) Sa. (12) 


If we make small variations 5u,, apply Green’s lemma, and make some interchanges 
of visite suffixes, we ‘i the equations of equilibrium 


CU, dU, ou,  dU,du 
=— (AAS, ' or o . 2 s 
aa wt He) +69 (Ur + Uz) + po (ue Ox;OX, Ox, Ox; Ox, Lae x Ox, Ox: at) = 


(13) 


The terms containing U, are expressible as 


al C) aU, 
sa (rate) 5" ay + oy, (ous Ox iz): (14) 


and when we use the a 


(Po, Uo) 
Hx, 4) X) as (5) 
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these reduce to the terms depending on Xo, Yo, Zy in the equations (esp. (20)) 
previously found by the Eulerian method (Jeffreys, 1929, p. 163). 
We take first u;=u,’. In calculating W for u,’ we must replace U, by 
Y = U, + w(xi + x3). (16) 
We then impose a rotation as in (2). In this the internal reactions depending on 
the mutual attractions and internal stresses do no work and do not alter W. The 
initial pressure, as modified by the inclusion of the rotation terms in Y, still does 
no work, but the extra term added to U, does some. ‘The contribution is the 
difference between [ff py}w(&7+&)dr in the positions specified by x,;+u,’ and 
€; and hence is 
> }a* | | | p(&;— x; — u, (€; + x; + u;')dr. (17) 
i=1,2 PF 
Evaluating this to the second order and remembering that the axes are principal 
axes of the body in its standard position, we find that it is 
rr ° 
— Jw(C — A\(P +m?) + ( J | pot — xy(lx + mx)dr. (18) 
JJJ OX; 
The whole displacement has been separated into the displacement u;’ followed 
by a rotation, and so far it has been arbitrary how much of it is interpreted as 
belonging to u;. We can remove the ambiguity by introducing the condition 
that u;’ in the shell makes no contribution to the angular momenta; this is 
equivalent, to the first order, to 
fff PoE itemlle XAT =e, ( 19) 
shell 
Dropping exact derivatives with regard to the time we have for the terms in 
En — &m 
w ff { pot(u; + lxs)( ti, + tnx,) — (us + mexs)( tt, + lx) }dr 


=w fff po( lin — lm)x2dr +w fff Po {(uj ie, — ugtt,) + x4(lig — mia’ 
+mu;—lus)}dr. (20) 
We can replace li;—miuj by — lus + mut, and then by (19) we can further 
replace the resulting integral in the shell by 
f f f Pol — ixsus + mx,us)dr. (21) 
Then (20) is replaced by 


w(A—4C)(lin— Im) + | ff po(ujieg — usta; )dr 


+0 | | J potty = {se9( tnx, — ix,)}dr + 20 J | poxe(tnut—lui)dr, (22) 
shell core 
The transformation for the shell makes the second integral equal to that due to 
an addition to the gravitation potential. If we extend this integral to the whole 
Earth, the core integral is replaced by 


w» fff pata vnu; — ius) — us rns, — lg) }dr, (23) 
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‘The integrals through the whole Earth containing wit; and w*u; can be taken into 
account by replacing U, by 
U! = U, + w2xg( lx, + mx,) + wxy( rx, — Ly). (24) 
We can write to the first order 
U;= Uj + Ups (25) 


Uy; = (1x3, mx, — lx, — mx). (26) 
Then 


L [ff pole + at + @)dr=3 [| poll? + h2)(x} + x3)dz 
+ [ff polteritgs + Siy’2)dr. (27) 
The first term is LA(P + m?). 


If 
U = 3(¢,%, + €2%2), (28) 


| [ | atta = dr= —(C—A)(Ie, + me). (29) 


Then in all (integrals being through the whole Earth unless the contrary is stated) 
L= SA(P +m?) + w(A —4.C)(lin— Im) — }w*(C — A)? + m*) — (C — A)(Ie, + mcs) 
+ [Pf polie’ttos + 3i2)dr + [ff pocolutig — usaet)dr 


aU; 
+ ! <—) dr 
{ff Po; Ox; 
Ou; 


+ | [ | pto {aca(siess’ — iss!) — u04( six, — ix.) }dr — 3 | [ | (Ad8,4-+ 26%) ier 


vd 


core 


ie dU, OY Ou; , Ou, 
+ [ff eon? (ua 5 +4202 tah +3 nod =) dr. (30) 


— — > 
OX; ™ Ox; OX}, m Ox), Ox; 


The choice of axes and angular velocities in problems of deformable bodies 
has been discussed often (see especially ‘Tisserand, 1891). The usual choice is 
to make the axes rotate so that the actual angular momentum at any moment 
agrees with that of a rigid body with the same inertia tensor as the actual body 
and rigidly attached to the axes. If the axes are initially principal axes, they do 
not in general remain so throughout the motion. 

We have defined /, m so that they give the actual angular momentum of the 
shell, not of the whole body. It is this choice that reduces the terms in usu; to 
the work done by a potential as in (22), and this is part of the standard treatment 
of problems of elastic rotating bodies. Its convenience is that the u,’ motion 
can be expressed in terms of the normal coordinates of the shell vibrating by 
itself. The longest period of free vibration is about 40 minutes, while we are 
concerned with periods of about a day or longer. Hence iz’? can be neglected 
in the shell in comparison with the elastic terms in u;’. The same applies to 
w(ujits —usu;). The term in i’; vanishes in consequence of (19). 

This is not true for the core. It is obvious that if the core boundary was 
exactly spherical, rotations of the shell would not affect the core at all, so that 
u/=—u,,; and this is in fact approximately true for slow speeds. But the 
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slowest gravitational oscillation of a fluid sphere with the density of the core 
would have a period of about 1 hour (Lamb, 1932, p. 457). A longitudinal wave 
passes through the core in about 10 minutes. Hence in the core we can neglect 
the terms in #; that arise from compression and from motions normal to the 
level surfaces. We cannot however neglect those for motions along the level 
surfaces, since for them the free periods are very long. 

2. Lemmas on small oscillations.—The following principles are often used in 
the theory. 

2.1. Quasi-statical coordinates.—If the Lagrangian function is given by 


2L= 694s + €r39r4s . b,.9r4s + 2¢,9, (r, S=1, 2... n) (1) 
where g, does not appear in L for r=k+1 to n, we write 


- 2 


L'=L-} ean. 2 
bs fet "Fs ( ) 
Lagrange’s equations give 
oL 
~" sxk+1,...%, (3) 


which determine the g, (s=k+1 to n) as functions of the g, and g,(r=1 to k). 
Consider small variations of the q¢,(r=1 tok), which we shall now denote by q,, 
i running from 1 tok andj from k+1ton. Then if (3) still hold for the varied 
q; we have for all variations of q; 


oe aids re =0, 6L’=8L. (4) 
0g; oq; 

Hence if we use (3) to eliminate g; from L’, L’ can be used as the Lagrangian for 
the g;. This result is exact. In the formation of L’ terms quadratic in the gq; 
disappear and linear terms in q; are halved. ‘The result will ordinarily be 
used as an approximation when some of the free periods are very short compared 
with that of the motion being considered. It would, for instance, provide an 
easy way of estimating the correction of pendulum observations for elasticity of 
the rod (Jeffreys, 1956). 

2.2. Exact derivatives with regard to the time.—We often make use of the 
following principle. Dr T. J. A. Bromwich mentioned it in lectures, and it is 
probably well known, but we have not found it in print after a good deal of search. 
If L contains a batch of terms expressible as 


d Fg. 
FS (I++ Iw t)= at a aq, 2” 


the contribution to Lagrange’s equation for q, is 


()- af ae 





dt\0,) ~ dq,6t 24,09," 
all terms of which are immediately seen to cancel. Thus such a set of terms 
may be omitted from the Lagrangian. 

3. Now consider small variations of u;' in the shell. We write 


U,=U,+U, (1) 
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corresponding to contributions due to u; in the shell and core respectively. 


hen an alteration of u,’ in the shell alone does not alter U,, but by reciprocity 


wr ' ,oU, ,oU rf OU, 
ff Poll; ay dr = | {J Pot; * dr + ff Po; Ox; dr +2 {IJ Pot; Ox; dr. (2) 


shell core shell 
The first integral is quadratic in shell displacements. We also put 

Uj =U, + wong {(tn + col xy — (1— com) xg} = 9c}, + 63%). (3) 
Then, all integrals being through the shell, 


. ” au, Q 
3L= | {{ podu; (Se ats “) dr— [ff (AS5in-+ 266) 5—— dude 











Ox; O%; 
a 
+ en —}— "du -} _ 
\] Po h NjOX p. 2 Ox; OX, : “oO : Ox, k 
ov ou’, oY a 
ig a bu; + $u;—— =— uy, dr. (4) 
* OX, Ox; OX, OX; 


This gives, by application of Green’s lemma and some interchanges of dummy 
suffixes, the equations of equilibrium (1(13)), with U, replaced by ¥ and U, by Uj. 

Now suppose that u,’ and its variations are restricted to satisfy the equations of 
equilibrium. Multiply these by }u,’, subtract from L, and apply Green’s lemma. 
Then the contribution of u,’ in the shell becomes 


1 {[ pom’ (U;+ U,)de—3 {| heu/(AIy +2 nej,\dS. (5) 
OX; 


The halving of the terms in U/ and U, is an instance of the general theorem of 2.1. 
The surface integral depends only on surface values of u; and its derivatives. 
It represents work done by the stresses over the boundaries and can be simplified 
by means of the condition that there is no tangential stress. 

The core integrals in lines 3, 5 and the second in line 4 of 1(30) can also 


be simplified. The part arising from U, has already been included in the shell 
integrals. Green’s lemma gives 


: ine 
| ipo’ {Ui + ui, — +}U,}dS 


k 


or ' ou,’ op oF du, 
sa ‘aly | — + gy! 8 al gel iMbe 
ne U; + 3, Ox, +3 u,} (+. Ox; + Ui; Pe ar 2 Jf foo; Ox; Ox, dr (6) 
the surface integral being over the core boundary. We write u’, for the normal 
displacement ; and if gravity at the core boundary is g,, 


or 
Ox, ia —1hgi, (7) 





whence the surface integral is 


JJ pown{U;— des wn + 4U,jaS. ee 
The volume integral would vanish for a uniform incompressible core and can be 
treated as a correction. 


The normal to a boundary is as usual always taken outward from the region 
considered, so that at the core boundary J, has opposite signs according as we are 
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considering a shell or a core integral. With these simplifications 2L is reduced to 
2L = A(l? + m®) + w(2A — C)(lin — Im) — w?(C — A)(P + m2) — 2(C — A)Ic, +: mes) 

+ Sy Pol i;'? + 2i4;'it,,;,)d7 + 2 ff pow(ujits — ugit; dr 


core 


+2 j i} J poco {24(ritu! — lu’) — ui( tnx, — lx.) }dr 


core 


+ | J J pot’ = (Uj + U, dr — | | 1,14; (AAS 4 + 20,4 S 


shell shell 


i | J pott,(2U; —g,u), + U,)dS 


core 
~ fff (20; + ue = +U,) (9S + ou’ 28) dr 
core 7 
— [Jf pow Pott; ‘s m4 — star {fj AA*dr. (9) 
core core 
4. Form of solution for the shell_—We take briefly 
U, =cK,=cr’S,, (1) 


where ¢ is a constant and K, a solid harmonic of degree 2. At present we ignore 
the difference between U, and U}, which can be adjusted later when required. 
We also take 





U,+ U,= KceKg. (2) 
K is a function of r and is a number. We introduce the variable 
g =r/a, (3) 
where a is the outer radius, and write 
xu; =q=q(&)Ko. (4) 
q(é) is anumber. K’(€) is discontinuous at the boundaries, and we have 
={K(Q-1}eK, €<1 (5) 
={K(1)—1}(a/r’cK, €>1, (6) 
and at the outer surface 
1), 
[ = | = —47fpol Wk (7) 
This leads to, at €=1— 
[ek 
g + 5tK(1)— 1} | =47fpo(1)9(1). (8) 
We can also take 
u.=P(™)'cK, & 
ind SS, eating (9) 
U,=QcK, E<a. (10) 


Then we have 


aU 2 ,. fU\ _ (aU, K q(x) 
(=). -3 = PcKy; ’ (=) = (3 Mi: + 47fpo(a— reed (11) 
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au. 2 a aU, om ou, q(x) ‘ 
(=). = precy (=). a (= a + 4nfp («+ ) — K, (12) 
and by comparing the values of U, and 0U,/ér at €=« + we have a pair of equations 
for P and Q, the solution of which is 


5Pc= —acK'(x+) + 47fpo(x + g(a) (13) 
sQc=acK (a+) +5{K(2)— 1}e— 4nfpa(a+)a(a). (14) 
Also 
ac| Se | =42fale)inde+)—eol2—)} (15) 
For the displacements we take 
= FO st + Oak, (16) 





This is of the form of the statical solution and has the property of orthogonality 
with rotations used in Section 2. F and G are numbers. Let accents denote 
derivatives with regard to €. Then 


q(E) =2F(E) + &G(E). (17) 


Again, with /,,=x,/r, we have 


1,(AA'3;, + 20€%) = - Te Ke oF + FE+ Ge) 


i [A(2F’ + 8G’ + 5€G) + p(2F" + 2€G’ + 4€G)). (18) 
Since there is no tangential stress at the boundaries 
EF’ +2F+@G=0, €=1,4% (19) 
and the elasticity contribution to 2L is 
— @EM[A(2F’ + &G' + 5€G) + u(2F’ + 2&G’ + 4€G)]9(é), (20) 


where 


M= ff S3dQ=4n (21) 


when K, is the solid harmonic xz or yz. 
For €=1, (18) is correct. For =a, it must be taken with the opposite sign. 
For integration over a sphere we have the identity 


(Gey as =") || Kzas. . (2a) 


U,+ U.=c{1 + P(a/é)} Ko, (23) 


whence the volume integral in 2L is 


Now 


Matc | Pe (oF + 22G — 36P ede. (24) 


(18) is already expressed in terms of boundary values of F, G and their first 





sa eee: ‘ee 





sitar sae ‘ 





No. 2, 1957 Theory of nutation and the variation of latitude 151 


derivatives. (24) still contains internal values of F and G, and the corresponding 
expression requires a numerical integration. ‘Takeuchi solves the differential 
equations for F, G and K (the last being Poisson’s equation) in terms of their 
values and their first derivatives at€=a+. His functions, however, are multiples 
of those used here; in particular he takesc=1. In view of the pair of boundary 
conditions (19) we can eliminate two of the constants; and on account of the 
importance of q(1) and g(«) it is convenient to eliminate two others in favour of 
them. Further, c is a datum, and we have the relation (8) connecting K’(1) and 
K(1). One constant remains arbitrary, and on account of the prominence of 
P it seems best to take K’ (« +), which is wholly determined by the displacements 
of the core. When this is done the internal values of u,’ in the shell will have 
been eliminated. 

If we indicate Takeuchi’s functions by the suffix 7, the correspondence 
between his and ours is 


F(é)=ca®F7(€); G(E)=ca*G7(£) ; \ (2s) 
q(é)=ca(2Fp+Gp);  K(£)=4nfa?Ky(€). 
Since his equations are all homogeneous the factor a* can be omitted. We write 
c/4nf =D and take the correspondence to be 
F G q DK 
Fr Gy qr Ky 
D and p are densities, and where they appear in the numerical solution they must 


be interpreted as (D, p)/(1 g/cm’). 
Then (8) can be written 


D{K"(1) + 5K(1)}=5D +p(1)q(1). (26) 
Also 

2F(1)+G(1)=4(1), (27) 
F'(1)= —q(1), (28) 

and with «=0°5447, 
2 F(«) +0°:2967 G(x) =q(a), (29) 
0°5447 F"(x)= —q(a), (30) 
DK'(«) =cK p(x). (31) 


Dr Takeuchi very kindly gave us a copy of his matrices for his Model 2 
connecting F(é), G’(é) ... A(€) with their values at =a, which are not given in 
detail in his paper. These were transformed to the four constants used here, 
by means of the six equations above. ‘There are discrepancies between the last 
lines in his (198) and (205); the second entry should be 0-072887, the fourth 
0°23256 in both cases. 

The solution is as follows : 


F'(a) = — 18359 q(«), 

F'(x) = 24°6632 (1) — 10°2459 g(x) —0'00013 DK’(«) + 18392 D, 
DK'(«) = DK'(«), 

F(a) =0°31942 g(1) +0°18614 q(x) —0°00062 DK’(«) + 0°03546 D, 

G(«) = — 2°1524 g(1)+2°11499 9(%) +0°004145 DK’(x) —0°23895 D, 
DK(«) = 0-00078 (1) +0°3555 9(«%)—o°'1110 DK'(a) + 1-1020 D, 


(32) 





II 
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F(1)= —4(1), 
G'(1) = 3°50477 4(1) — 2°09580 g(x) — 0°00325 DK'(x) — 0°1476 D, 
DK'(1) = 1-08059 q(1) — 0°59626 q(x) + 0°02873 DK’(«) — 0:22403 D, 
F(1)=0-12191 g(1) —0°06591 g(x) —0°00058 DK'(x) + 0°04254 D, 
G(1) =0°75622 g(1) + 0°13178 g(a) +.0:00118 DK'(x) — 0:08488 D, 
DK(1)=0-38367 9(1) +0°11938 9(«) —0°00578 DK’(a) + 104472 D. | 

Some features of the solution are worth special mention. The coefficient 
of D in K(€) never deviates far from 1. _ If motions of the core and the boundaries 
were prevented, the displacements of the shell would have little effect on the 
potential. The coefficient of g(1) in DK(«) is very small. With a homogeneous 
incompressible shell and core, and with q(«) =o and no external field, K(£) would 
be constant. This is the problem of distortion by loading over the outer surface, 
with a constraint preventing core motions. But here the coefficient of g(1) in 
K(€) varies by a factor of 50. In a uniform sphere the effects of compressibility 
are usually not large, but it appears that in this problem they are. The large 
coefficient of g(1) in G’(«) corresponds to a large disturbance of density just 
outside the core, and it appears that the reduction of potential due to this cancels 
most of the increase due to the elevation of the outer surface. The coefficients 
of g(a) and DK'(«) in DK(€) decrease greatly outward; this is mainly due to the 
factor €~ in the contribution from P. 

Extra figures were calculated in some of the coefficients because of cancel- 
lation in the checking. 

The solutions for other values of € were calculated on the EDSAC machine 
at the Mathematical Laboratory, Cambridge, by Miss M. Lewin (now Mrs Mutch). 
As Takeuchi’s last two figures are mainly guard figures they are rounded off at this 
stage, and the solutions are in Table I. ‘The values of € are as follows: 


fo=%=0°54474, €=0°59184, £2 = 063893, = 0°68603, 9 £4 = 0°73313, 
£, = 078022, 6 = 0°82732, 7=0°87441, g=O'92I51, £y=0°96075, 
£10 = 1°00000. 


, (33) 





The intervals for r, measured from the outer surface, are 250 km to 500 km, then 
300km to 2g00km. ‘The data assume a discontinuous change of properties at a 
depth of 500km. ‘The rows in the matrices give F’, G’, DK’, F, G, K for each €. 

In the integration we need only F and G, but some of the coefficients vary too 
rapidly for accurate integration by the Gregory formula. Since, however, F’ 
and G’ were also available, the interval was halved by the formula (Jeffreys, 1953) 


F(dh) = 3f(0) + 3 f(A) — ght f'(h) -f'(0)}. (34) 
This formula is more accurate than the Newton—Bessel formula with second 
differences and can be used at the end intervals. When this was done integration 
by Gregory’s formula was accurate enough. 
After integration it was convenient to use (13) to eliminate DK(a+); with 
the values adopted po(« +) = 5-53 g/cm’, 
DK"(%+)= —9:1787 DP + 10°152 g(a). (35) 
Then the surface integrals give in 2L (lines 4, 5 of 3(9)) 
— Ma* x 10" 9(1){8-415 9(1) — 3-168 g(x) —0:00670 DP — 0-7698 D}, (36) 
— Ma? x 10" 9(x){ — 3-157 9(1) + 1°644 g(a) +0:0631 DP—o-1047 D}, (37) 





2 
BS 
% 
; 
7 





Py wocinnntscpaagt 


No. 2, 1957 Theory of nutation and the variation of latitude 153 


and the volume integral (line 4) gives 
Mea* x 10!*(0-7401 g(1) D+0-1003 q(x)D +0:0164 D?P+0-0744 D® 
+0°0068 g(1) DP—0-0621 q(«) DP+0-0017 D?P*. — (38) 
C.g.s. units are understood. 

We have two useful checks. (1) and g(x) are Lagrangian coordinates, and 
by a reciprocity theorem the coefficients of g(1) g(«) in the two surface integrals 
should be equal. ‘They agree to about 1 part in 300. Also D and DP represent 
external forces, acting on the shell, and if g(1) and g(«) are put zero the terms in 
D and DP will represent work done by the external forces with the boundaries 
fixed. Then the amounts of work done by D on the displacements due to DP 
and done by DP on those due to D should be equal, by another reciprocity theorem, 
and the coefficient 0-0164 is in fact the sum of two of 0-0082. ‘These comparisons 
are checks both on the present work and Takeuchi’s. 


TaB_e I 

q(1) q(x) DK"(a) D 
Y, —o0'069 —1°296 0*0002 O°O14 
15°198 —8-770 — 0009 1°056 
0°486 —o'468 o-611 —0°042 
0°317 O°1ls — 00006 0°036 
—1°'238 1663 0°0039 —0°173 
0°020 0°343 —0°074 I‘IOI 
le —o'211 —o'892 00003 0014 
9°829 —6-918 —oO°’orl 0°631 
0°723 —o'620 0°388 —0°073 
0°310 0'065 —o0-0006 0°036 
—o'662 1293 0°0034 —0°134 
0°049 0°316 —o'osi I*I00 
Ys —0'297 —o'661 ©0002 0'016 
6613 —5°296 —o'o1o 0°390 
0°847 —0°634 0°255 —0°097 
0°298 0029 —o0:'0006 0°037 
-—0'282 1°007 0°0029 —oO'lil 
0086 0°287 —0'036 1°094 
Y, —0°371 —0'503 00001 o-o18 
4°610 —4°054 —o-008 0°249 
0'922 —0°593 0'173 —O'lI7 
0°282 0°00! — 00006 0°038 
—o'0o2I 0°789 0°0025 —0'096 
0128 0'258 —0'026 1‘089 
 £ —0°437 —0°394 ©0000 o°01g 
3°369 —3°139 —0'007 o'161 
0'972 —0*540 "120 —-0°134 
0°263 —0'020 —0o0'0006 0°039 
0°163 0°621 00022 —0'087 
O°172 0°231 —o-'o1g 1083 
Ye —O'514 —0'310 —0*0000 0018 
2°599 —2°467 —0°005 0-099 
1°009 —0°493 0°085 —0°149 
0°240 — 0036 —0'0006 0°040 
0°302 0°490 00019 —o-o81 


0°219 0°207 —O°O15 1'076 





154 Harold Jeffreys and R. O. Vicente Vol. 117 


q(1) q(x) DK'(«) D 
Y; —o'589 —0'239 —0*O00! O'O17 
2°163 —2'032 —0'005 O'O5I 
1°037 —o°461 0-062 —o'164 
0°214 -~0°049 —0'0006 0'040 
0°413 0°385 00016 —0°077 
0°267 O'185 —o-oll 1°069 
Y, in —0o°667 —O°175 — 00001 0'016 
2°105 —1°832 —0'004 0°003 
1°087 —0°454 0046 —o'180 
0184 —o'058 —0-0006 0°041 
o’511 0°296 00014 —0:'076 
0°317 0°163 —0°009 1°061 
Y, out —0'603 —o'184 —0'0002 0'028 
2°949 —2°177 —0°004 —0°097 
0626 —0'531 0046 —o'1go 
0184 —o'058 —0'0006 0'041 
O’511 0°296 0°0014 —o'076 
0°317 0163 — 0009 1°061 
Yo —o-’788 —0°094 — 00001 0°O17 
3°090 —2'072 —0°003 —O'1l3 
0856 —0o*560 0'036 —0°205 
O°157 —0'064 — 00006 0°042 
0°628 0°213 0°0013 —o'o80 
0°346 O°142 —0'007 1°053 
Yio — 1°000 —0*000 —0*0000 — 0000 
3°505 —2°096 —0°003 —0'148 
1081 —0'596 0°029 —0'224 
o°122 —0:'066 —0:'0006 0°042 
0756 0°132 00012 —o'085 
0°384 o'11g --0°006 1045 


5. Models for the core—'The actual core is too complex for detailed treatment, 
but useful simplifications are possible. It has appeared probable that so long as 
the moments of inertia of the core and the ellipticity of its boundary are correct, 
the detailed distribution of density is of secondary importance. But if the 
Radau approximation is valid, the mean moment of inertia is enough to determine 
the other quantities. Now Bullen (1942) found that extreme hypotheses gave 
for the core 

I 
We =0°389 and 0-373 

in comparison with 0-400 for a homogeneous core. Now it has been shown in 
a former paper (Jeffreys, 1951 a), in which K = corresponds to J/Ma?=1}, that 
the ratio e/m for a Wiechert model varies over the whole range of possible core 
radii by only 2°5 per cent. For J/Ma*=0-4, the body is homogeneous and the 
ratio does not depend at all on the radius that we might choose to be that of the 
core. ‘Thus it seems likely that for the Earth’s core (in which the inner core 
would play the part of the core of a Wiechert model, and the potential due to the 
figure of the shell would be added to that due to rotation) the Radau approximation 
will always give solutions valid within 1 per cent. 

Through most of the core the variation of density is mainly due to compression. 
If this held for the whole of the core a smooth distribution of Laplace’s or Roche’s 





% 
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type would fit fairly well. It is desirable also to have an estimate of the effect of 
the inner core, but for simplicity we shall treat this as a point mass at the centre. 

Bullard’s (1946) density distribution makes the moment of inertia of the core 
0-107 of that of the whole Earth, and the ellipticity 0-002567. The dynamical 
ellipticity is 0-002562; this was computed from his solution by numerical inte- 
gration. We shall adopt these values in all models. The reason for the close 
agreement is that the dynamical ellipticity is not a weighted mean of the internal 
ellipticities in the usual sense. In computing it we compare integrals of the 
forms fpd(r’e) and fpdr’, whose ratio would be e if e was constant. With 
variable e, the smaller e near the centre reduces the first integral, but this is 
compensated by the fact that the consequent reduction of r’e at intermediate 
depths has to be made up by increased values of its steps as the boundary is 
approached, so that the dominant part of the integrals still comes from values 
near the surfaces. 

We consider two extreme models for the core. In the first the core is treated 
as homogeneous and incompressible, with a point mass added at the centre. In 
the second we take a Roche density law 


p=hy— kh, & (1) 


where a,=azx is the radius of the core and £,=r/a,. ‘This makes the mean 
density 





p=ky— 3h, (2) 
and the moment of inertia is that of a homogeneous sphere of density 
Po= ky — $hy. (3) 
Bullard’s values for the density lead to 
p=10°70, p,=10°28 g/cm’, (4) 
whence 
ky = 12-91 g/cm®, k,=3°68g/cm’, k,/k, = 0-285. (5) 
For a distribution of density p, with a central particle of mass m, the mean density is 
* m, 
= y ae 6 
P=PeT ina’ ( ) 
whence 
“i_ =0°04 (7) 
5 =O°04I. 
4mp,a} 


We note that a 28 per cent variation of density in the Roche model is equivalent 
to putting 4 per cent of the mass in a central particle, so far as the hydrostatic 
theory is concerned. 

More can be said; for Radau’s approximation is good for any Roche model, 
and it is easy to show that for any model satisfying it the surface ellipticity of a 
free body is given by 


e I 


(Ga) (8) 


to the first order in the small quantity 5//2Ma?—1. But for the central particle 
model we can easily show that the same is true. Also to the same order the 
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dynamical and geometrical ellipticities are equal. ‘There can therefore be no 
inconsistency in adopting Bullard’s ellipticity of the core, namely 0-002567, in 
both models. 
For Roche’s models we take the undisturbed density over a level surtace to be 
Po=ky— ki €, (9) 
where a,é, is now the mean radius of the surface, and the corresponding geo- 
potential is 





Y= — 4nfai(bhofe — aphi€2)- (10) 
By Radau’s approximation the ellipticity of a stratum of uniform density is 
( k ) 
ema I~ BE (1-H) f (11) 
: 0 
where e, refers to the boundary. Then the equation of the stratum 1s 
r/a,=£€,{1 + e(4 —cos* 6)} (12) 
and we put 
W,=17(4 — cos? 6) = 4(x? +.y? — 227). (13) 
Then to orders k, and e, 
P wh 
p=hy-k 4 +—>W, (14) 


1 2 
Y=- jaf (ho? th.) + $7f(Ro— =k, eV o- tafkyes W. (15) 


We shall not need the ellipticities of surfaces of equal density again, and shall 
therefore denote the ellipticity of the core boundary simply by e. 

«6. Solutions for central particle models: statical case.—In this tase the volume 
integrals in 3(9) disappear; and 


: ax\* , 
U.= tnfoata) (=) K, (1) 
whence 
DP= $p.q(%) = 2'0564(x). (2) 
Also if the normal displacement is 
u, = 4(%)K/aa, (3) 


| | lip," {2 U,+U, +4; a dS 
: 


=f f pcug(2U, + U,—ygyup)d 





= 4nfa®. Ma. «5p, (2p —2,q(x) — = 5 a()) q(x) 
4770 
=0°1677 x 10'*Ma*{2Dq(«) — 1°51 g(a)}. (4) 


We add this to 3(36) (37) (38) to give the potential terms in 2L. We take first 
the statical case. We have in all 


2L . 
i Ma? x 108 = ~ 8415 9(1)? + 6-353 9(1) q(x) —2°147 9(x)? 
| 
+1° 
whence 1°510 g(t)D + 0°574 q(x)D + 0:0744D (5) 


W(1)=0°317D, 9(x)=0-603D, K(1)=1:289, F(1)=0-0444D. (6) 
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Love’s numbers A and k are 


gq) ,_ 41) 
k= K(1)—1=0°289 (8) 
and the Shida—Lambkert number is 
= 
l=4p D =? 082. (9) 
For comparison ‘l'akeuchi’s values are 
h k l 
Model 1 0°587 0°290 0-068 
Model 2 0610 0:281 0082 


Agreement is good, but as the present model for the crust is nearer to Takeuchi’s 
Model 2 than to Model 1, it is rather surprising that two of the numbers agree 
better with Model 1. 

The factor M has different values for different second harmonics, but cancels 
in the statical case, which is applicable to the semidiurnal and long-period tides. 

7. Central particle model: dynamical case.—When |, m are not zero t,' can no 
longer be neglected in the core. ‘The statical solution for the shell must be 
modified by restoring U}. 

u’. can be broken into two parts, 








Uj) = Vj, + Vay (1) 
where 
L,za, m,za, (1.x+my,y)e, 
Uy, = ~ * C a a ’ (2) 
1 1 1 
De Iza, mysa, (1,x+myy)c, (3) 
“ar Cy ’ Cy ’ a, ° 3 


This makes u,’ linear in the coordinates; this is justified because the obstruction 
due to the central particle will be of no importance since u,’->o at the centre. 
The normal displacement at the boundary is 


5 
uy, = ae, te (4) 
where 
s #+/f # 
= (5) 


m= t+: 
R? at ct 
u;, must be continuous with the normal displacement of the shell; since this is an 
elastic displacement R can be replaced by a,. 
We introduce the constants 


Ay= ff pyx?dr, Az = Sf Podr ; (6) 
C,=2Ay, A;= Ay + Ags, F=2" Ay=2- Ay. (7) 
Then 
F? = 4A, Ag3=C,(2A,—C;). (8) 
The ellipticity of the core being small, we have nearly 
C,=A,(1 +e), 24,—-C, = A\(1 —e) (9) 





158 Harold Jeffreys and R. O. Vicente Vol. 117 


and hence 
F=A,(1 + Ofe?)). (10) 
We find vm 
[| | pig’ dr = Ay(lE + B+ ant + snd) + 2(Cy — Ay(lile + rnyrng) (11) 
2 | | | pityglt! dr = 2F (Il, + mmy) (12) 


2w fff Po( tits — Usit; dr = wC (Lm, — ml, +1,m, — ml, + l,m, 


— Me], + latit, — Maly) 

which is equivalent to 
The term in /,m, can be dropped, since it is of the order of the kinetic energy of 
the elastic displacements. But those in /,m,—m,/, are of the order of those 
arising from the correcting terms in U;, and must be retained. 

2 | | | pow {sg( seu; — lus) — us(tnx, — Ixy) }dz = 2 Fex( mal, — Im,). (14) 
The core integrals arising from the last three lines of 3(g) are the same as in the 
statical case except that U, has to be replaced by U}, and with this change 6(4) 
will hold. 


We take 4,;=0'107C, C,=A,=eA,, e=0-002562; 


@,/a=“%=0°5447; w=7'2921 x 10%/sec; C=8-05 x 10! gm cm*; 
A=C(1—H), H=0-003273; 


a= 6°371 x 10°cm. 


Then Cw = 4°28 x 10° gm cm?/sec? ; 
M=7, Ma*=2-166 x 10% cm’, 
Ma* 


10!2 ——. = 5061. 
Cw ° 


We keep H and e explicit for the present; this facilitates some checks. We take 
also 
Q(i)Ky=ryxz+rgyz 
Q(a)Ky=s,xz+syyz; le=4s,, my=4s5, 
DK, =d{xz+dyyz 


w® (fe m =f 
= A (S42 +l) ase (S —_ +m) yz} (15) 
and w*/4mf = 0006345. 


It is convenient to use 2L/Cw? instead of 2L. ‘Then 
a = — 425°Q(rj +13) + 321°5(145; + 7952) — 108-7(s? + 82) + 76-4(r,d) + ry) 
+ 29°05(s,d; + 8.d3) + 3°765(d/? + d;?). 
+ A [(1— AYE + in2) + 02(1 — 2H)(bin— in) — (2 + m2) 


+0°107(/2 + m:) + 0-107 e(l,i, + 159) 


+0°214(Ih + tintin) +.0°107(1 + €)co( lyri, — myl, + syrin, — Syl) 
+0°2140( ml, — lm,) — 2H (Ic, + me,)). (16) 
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The terms in el,s,, wes,m, can be dropped in comparison with those that do not 
contain e. 
We put Cy+tg=k, l+im=l, 1,+im,= ,. (17) 
Derivatives of r,, rg, 5;, 5, do not appear in the Lagrangian; hence these 


coordinates can be eliminated. 


We find 
(71, 72) =0°3774(51, 52) + 0°0897(d;, 3), (18) 


» an, OOBEIIS,. : 
(51, 52) = 0°607(d;, dy’) + ———.(m,, — 1), (19) 


and the modified Lagrangian is given by 


L’ : , 
~ = (1 —H+0-00100)(/? + m*) + w(i — 2H + 0:00200)(/m — Im) 





— w*(H — 0:00100)(/? + m?) 
+ 0:10706(/2 + m®) +0-21441(Il, + min,) +0" 107w(1 + e)(/,m, — m,l,) 
+0°2140(ml, — im;) + 0°00041w(lm, — im) 
+ (—2H +0°00200)(c,/ + cm) + — (cym— cl) 


4 SOCOE (catty — Cah) + POP? (cz + 2), (20) 

We form the equations of motion and assume a time factor e”*. ‘Two speeds 
of free motions are found to be o and —w. ‘This gives a valuable check, since 
the speed o corresponds to a statical displacement of the core without moving the 
shell and that of —w gives a steady displacement of the axis in space. A further 
check is that for y= —w the rate of motion in space of the axis, given by (y+), 
is found to be the same as for a rigid body. ‘This represents the precession, 
which in all previous models has been found to be independent of the internal 
constitution. 

For the other free motions we find 


z =0°00255 ; . = 3924; = —0°99931, (21) 
a  —=- ; sen P 
tale a Sg a Tee (22) 


The first of these is the variation of latitude. ‘The period is shorter than any 
found from observation. It would be lengthened by the effect of the oceans; it 
remains to be seen whether compressibility of the core would make much 
difference. ‘There is very little motion of the core, since apart from the effect of 
the ellipticity the displacements due to ¢ and ¢, nearly cancel.* 

The second movement is not observed. The core displacements are opposite 
to the general rotation and much larger. 


* A former solution (Jeffreys, 1951 6) by the method of undetermined multipliers gave 463 days 
for the period. A homogeneous core was assumed. This differs from the present model only 
by the central particle, and a great difference was not to be expected. The difference has been 
traced to a misplaced 2 in the earlier solution, which when corrected leads to 390 days. 


' 
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For the forced motions we write y+-w=n. As a standard of comparison we 
take the value of ¢ for a rigid Earth and call this ¢. We have 


Hk= —nl,(w—n+nH), (23) 


and the equations of motion become 
n 
(w —099773n)£ + 0°10720(w —m)C, = (w — 0-996730) (: — 0°3055 “) for (24) 
0:10720n£ + (0:000214w + 0°10706n)l, = 0061 1n (: —0°3055 =) Co. (25) 


The coefficient 0-3055 is 0°00100/H, and is a consequence of elasticity of the 
shell, not of fluidity of the core. 
Results for the most interesting speeds are as follows: 





n/w Clee Clee ‘Tidal component Nutation component 
— — 10768 —0'509 OO Secondary fortnightly 
7 
_ * 10895 —o'819 
_ a 0°9964 0°034 K, (companion) Principal nutation 
° 1*0000 0°000 K, Precession 
zm 10036 —0°034 K, (companion) Principal nutation 
: (correction) 
183 1°0350 —0°342 P Semiannual 
= 1°0269 —o'458 O Fortnightly 


For n/w= — 1/447, ¢/Cy and ¢,/%, become infinite and change sign; and for 
n/w = — 1/479, ¢ vanishes and ¢, does not. This is the explanation of the peculiar 
behaviour between — 1/183 and — 1/6800. 

On the whole the results are fairly similar to those found for a rigid shell. 
The principal difference is that for a rigid shell, unless |n/w| is small of order e, ¢, 
is nearly —{ and ¢/{, is about 1-12. This is not true with an elastic shell. Even 
for O and OO {, is only about —o-5 %, and for the principal lunar diurnal tide O 
the predicted amplitude of { is under 3 per cent in excess of the rigid body value. 
H. R. Morgan (1952) has made a thorough analysis of the observations and has 
finally decided that they will not determine any significant correction. 

For a rigid shell the correcting factor to the principal 19-yearly nutation was 
09940. The result here is 09964. For comparison, Newcomb’s standard 
value of the nutation in obliquity is 9"210. The rigid body value based on 
Spencer Jones’s values of the solar parallax and the lunar inequality is 
9":2272 + 0”-0008 (p.e.) Rabe’s parallax and my value of the lunar inequality 
lead to nearly the same results. Recent observed values have mostly been 
about g":207. Corresponding to Spencer Jones’s calculated value, the amplitude 
in longitude would be 6’-872. A determination by E. P. Fedorov (1952) separates 
the corrections to the amplitudes in obliquity and longitude, and indicates that 
the former should be 9’-190 to 9’:200, which would agree with a correcting 
factor of 0-996 to 0-997. In longitude he gets 6-874 to 6”-886, larger than the 
rigid body value. But theoretically both components should be reduced, that 
in longitude by the larger factor. There is therefore still an inconsistency. 
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8. The bodily tide numbers for diurnal tides.—To get these we have to express 
fy in terms of k/w* and then express r,+ir2, s;+%s, in terms of c,+ic,. The 


results are as follows: 


n/w 
I 


~ 13°7 
I 
183 
I 
6800 
° 
I 


6800 


183 
I 


13°7 


h 
0°590 
0°523 


0°490 
0°492 
0°494 


O°S355 


0°584 


0°244 


0'218 


0°205 
0°206 


0°207 


0°231 


0°242 


l 


07082 
0°084 


0086 
0086 


0°086 
0082 


0082 


lis practically the same as on the statical theory. 


1—h+k 


0°654 
0°695 


e715 
O'714 


0°713 
0°676 


0°658 


1+h—gk 


= 


= 


224 


‘196 
“182 
"183 
“184 
"209 


‘221 


h and k are substantially less, 
except that A approaches the statical value for the tides O and OO. 
values of 1—h+k and 1+h— 3k are 0-704 and 1°152. 


The statical 
These are the factors 


associated with observations of the bodily tide, and it appears that K, and its 
companions are likely to be nearly the same as on the statical theory, but there is 


some hope of detecting a difference for O. 
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THE THEORY OF NUTATION AND THE VARIATION OF 
LATITUDE: THE ROCHE MODEL CORE 


Harold Jeffreys and R. O. Vicente 
(Received 1956 December 22) 


Summary 

The theory given in a previous paper is modified by the use of a core with 
the density a quadratic function of the radius, the variation of density being 
wholly due to pressure. Three pairs of cubic terms are introduced into the 
displacements. The period of the Eulerian nutation would be 395 days, and 
the amplitude of the 19-yearly nutation would be multiplied by a factor 0-9975. 
Allowance for the effect of the ocean would increase the former period to about 
430 days, and in both respects there is now tolerable agreement with observa- 
ations. Correcting factors are given for the other principal nutations. 
Those for the fortnightly nutation do not differ much from those for the 
central particle model. There are substantial changes for the semiannual 
nutations, owing to changes in free periods that occur in their neighbourhoods. 
The 19-yearly nutation in longitude remains anomalous. 





1. Ina previous paper (Jeffreys and Vicente, 1957) we have used ‘T'akeuchi’s 
solution for the shell; the core was replaced by homogeneous incompressible 
fluid with an extra particle at the centre, adjusted so that the mass and moment 
of inertia of the core agree with those used in Bullard’s theory of the figure of 
the Earth. In the present paper we consider as an alternative a core witha 
distribution of density of Roche type, 

Po = ky — kyr? /a/, (1) 
the variation of density being considered due wholly to compressibility. The 
displacements are no longer linear functions of the coordinates. In a previous 
paper (Jeffreys, 1950) a suitable set was introduced, containing cubic functions 
of the coordinates. ‘This is somewhat modified in the present paper. 

The change of density at any point is —u,(dp)/0x,) — po(Au,/Ax;). If the core 
boundary is not deformed and there is no change of density at internal points, 
the gravitational potential due to the core is unaltered. In the central particle 
model the(/,, m,) displacement is along level surfaces and satisfies these conditions ; 
but since in the Roche model the ellipticities of the level surfaces are not all equal 
this might need modification. It is found that independent displacements 
satisfying ie am, 
m5, = O35 9 =o (2) 


are, to order k,, 


6k, , 6k, *\  \,ek,2? 
oyna {rte(r— St + tz) + Be am) t=1,2 


35k) | 35k oat] © ® yay (3) 
6k, 6h, r2\_,, ek,3® 
y= —(4x+myy)<1- —~-—= +—5)+#=5} i= 
neem vy e(: sah, * show) +E Re} ie 
ou= (a(1-+6)+ 2008} (A?) i=1,2. 
i] 


03; = r(1—e)+2e2*)} i=3 
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A further simplification may be made, for we are neglecting k? and e? and may 
therefore naturally neglect kye. ‘Then v,; reduces after all to the form for a 
uniform core. In v3; we have already dropped ,, since vs; is a consequence of 
non-uniformity in density and its reaction on /,, m, will contain a factor k,; we 
retain e, however, because /,/l,, m,/m, are not small in some motions. These 
two parts generalize * : (/,, m,) displacements of the homogeneous case. 

In addition we must generalize the (/,, m,) displacement to include some 
cubic terms. If the change of density is 


I r? 
Ap= a Pi + Peas Ky, (5) 
and the normal displacement at the boundary is 
v= 12K, (6) 
a 


there is effectively an extra surface density 
a), 
ao =(ky—k,)u,, =(Ro— Ry me Kg. (7) 
The change of external potential is found to be 


a 5 
U.=47f (bo hi)a(a) +4014 hed (4) Ky (8) 


and that of internal potential is 
— a ; 
U.= af {ky k,)q(«) + (: at ) Ait (: te =) ps} Ky. (9) 


If the displacements are 


r? OK, 5x; - 
wm (B45) ae + HR oo 
we find 
: K Pr 
Ap= - (ho(4y + 58)—48hy} + hy(8y + 78) Ke (11) 
1 
x; r?) Ky 
w= Fu= 4 af+ (ay +8) b=. (12) 


For the choice of component displacements we have the following considerations. 
q(x) is a shell coordinate and may naturally be adopted also as a core coordinate. 
We can take a linear form to give this. We take two other forms, both giving 
no normal displacement at the boundary, of which one contributes nothing to 
du;/Ox;, but gives a radial displacement at internal points, while the other con- 
tributes to ou;/dx;, but gives no radial displacement at internal points. ‘These 
are all essentially similar to the elastic displacements in the shell and their products 
by e can be neglected. 
We take (writing only one harmonic) 





Vai = 29(%) 5» (13) 
pe Se a \ eke stn 
Uy=l, (1-35) Ox; er .), (14) 


(15) 
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Then 


w= Ke (16) 

OX; 

r? a” 
qKy= xju;=TUu,, = {yay an( - =) bK (17) 
1 
py, = 2k, g(a) + 4hyly + ORol; = 2k,g(«) +0, (18) 
— 4h, 1, — 6k,/;, (19) 
k : 
EPit+ =Pe= = g(a) + gap Rala+ 3 (Ro — ZR, ls. (20) 
Then, up to cubic terms, 

U,) = Vy; + Vg; + Vai + Vg; + V5 i: (21) 


We need also an expression for A, the bulk-modulus. Since the variation of 
density is attributed entirely to compression, we have 


dp_ dV, k(n \? hy 
nap =p =hafatp(1- fe) (1-274), (22) 


where &,=r/a,. (A/p)"* varies from 8-7km/sec at the centre to 6-2km/sec at 
the boundary. The seismological values are 11-3 and 8-1 km/sec. The difference 
is due to the use of a Roche model, which takes no account of the inner core. 
With this model, subject to the mass and moment of inertia of the core, the above 
expression for A cannot be altered. ‘The correct values of the compressibility 
are about 0-6 of those given by the Roche model, but since we have already the 
solution for the incompressible model we may expect that interpolation will give 
a result not far from the truth. _In integration we must use the full expression 
for A, since, though k,/k, is treated as small, 2-6k,/k, cannot be. 
2. The potential terms.—We write as before c= 47fD, and find 


| | eow,(2U, ~g.m, + U,)dS 
= 4nf(ky — k,)Maiq(ak2D— 2 (ky —*hk,)q(a)+ 3 aR 
+ 5( sean (1) 
~ {{f(2er+uS +0.) (9M +2) 
= — 47fMa}[2D{F * = q(a)+ 2 Ry I,+$(Rg—- “ R,)l;} 
+ 5 Rokig*(%) — x2 Roki q(a)l, 


+ ARo(Ro i $hy)q(a)l; + sorkoki li — 63 is olals +% 5 Roll, (2) 
re LOY ou’ 2 
-|| | pot — s# dr = — Saf Maily| (12 — Skok, 9(2) (3) 
Aj k 
+ 2l, & hi — ie hoki), 
: ou 
-(ffa ( = ) dr = —4nfMa’. 26522, (4) 


where k, and k, have been given their numerical (c.g.s.) values in the last expression. 
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We have also 

DP = 2:056q(«) + 1°723/; + 0°0935/4, (5) 

DK'(«+)= —9:179DP + 10°152q(«). (6) 


The shell terms given in 4 (36) (37) (38) of the previous paper must be expressed 
in terms of g(1), g(«), 4, and /; and added to these core terms. ‘The total is 


2L = — Ma® x 10"[8-4159°(1) — 6°3539(1)q() — 1°5369D4q(1) 
— 0°02339(1)/, — 0°001269(1)/, — 0°5740Dq(x) 
+ 2°13897(«) + 0°0373q(a)l, + 0°-629q(«)l; 
—0°0167DI, + 0:0239/7 


—0°389DI, — 0:2368/,/, + 4:58/2 — 0:0744D"]. (7) 

This is stationary when 
g(1)=0°324D; g(a)=0°618D; 1,=-—o0-129D; 1,= —0-0025D. 
These lead to the bodily tide numbers, for the statical solution, 
h=0°598; k=0:273; /1=0-082. (8) 

In comparison with the central particle model, g(1) and q(«) are slightly increased, 
but /, and /; make small negative contributions to k, which is slightly diminished. 
/ is unaltered. 

However, it is easy to see that the differences can represent nothing but minor 
errors of approximation. We are considering similar constitutions of the shell ; 
and the core in the two models has the same mass and moment of inertia. Suppose 
the same displacements applied to the shell in both models. ‘Then the potential 
acting on the core from outside it is the same for both, and therefore, by the 
Radau approximation, the level surface of mean radius a, is the same in both. 
But this implies that U,, outside the core is the same and therefore the equations 
of equilibrium of the shell hold; and so do the stress conditions. 

By using the condition that ‘ty is a function of py we can prove that the sum of 
(3) and (4) is exactly 

“cof eu. epg \?d¥ 

- JJ) (ross +438) Gee ) 
and in (2) the terms U, + u,0 ‘V’/dx;,, give the change of potential at a given particle 
in its new position. Hence where an equilibrium theory is applicable the usual 
theory of the figure of the Earth can be adapted to the present treatment. 
Examination of this theory showed that the terms in the ellipticity that contain 
ki/k? have coefficients of the order of 5. Hence neglect of these terms may 
lead to less inaccuracy than might have been expected. 

As an additional check the statical case has been treated by solution of 
Clairaut’s differential equation for a Roche model core. The elevation of a 
surface of constant density is taken as «K,/a, and if the central value of «€ is €, the 
result is 

€,=0°2604D, g(a)=0°604D, g(1)=0°319D. 
The results agree with (8) as closely as need be expected. 

Thus the whole of the equilibrium conditions are satisfied to the first order 
in the departures from homogeneity ; and any difference must be due to the neglect 
of square terms. Since corresponding terms have been retained in the treat- 
ment of the central particle model the solution for this is the more accurate. 
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3. Terms containing the velocities. —We first define the following core integrals. 





[| | po(x?-+ 2*)dr= Ay, [| [oa+y")\dr= Cy, { | { nar oe 

[J fro =A te o 
A, = ala eee (2) 
4n={{ | pot Zade=$ rai(hy Fh), (3) 
Aun [ || 08 5 dr= 5 mallhy— fh), (4) 


eee Ce 
C,-4,\= vd pute’) = ferent (hy — Wye ae? — $F aE} 5) 


In the second part of (5), the integral of the product by , is zero, and that of the 
product by &,€? is of order k7. Hence 


and similarly to the same accuracy 
Cy—-Ay=eAy, Cip—Ay2=eAe. (7) 
Also, if we neglect e, 
I a T= 3441p (8) 
{| | poze — Ay, | | | Po ~ —_ dee Aw (9) 


224 r2z4 222 
{| [ro qt Aj», i | po p arm igA 129 | | o— =e dr = Ars. (10) 


{{ polt,2dr = Ay(T? + 2) + 2(1 + 4¢)Ayy(lyly + ritytng) + (1+ 8e)A,o(i2+ m2), (11) 
| | 2poithit,,dr = 2A,(il, + nnn.) + 2(1 + $e)A, (lly + minting), fats 
{lf 2pqe{2(rnu; — Iu’,) —u.(ax — ly) }dr 

= 2wA,(til, — Im,) + 20A,,(1 + $¢)(tnl, — lms), (13) 
if 2pyco( ui it — ubit, dr — co 1 + €)[A,(l,2%, — mi) + Ay(1 > “) 


* (Lmhz + lym, — my, - msl;) +A,(1+% Fe)(1573 — msls)] +wA,(m,s,— 1,59) 
+ wAj, (mss, — 159) + 2w(A, — $A,,)(m,1,— im) + 2w(A,, — 3Aj9)( msl, — 1,4) 
+ $wAjy(tiny]s— Lyms) + $wAjo(tingl — Lyms), (14) 


where the usual: approximations have been made; the sign — indicates that a 
derivative with regard to t has been dropped. 
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There is evidently an advantage in adopting a new pair of variables 
(Uo, me) = (lym) + (1+ $6) (lym). (15) 
Then the terms in 2 that contain the velocities or C— A are 
= A(l? + m*) + w(2A — C)(lm — Im) - w(C — A) P+ m2) — 2(C — A)(le, + ey) 
+ Ay(E5 + 12) + 2A (lig + rinsing) + 2A (sil, — img) + (1 +8€)(Aqg— A3,/A,)(23 + 12) 
+ Ayo 1 + €)(Igtitg — mle) + $e A 4(Igritg — mel) + «(A yg — A?, Aya + Ze) 
x (Isms — msls) 
+ Ay (tings, — los) + 200(A, —5Aya)( tingly — Lom) — '$eo( Aya — Abs Ay)( tingly — iy) 
+ 52g (titgls — Igmms) + $e0( Aya — Af, /Ay)(tingls — |g). (16) 
Then in all 
=(1—H)(I2 + m®) + (1 — 2H)w( lin — im) — w* HP + m) — 2H(c,1 + cgm) 
+0°107[/2 + m2 + 2(Il, + ntitg) + 0°047630( 1 + 8e)(/2 + 22) + 200th], — long) 
+ (1+ €)cu( tingle — mel) +.0°5584.4e00(lgring — mel) 
+ 0°04763(1 + Y2e)(Igrtg — mals) + (rings, — Uys») 
+.0°04546:(ritgl, — J,m) — 0° 133400( msl, — i,m,) 
+ 0°83766:0(titgl;, — ],m;) + 0°05716:(tigl, — /,m;) 
+ w®[ — 425°9(r} +77) + 321°5(145) + 7252) + 77°78( dir, + dor) 
+ 1°179(r,/; + rgms) + 0°0638(r,/, + remy) + 29°05(djs, + dys.) 
— 108-2(s? + s3) — 1-888(s,/, + sym,) — 31-83(s,/; + Syms) 
— 1:210(/? + m2) + 11°98(/,/, + mym;) — 231°8(22 + m2) 
+0°845(d;)/, + dgls) + 19°69(d;/; + dgms) + 3°765(d," + d3”)) (17) 
115 Yo Sy» Soy Ly, my, ;, m; are quasi-statical coordinates. 
We find in turn 
(71, 72) = 0°3776(5;, 82) + 0°09131(d;, d;) + 0-001 4(/;, m;) + 0°0000075(/,, m4), 


(S15 $2) = 0°6147(d;, dy) — 0°3303(Is, me) — 0°01 96(/4, ms) mn 
+0°001126(m,/w, —/,/«), (19) 
(/;, ms) = 0°0001 196(m,/w, — I,/) +0:00001349(m3/w, — s/w) 
+ 0:00093(d;, d3) + 0°0278(/,, m,), (20) 
(1,.m,) = 0-002 10( tng /w, — Iy/e») — 0°00693(1H1/w, — I,/«2) — 0°139(d i, dy). (21) 


Forming the modified Lagrangian and eliminating dj, d, by 





(d;, d3) =o-006345(Sne2 + = ) +(i, m) 


w 
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we have 
md =(1-H+ 0:001018)(/2 + m?) + (1 — 2H + 0:002036) (lin — Im) 
~(H —0-001018)w*(/? + m?) — 2(H — 0-001018)(c,/ + com) 
0:002036 
+ 
147) 





(c,1n — cyl) +. 0°1070679( 2 + m2) — 0-0000289(Ie/ 5 + ritgtits) 
+ {0:005096(1 + 3e) + 0:0000489}(/? + 12) 

+0°2 1443(Il, + mm) + 0°21443(/,m — mg!) + 0°107(1 + e)w(lgmn, — mels, 
+ 0°0597 5ew(/,m, — Mel) + 0°005096u(1 + Ze)(1,%i, — msl) 


+ occ (cy1M3 — Coly) + 0°00001245(lr2, — mls). (23) 
We put 
l+im=(e™”, I,+img=C,e""', 1, +i1m,= Ce". (24) 

The equations of motion lead to three free speeds 0, —w, —w, as was to be 
expected. The other three are near these values, namely 

y=0'002531w, fg=—0'ggg940l, f,= +0-014¢ (25) 

y= —-w+n; n= —0'00403w, (,= —9°326f, (,= —18-4¢ (26) 

y=—-w+n; n=+0:00679w, f6=—9:333f, f,=+87°8¢. (27) 
The first makes the period of the Eulerian nutation 394-9 days, slightly longer 
than for the central particle model. ‘The slip at the boundary is essentially 
£,+€,=%,+2l,= —0-g9g0f (with a small correction due to ellipticity of the 
boundary). 

The second corresponds to the second free period for the central particle 
model. mn is nearly doubled. The third is a consequence of the variation 
of density in the core and is a completely new feature. {&, corresponds 
to ¢, of the central particle model and has nearly the same value. The striking 
result is that the slip in the second is much increased, and in the third is reversed. 
The cubic terms in the displacements are therefore of great importance. The 
term 0:0000489 (/}+ m3) in the modified Lagrangian has arisen mainly in the 
elimination of /,,m,, and the results are very sensitive to the difference between 
this coefficient and that of w(/3;%m3—mgl,), which otherwise is a multiple of e. 
Physically the explanation is that J,,m, includes a radial displacement, and 
consequently would imply a change of angular momentum, unless it was com- 
pensated by a change of transverse velocity. This can be verified by considering 
the behaviour of (17) when k,->0. All the J,,m, terms not containing the 
velocities tend to zero. ‘The matter tends toincompressibility and the coefficient 
of [+m tends to infinity. Thus /;,m;->0, and then the equations of motion 
for 1,,m, give 1,=0, m,=0. Thus [s,m are essentially a consequence of non- 
uniformity of density. 

The solutions for the principal forced nutations are in Table I, {, being the 
value for a rigid Earth in each case. The values are means of our solutions 
and a set made at the Mathematical Laboratory, Cambridge; the differences 
were within the possible effects of rounding-off errors, 
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TABLE I 
n/w C/Lo C/o l3/Lo Slip/Zo 
—1/13°7 1°067 —O'415 +O'274 —o'29 
—1/183 1°142 —I'312 — 2°03 — 1°89 
— 1/6800 099889 00104 0°047 0°023 
° I ° ° ° 
1/6800 1‘00120 —O°Ol12 —0°042 —0°023 
1/183 0°9707 0:258 -4°61 —1'o 
1/13°7 1°0266 —o'460 0482 — 0°32 


The slip is €g+(1 —A,,/A,)é3 or, nearly, C4 +23. 


We write /,+im,=(,, 1,+im,=¢,; the solution for the remaining variables 
is in Table II. 


TaBLe II 
niw Qa)/D g@®/D %4/D {,/D h k lL 1-—h+k 1+h—#k 
- — 0°324 0°617. —0'126 —0'002 0°'597 0'258 0°070 0°661 1‘210 
=" 0°385 0°779  —1°433 0°002 «0710: 0298 0°'072 0°588 1263 
o 0'299 o°552 +0°799 o'00! O°551 0°244 0°082 0°693 1°185 
* 0308 0°574 +2°922 0°027. =0568 ~0°264 0°084 0:696 1°172 
= 0°327 0625 —0'169 —0'003 0°603 0:261 0:078 0°658 I°211 


The values of 1—~h+k and 1+h-— $k for semidiurnal and long-period tides are 
0675 and 1°188. On the whole the behaviour is very similar to that for the 
central particle model. The chief differences are an appreciable increase of 
yielding for all the K, tides and a considerable change in some of the variables 
for the P tides. The latter are connected with the change of the free 
period near n/w=—1/183 and the introduction of a new free period near 
n/w= +1/183. The theoretical amplitude of the semiannual nutation in 
obliquity is 0”-55; if this could be measured it would give useful information. 

4. The Roche model used certainly overestimates the compressibility, and 
better theoretical values will usually be found by interpolating; if x denotes an 
estimate for the central particle model and y that for the Roche model, 0-4. + o-6y 
should be better than either. This will not apply, however, to the semiannual 
nutations and P tides, since the effects of changes of the free periods 
cannot be interpolated in any such simple way. 

The correction to the 1g-yearly nutation is reduced, and the correcting 
factor should be close to —0-0022. ‘The coefficient would be about g”-207, in 
fair agreement with most observed values, though the difficulty about the ratio 
of the amplitudes in obliquity and longitude persists. 

Both models are imperfect in so far as they do not allow for the, size of the 
inner core. Allowance for this would introduce much greater difficulties. For 
a homogeneous incompressible core, and with some accuracy for a core deviating 
somewhat from this condition, the solution depends on the differential equation 


GQ FQ PQ 
a 4 aa escinnie a 
y (Fa + = +(y?—4w a3 =O. (28) 
For w=o solutions can be built up from 1/r and its derivatives, and have no 
singularity except at r=o, But when |y|<2w this procedure leads to functions 
12* 
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with infinities on the cones z= + (4w?/y? — 1)"?(x? + y®)"?, and these would extend 
right through the fluid. The discovery of non-singular solutions between two 
spheres would resemble the treatment of oscillations in an elliptic lake; they will 
be of exponential type over some ranges of latitude and oscillatory over others. 
Whereas the effects of an inner sphere, in potential problems, are of the order of 
the fifth power of the ratio of the radii, they may here be of the second or third. 
Since the radius of the inner core is 0:36 of that of the main core the effect is 
probably only of the order of a tenth of that of the core as a whole. 

5. The effect of the ocean.—Larmor (1915) and one of us (1915) noticed at 
almost the same time that the effect of the ocean on the period of the Eulerian 
nutation is not negligible.* It could be allowed for in the present problem as 
follows. The period being long, a statical treatment is assumed permissible. 
The elevation of the ocean surface can be introduced as an additional coordinate, 
contributing to the potential energy, and a new numerical solution can he made. 
However, this hardly seems necessary. The effect of the fluidity of the core by 
itself would be to reduce the period by about 30 days, and the increase due to the 
elasticity of the shell is about 120 days. Thus elasticity is much the more im- 
portant. But for an entirely elastic body the effect on the period depends wholly 
on Love’s number k, and we can use any model that reproduces this correctly. 
We therefore consider a homogeneous incompressible Earth of rigidity k, covered 
by a shailow ocean. The theory for the solid is well known; for the ocean we 
add the conditions that the free surface is an equipotential and that the normal 
stress is continuous at the interface. We denote the normal displacements at 
the interface and the ocean surface by 9(1)K,/a, q(1+)K,/a, the densities by p 
and p(1 +7), and write 

° : be a 
U,=hk,K3; zi +)e =, 
We are led to the equations 


. qiit) _ ake, 
¥ + tate yey 2 


3 ak 
ay at) eae t)=- 
poe 2 &_ Wit) + 78(1) 
~ Sak, 1+ ’ 
Here 1+7=5'53. For the rigid Earth, o= 0, g(1)=0, and the solution is 
q(t +)=1°1215 ak,/g, 
whence k=0-1215 and the free period given by the Larmor-Love formula is 350 
days. ‘Thus in this case the ocean lengthens the period by about 45 days. 

We have found from our theory, with no ocean, k=0-289. A uniform solid 
Earth would give k= 0-290 if = 0-439, and the free period would be 441 days. But 
it is better to proceed as follows. ‘The theory, with no ocean, leads to a free 
period of 392 days. This would correspond, inan entirely elastic body, to k = 0-208, 
and the corresponding rigidity would make o=0-652. We insert this in the 
above equations and get 

q(1)=0°3266, g(1+)=1°3009; k=0-301 
and the free period is 449 days. Thus the lengthening due to the ocean is 57 days. 


j * The account of Jeffreys aimed at an estimate of the damping due to viscosity in the ocean; 
this has of course long been superseded. The lengthening of the free period came out as a by-product. 
Larmor used an equilibrium theory, which cannot be far wrong for a 14-monthly tide. 





and 

















sta rtesasallll 


Pai ie Ae th alceb nlc Rt: es 


rs i ter cet Ai 


No. 2, 1957 Theory of nutation and the variation of latitude 171 


Larmor made various rough corrections, allowing for the fact that the ocean 
does not cover the whole surface and for elasticity. The latter has already been 
allowed for here. The oceanic effect is likely to be about 3 of that for a complete 
ocean, namely about 38 days, and the free period will be about 430 days. The 
uncertainty is probably about +5 days. ‘This is sufficiently close to the observed 
value. 

The forced nutations correspond to diurnal tides. There is no satisfactory 
theory for these in the actual ocean. For an ocean covering the Earth, and for 
exactly diurnal tides (speed w) there is no change of level at all, but there are tidal 
currents. ‘These may contribute to the angular momenta. We may note, how- 
ever, that the effect on the nutations is likely to be of the order of the ratio of the 
moments of inertia of the ocean and the body of the Earth, say 1/1000. When 
the precession is used to compute (C—A)/C an error of this order may be 
introduced. But for the nutations this factor is likely to be the same within a 
factor of order n/w. Hence if the computed value of (C—A)/C is used the 
correcting factor due to the oceans for the nutations will be of order 4; n/w 
and utterly negligible. 

Attention needs to be called to the fact that ¢, does not correspond to the 
nutations given in the standard accounts. Most of these proceed by writing down 
the equations of motion and neglecting second derivatives with regard to the time. 
In our definition of ¢, these derivatives have been taken into account; thus the 
standard treatments do not quite give the theoretical motion even for a rigid 
Earth. 

Oppolzer (1886), followed by Woolard, defines the nutation as the motion of 
the instantaneous axis of rotation. ‘This motion is almost the same as is given by 
the approximation, but then the motion of the axis of figure relative to that of 
rotation has to be added as a correction when we are considering displacements of 
stars relative to the Earth. Woolard (1953) gives the coefficients of the fortnightly 
terms in obliquity and longitude, for the axis of rotation, as 0”:0884 and — 0”+2037. 
The corresponding values for the axis of figure are 00948 and —o”:2211. 

It is worth while to have some indication of how much of the difference from 
the rigid body theory of the nutations is due to elasticity of the shell and how 
much to fluidity of the core. This can be obtained from (24) of our previous 
paper by omitting slip at the core boundary, thus taking ,=0. ‘The free period 
becomes 440 days. For the nutations we should have 


C=(1—0°3045 m/w) Ly. 


The effect of elasticity is therefore negligible for the 1g-yearly nutation, as 
Poincaré noticed. ‘The extreme value for the periods considered here are: 


n/w £/Lo 
I 
—_— 1°022 
13°7 
I 
—— 0978 
13-7 97 


For both components the effect of the core is the remainder, about 0-05 in each 
case. 
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6. The calculated value of the 19-yearly nutation.—Spencer Jones’s data on 
the solar parallax and the lunar inequality led to (Jeffreys, 1948) 


E/M =81:278 +0°025, H=0-00327360 + 000000069. 
For a rigid Earth these give for the constant of nutation 
N=9":2271, 


in good agreement with Spencer Jones’s 9”:2272. (It should be noticed, however, 
that this depends on using a steady-motion theory ; allowance for the nH term in 
7(23) of our previous paper would decrease this by about 0”-oo10, and there would 
be a proportionally larger decrease in the longitude term. We shall however 
neglect this, since the effects that concern us are in the second decimal.) Rabe’s 
parallax (1950) and the lunar inequality from Spencer Jones’s data give 


E/M =81-356; H=0-00327468; N=g9":2242. 


The change of the solar parallax produces more effect in H than in N. 
Applying the corrections ¥0-0022 to ¢/f) for n/w= + 1/6800 we get for the 
nutation components in obliquity and longitude: 


a hee _ Rabe’s parallax 
Jones’s parallax 
Obliquity g"*2120 g”-2090 
Longitude 17"*210 17"*204 


Multiplying the longitude terms by sine to give actual displacements gives 
6"-850 and 6”-848. ‘The obliquity values by themselves would be in agreement 
with nearly all observational estimates; but Fedorov’s discussion, which is so 
far the only one that separates the corrections in obliquity and longitude, gives 
the values 6”:886, 6":874, 6":874 for longitude term from three series of obser- 
vations. ‘Thus the longitude term remains in disagreement. 

7. Takeuchi’s model for the shell is close to Bullen’s Model B; and it may 
be asked whether permissible changes in the shell model would have appreciable 
effects, in particular because if the temperature distribution is not adiabatic 
there will be corresponding changes in the distribution of density and propor- 
tional changes in the elastic moduli. This appears to be answered by a recent 
paper of Bullen (1956), who shows that for his Model A the greatest possible 
correction to the density at any point is only 0-07 gm/cm*. The effects of 
elasticity of the shell are therefore reliable to about 1 per cent. 

8. It has been suggested that electromagnetic interaction between the core 
and the shell may have an appreciable effect. ‘This is not considered here, since 
the object of this work is to test whether the mechanical properties inferred on 
other evidence suffice to explain the 14-monthly period and the nutations quanti- 
tatively, and thus to provide a comprehensive check on the theory of the mech- 
anical properties. If a definite discrepancy had been found it would have to be 
considered as a starting-point for further investigations; but as the agreement 
seems to be within the uncertainties the presumption is that magnetic effects 
are small. 


St. John’s College, Faculty of Sciences, 


Cambridge: The University, 
1956 December. Lisbon. 
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THE VARIATION OF DECIMETRE-WAVE RADIATION 
WITH SOLAR ACTIVITY 


C. W. Allen 


(Received 1957 January 19) 


Summary 


A statistical method is used to segregate the quiet component from the 
slowly varying component of solar decimetre-wave radiation in the period 
1947-54. For this purpose the radiations at frequencies 2800, 1200, and 
600 Mc/s have been correlated with sunspot numbers, sunspot areas, and 
faculae. The mean lives of the various radiations and activities have been 
estimated and compared. ‘There is an increase of life with decreasing radio 
frequency. The life of 2800 Mc/s radiation is about the same as sunspots 
but measurements of the latter show some anomalies. The slowly varying 
radiation flux is proportional to frequency in the range studied. ‘The quiet 
sun flux has a small but significant variation with solar activity, the relative 
change being greater for smaller frequencies. The possibility that this 
variation may be associated with uncorrelated local sources, such as promin- 
ences, is not entirely excluded. 





1. Introduction.—The radiation from the Sun in certain wave-length bands 
is clearly composed of (i) a quiet component, and (ii) a slowly varying component 
(x). This is the case, for example, with decimetre radio waves and with the 
ultra-violet radiations that produce the ionosphere. ‘The regular observations 
give a measure of the total radiation only and thus call for a method of segregation 
in order that the two components may be studied individually. In the present 
paper an attempt is made to segregate the components of decimetre radiation 
during the period 1947-1954. 

The quiet component comes steadily from the whole solar surface and the 
main question concerning it is whether and how it varies during the sunspot 
cycle. It is thermal radiation and any variations would represent interesting 
changes in the normal solar atmosphere. The slowly varying component comes 
from the local active areas which can be detected by the various visual activity 
phenomena of sunspots, faculae, etc. We can recognize and measure this 
component by its correlation with the activity. We are interested in its intensity 
in relation to the various measures of solar activity, and also in the question of 
whether the activity and the radiation appear over the same period of time. 
We are not concerned with the more rapidly varying features—noise storms, 
bursts and outbursts. 

It would be possible in principle to detect and measure all the local emission 
sources by the use of eclipses or high resolution interferometry. This is difficult 
to do comprehensively at sunspot maximum, but it is hoped that in due course 
the conclusions of this paper will be checked by adequate observations of this 
type. 
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Our method of segregation was first used by Pawsey and Yabsley (2). Later 
it was used by Christiansen and Hindman (3) to demonstrate sunspot cycle 
changes in decimetre radiation. The method can be explained by reference to 
Fig. 1, where values of decimetre radiation flux are plotted against simultaneous 
measures of solar activity. Such plots can be made for any phase of the sunspot 
(sp-) cycle and extrapolating the resulting curve to zero activity can give the 














Decimetre radiation flux —> 











Measure of solar activity —» 


Fic. 1. 


quiet component. However we must allow for the following complications of 
which (g) has been regarded as the most important (4, §). 
(a) There are two regressions representing the linear solution for the points 
and the true solution may be anywhere between them. 
(6) ‘The curves may not be linear. 
(c) There may be changes in the measurement standards during the sp-cycle. 
(d) The scatter of points may be due in part to the progression of the sp-cycle. 
This would make the slope steeper as can be seen by the broken line in 
Fig. 1 representing the slope for the whole sp-cycle. 
(e) The slope might be affected by measurement changes associated with 
directivity and the position on the Sun’s disk. 
(f) Some local emissions might be entirely uncorrelated with measurable solar 
activity. 
(g) The local emissions and the corresponding solar activity may not occur 
throughout the same period of time. 
‘These complications were taken into consideration when selecting the data 
and the form of the correlations. 
2. Data and procedure.—The radio data used were: 
28 The daily flux at 2800 Mc/s as published by Covington (6,7) and 
multiplied up to the end of 1953 by 1°18 in accordance with (7, No. 104). 
12 The daily flux at 1200 Mc/s measured at the Sydney Radiophysics 
Laboratory (7, 8) and corrected for 1947 by the factor 1-9 (9). 
6 ‘The daily flux at 600 Mc/s also from the Sydney Radiophysics Laboratory 
(7, 8) and corrected for 1947 and 1949 by the factor 1-2 (9). 
These are the only decimetre-wave data available over an appreciable part of 
the sunspot cycle. Monthly means have been derived. 
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The solar activity data used were: 

R_ Daily and monthly Ziirich sunspot numbers (7, 10). 

A Daily and monthly corrected sunspot areas averaged from Washington (11) 
and Greenwich (12, 13, 14). 

F Daily and monthly corrected facular areas from Greenwich (12, 13, 14). 
Since faculae are seen only near the limb, the daily values were converted 
to the Sun’s centre by adding the eastern faculae of 5 days earlier to the 
western faculae of 5 days later. Data were supplied by the Greenwich 
Solar Department for this purpose. 


The reason for using three different measures of solar activity was that the 
complications (a) to (g) take on a somewhat different form in each case while the 
deduced segregation of the quiet and slowly varying radiation should remain the 
same. The repetition of the analysis for three different activity variables provides 
a much needed check on the results. 

As far as possible the six sets of data, which we will label 28, 12, 6, R, A, F, 
were treated alike. 

In order to examine the material the values were plotted out on daily and 
monthly charts. ‘The 1951 daily chart is shown in Fig. 2, and the monthly chart 
in Fig. 3. The monthly means were subject to a 13-month smoothing (with 
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Fic. 2.—Daily variations during 1951. 


half weight for the extreme months as adopted for Ziirich sunspot number 

smoothing) and the values correlated were the departures from the smoothed 
curve. 

In order to minimize the regression ambiguity (a) it was desirable to have 

a good correlation and a wide range. Monthly values were considered, but 

_ their range is too small for these to be suitable for the main correlation. Daily 

values are better but they require some averaging to cut down the scatter and 

to keep the numerical work within reasonable bounds. ‘The daily plots (Fig. 2) 

show that the maxima and minima (which are to supply the key points for the 
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Fic. 3.—Monthly variations. 

correlation) last about 7 days. Consequently weekly means were used for the 
main data. We could not afford to miss the main maxima and minima by taking 
random weeks and instead three weeks per month were carefully selected with 
the view to giving the most useful values for correlation—usually one maximum, 
one minimum, and one plateau in each month. ‘The central days for these 
7-day periods (Fig. 2) were chosen by considering together the daily charts 
for sunspot area, sunspot number, and 2800 Mc/s radiation. Care was taken 
to ensure that the sunspot and radiation data were given equal consideration in 
selecting the days. ‘The weekly values for the correlation were scaled from the 
daily charts by using a mask to obscure all but the 7 days concerned and then 
measuring the mean distance of the daily curve from the 13-month smoothed 
curve. 

The departures from the smoothed curve so found are designated S(28), 
S(z2), S(6), X(R), X(A), and X(F). They are zero in the mean and therefore 
in a convenient form for correlating. 

The various values of S were plotted and correlated with the simultaneous 
values of X. The main requirement is the slope dS/dX, since, as can be seen 
from Fig. 1, the mean quiet radiation in any year is obtained from the mean 
total radiation minus the product of the slope and the mean activity. The 
procedure has been to make a primary estimate of the slope and then allow for 
the corrections to it later. The primary estimate was made both by plotting the 
values and by considering the two regression slopes [S?]/[SX] and [SX]/[X?], 
where the square brackets signify summation. ‘The study of the plotted values 
was given most weight since in this way one could discard some badly correlated 
points on the positive side of the diagram. One can see from Fig. 1 that it is 
the negative values of X that are to be extrapolated for determining the quiet 
radiation. When using the regressions the slope estimate was placed between 
the regressions and nearer to the regression on that variable thought to be more 
accurately represented. Details of this are given in Section 3. Both weekly 
and monthly values were considered for estimating the slope. The various 
estimates were weighted as thought appropriate and a mean slope (usually for 
each year) obtained. Results are given in ‘Table I. 
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TABLE I 

Primary estimates of slope dS/dX 

Unit of S =10-* wm-? (c/s)~} 

Unit of X=1 sunspot number (for R) 

Unit of X=1 millionth of hemisphere (for A and F) 

28 12 6 
Period on on on 
R A F R A F R A F 

1947 0°87 0°033 0°053 | 0°30 o'021 O'l4 O'o10 
1948 0°73 0'037 0°052 
1949 o'72 0°033 0060 | O'19 O’oll 0°07) =: 0003 
1950 0°66 0°031 0048 | 0°20 0°009 0708 =: © 004. 
1951 0°75 0030 0°046 | 0°18 0°005 0°07 = 0003 
1952 0°50 0'033 0'037 | 0-16 0°007 009 = 0004 
1953 0°44 0°030 0°030 | 0°08 0°004 o'1Ir 0006 
1954 o'4I 0°026 0°025 | 0°09 o'10 
1947-49 | 0°77 0°0343 0°056 | 0°25 O°OI4 O'014 | O11T =0°006 = 0007 
1950-51 | 0°72 0'0306 = 0046 _| 0°19 0°007 O'OII | 0708 0°004 0°007 
1952-54 | 0°48 0°0320 §=0035_ || 0°13 0:006 o-oI12 | 0°09 =0°005_~=— 0006 
1947-54 | 0°74 0°0326 = o'051_ | O'21 0'008 0'012 | 0°09 0°0045 0°007 




















The correlation coefficients r =[SX]/([S?][X?])'? are given 


TaBLe II 


Correlation coefficients r for weekly values of S and X 


in Table II. 








28 12 6 
Period with with with 
R A F R F R A F 
1947-49 | 0°88 0°93 0°62 082 o69 ~=§=60°42 068 of0 8 062 
1950-51 0°94 0'94 o'7I 0°62 0°60 0°50 0°46 0°45 0°40 
1952-54 0°77 0°83 o’61 0°56 «©6058 )~—0 49 0°37 0°38 0°28 
1947-54 | 0°88 0°92 0°63 | 0°69) «©6062, 046 | O51 = 4544 























3. Complications and corrections.—The complications will now be considered 
in the order mentioned in Section 1. Where possible, correction factors are 
derived which must be multiplied by the primary slope estimate before it can be 
used as in Fig. 1 for finding the quiet radiation. The factors are collected in 
Table V. 

(a) Most of the scatter in the S, X correlations is associated with the imperfect 
relation between the radiation and the measures of solar activity. For the 


correlation of S(28) with X(R) and X(A) the scatter is small (‘Table I1) and it 
is reasonable to adopt a primary slope (which we later identify with m;) midway 
between the regressions, namely ([.S?]/[X*])"*. This implies that the random 
effects causing the scatter are associated equally with the radiation and the 
activity. The error factor from this assumption must be less than 1—r and is 
not very serious. Nevertheless there is no way of overcoming this residual 
ambiguity and it was for this reason that the data were selected in such a way 
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as to correlate as well as possible. When S(28) is correlated with X(F) the 
scatter is much greater and we may blame the extra scatter onto the random 
variations in X(F). Hence for S(28)/X(F) it is appropriate to adopt a primary 
slope close to the regression on S, i.e. [S?]/[SX]. By a similar argument, since 
S(z2) and S(6) correlate rather poorly with X(R) and X(A), it is evident that the 
blame this time lies with S and we adapt the slope close to the regression on_X. 
The correlation of S(z2) and S(6) with X(F) is also poor. It is thought that 
the random scatter is associated rather more with S than X in this case and the 
solution is taken somewhat closer to the regression on X. Not much reliance 
can be put on the S(z2, 6)/X(F) slopes. 

(6) In order to determine whether the solution should be curved, one may 
compare the slopes at different stages of the sunspot cycle. The slopes from 
Table I have been plotted against activity in Fig. 4 and smoothed relations drawn. 
In general the slope is steeper for greater activity. The correction factor (d) 
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resulting from the curvature of S on X is simply the area under the slope curve 
(Fig. 4) to the left of the activity abscissa divided by the product of the ordinate 
and the abscissa. 

It should be mentioned that one would not expect a severe departure from 
linearity. The characteristics of sunspot groups do not change noticeably with 
the sp-cycle and the radiations from the groups must be additive. Therefore 
one would expect both the radiation and the activity to be proportional to the 
number of spot groups. 

(c) Changes in the measurement constants could have a serious effect on the 
present analysis which depends a great deal on differences from one year to 
the next. ‘There is no reason to mistrust the constants of the activity data which 
have been given in the same form for many years. The correlations lead to the 
impression that the 2800 Mc/s data is consistent throughout the period. ‘The 
correlations of 1200 and 600 Mc/s data with activity are not very good but one 
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cannot conclude that this is due to undetected changes in the measurement 
constants. Fortunately the change of this radiation with the sp-cycle is so great 
that the main conclusions are not likely to be affected by errors of measurement. 

(d) The variations due purely to the sp-cycle phase change are effectively 
eliminated by making the measurements from the 13-month smoothed curve. 

(e) Suppose Fig. 5 represents a characteristic difference in the daily curves 
of X and S such as might be due to the radiation having a higher directivity 
from the Sun’s centre than the activity. The slope dS/dX would be affected if 








Activity 
bod 


Pie.” P anal 
re 





Radiation 








Tew ewww Se eee wee ee eee eaewewewewewe Sew! 


Days 





Fic. 5 


peak days were chosen (as we have done) whereas if random days had been 
chosen any peaking effect would have averaged out. A study of the daily charts 
does not reveal any such characteristic differences. To test this, for each selected 
peak and dip of S(28,12,6) and X(R, A, F) I have measured the value on the 
central day and also the mean of the values three days earlier and later. Averaging 
the results gave the ratios of Table II]. There does not appear to be any 
significant change with the sp-cycle, nor do the peaks appear to differ systematically 


Taste IIT 


Ratio of departure from the smoothed curve three days before and after the selected day to 
the value on the selected day 








’ | | 
| 
28 12 6 | R A F All | 
Peaks | | 
1947-49 0°62 o62 0°54 0°62 0:63 0°49 0'60 
1950-51 0°63 (0'55 «=O | 0°58 o85 1:°0 082 
1952-54 0°62 076 I°0 | 0°67 0°68 0°53 0°65 
1947-54 062 o'60 0°98 061 069 0°70 0°66 
Dips | 
1947-49 0°62 082 o'50 | 0°55 O59 0°67 o'60 
1950-51 o-71 068 0°73 060 067 0°77 0°69 
1952-54 o80 1'oO 1'0 o-71 o61 o'50 0°67 
1947-54 0°67 0°79 0°70 0°59 0°62 0°67 0°63 
Mean 0°65 0°73 0°79 060 0°64 0°68 0°64 
Increase 1°06 1°00) 0°95 I'IO.) «61°06 104 1°06 
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from the dips. To obtain a mean, double weight was given to the dips since 
these are more important to the analysis. The smaller this ratio the more the 
measurement has been increased by the policy of selecting the peaks. The 
amount of this increase, given in the last line of Table III, is obtained by 
considering typical curves such as shown in Fig. 5. The correction (e) to allow 
for it is the (activity)/(radiation) ratio of these increases. 

(f) See Section 6. 

(g) The effect of any time difference in the radio emission and the optically 
observable activity is one of the main interests of this work and is discussed in 
Section 4 where a formula for the correction factor (g) is derived. 

4. Persistence and delay of activity or radiation.—The investigations of 
Piddington and Davies (4), of Waldmeier (5) and of Dodson (15) suggest that 
the radiations from sunspot areas disappear less rapidly than the sunspots 
themselves, and these authors consider the delayed radiation from old sunspots 
to be the main cause of variation in the non-sunspot component. However 
Piddington and Davies’ composite areas (4, p. 588) imply a delay of radiation 
variation of at least half a month whereas a study of Fig. 3 for 2800 Mc/s radiation 
reveals almost no delay. Piddington and Davies gave no method of determining 
the delay quantitatively but showed the importance of controlling the delayed 
data by comparing it with advanced data. Waldmeier made a quantitative 
estimate of radiation delayed by one rotation but did not prove the reality of the 
effect by comparing it with radiation advanced by one rotation. Had he done 
so he would probably have found that much of his effect was not due to delay. 
The present work is rather similar to these in aim and approach. However the 
factors determined are quantitative, and allowance is made for the persistence 
of both the activity and the radiation, which means that advanced and delayed 
data are compared on the same basis. 

The analysis depends on the correlations of S(28, 12,6) and X(R, A, F) with 
one another, both for simultaneous values and for values separated in time by 
a small number of 27-day rotations. We will use the notation X_,, X_,,... to 
represent the value of X at 1, 2,... rotations earlier than the occasion being 
considered. S and X are taken as the departures of the radiation or activity 
from the 13-month smoothed mean, and hence the summations such as [S], [X], 
and their uncorrelated products, will be approximately zero. 

The amount of solar activity X at a particular time will be made up of new 
activity Y plus the remains of activity that was new I, 2, ... rotations earlier. 
We can then write 

X_,=Y_,+aY_.+6Y_,+ 


X =Y +aY.,+6Y_,+ ; (1) 
Xi3=Y.,+aY +bY_,+ 
where a, b, c, ... are the factors by which the activity will be decreased in 
I, 2, 3, --. rotations. We will assume that there is no correlation between 
.Y., Y, Y.4, Ys... . Multiplying and adding (e.g. the weekly values 
throughout a year) gives 
[X7] = [Y?](1 + a®+ 6? +...) 
[XX_,] =[XX.,.]=[V*(a+ab+ be+...) >, (2) 


[XX_,]=[XX,.] =[V2(0+ ac+bd+...) 
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and the low order auto-correlation coefficients for adjacent and nearby rotation 
months become 


[XX_,]/[X2]=(at abt bc+...)(r+@+b +...) ; 


[XX g]/[X2]=(b+ actbd+...)(r+a2+02+...) (3) 


Now take the radiation S emitted from the activity Y to be mY, and the 
factor by which S varies in subsequent rotations to be a, B, y, ..., then 

S=m(Y+aY_,+BY_.+...). (4) 

The S auto-correlations are similar to (3) with a, 8, ... replacing a, 6,.... By 
multiplying the items of (4) with (1) and adding 

[SX] =m[Y?](1+ax+5B+... 

[SX_,] =m[Y?](a+aB+ by+.. 

[SX_.] =m[Y?](B+ay+ b5+.. 


. . . 
ee 
‘ 
—-— 
. 





(5) 
[SX ,,]=m[Y?](a+ ab+ Be+.. 
[SX ,.]=m[Y?](b+ ac+Bd+... 
The regression slopes become 
S on X=[SX]/[X*] = m(1 +ax+bB+...)/(1+ a? + b? 4+...) (6) 
X on S= [S?]/[SX]=m(1 + 2? + BF +...)/(1+ax+bB+...) |? 


and the mean of these is 
my = ([.S?]/[X?])!? = m(1 + a? + B2+...)22/(1 +a? + 527 4+...)"2, (7) 


It must be noticed that the scatter entering these statistics is associated with the 
non-correlation of Y with Y_,, Y_s,... only. Other random effects such as 
errors of measurements and variability of m are not included. However mr, can 
be identified with the primary estimate of slope mentioned in Section 2. 

The summations of (5) are mainly used in the form of ratios such as 


[SX_,]/[SX]=(a+aB+by+...)/(1+ax+bB+...), ete., (8) 


which are numerical ratios of various regression slopes and eliminate m and [Y°]. 
We now consider our main purpose which is to find the radiation for zero 
activity, i.e. to find S when X= —M, M being the smoothed sunspot activity. 


For a completely quiet Sun Y, Y_,, Y_., ... must all represent zero activity and 
will all have the same negative value Y_y. ‘Then from (1) 


Y_y=—-M/(1+a+b+...). 
The corresponding zero activity value of S is, from (4) and (7), 
S_y= —Mm(1+a+B+...)/(1+a+b+...) 
= —Mm,(g), (9) 


where 


_ 1+a+B+.. (1+a?+b?+...)' 
(g)= It+at+b+.. (1+07+f?+...)!/2" (10) 





Thus (g) is the factor by which the primary slope estimate m, must be multiplied 
to allow for delayed action, 
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Since the higher members of the series 1, a, b,..., 1, %, B,... cannot be 
determined from the correlations with adequate accuracy it is necessary to seek 
simpler series to replace them. The most satisfactory solution is 

I, a, b,c, ...>1, a, a, a, — 
1, & By yy... 1, a 2 08... f” (11) 
The relevant formulae (3) and (8) then become 
[XX_]/[A"] = a =[SX,4)/[SX] 
[XX_4]/[X7] = a? = [SX , 9]/[SX] 





[SS_,]/[S2]= « =[SX_,]/[SX] [ (12) 
[SS_4]/ [S]=22=[SX_4]/[SX] | 
and (10) becomes 
t= (F=2.8)". 9 


We see from (12) that both the auto-correlations (left side) and the regression 
slope ratios (right side) give values of a and «. 

The values of [XX_,], [X*], [SX_,], etc., obtained from the data are not 
precisely the same as the summations required for (3), (8) and (12) since the 
former are increased by some random scatter. ‘This extra scatter will not 
appreciably affect the ratios in (8) and the right side of (12) since they can be 
regarded as ratios of regressions on the same variable. However in (3) and the 
left side of (12) the [X*] and [.S?] from the data will be artificially too high and 
we must be prepared to increase the ratios by a few per cent. 

It is important also to allow for the fact that the selection of weeks was not 
a random one. For this reason [X*] was somewhat greater than would be 
obtained from a normal dispersion and greater than [X2,] or [X?,]. An allowance 
for such differences has been made by making all the correlations and ratios 
dimensionless, e.g. using [YX_,]/({X?][X2,])!” on the left of (3) or in (12), and 
[SX_,][X?]#?/[SX]LX= ,]'? on the left of (8) or in (12). 

5. Results.—The results that might be expected to come from the analysis 
are (i) the decay factors represented by a and «, (ii) the flux of slowly varying 
radiation, and (iii) the flux of quiet radiation. 

(i) The decay rates.—The values of a(R), a(A) and «(28) obtained from the 
various parts of (12) are in reasonable agreement and suggest that (11) is applicable 
for R, A, and 28. For F, 12, and 6 the decay is slower and disagreements in (12) 
show that (11) is not entirely adequate. Nevertheless the data are not consistent 
enough to justify the attempt to determine more than one constant for each 
phenomenon and the results have been averaged from (12) as well as possible. 
Monthly data have been analysed for this purpose as well as weekly, and the 
decay constants finally adopted are given in Table 1V. ‘The values of a and « 
refer to 27-day periods. For comparison with other data the decays have also 
been converted into mean-life in calendar months. 

The variations in the mean-life of solar phenomena have been demonstrated 
in the past (16, 17) by correlation delays. It can be seen from (12) that the 
correlation delay will be approximately half the difference of mean-life for the 
two phenomena concerned. Values are given in Table IV. ‘The delays with 
respect to R have been determined also by the earlier method (16) and these too 
are given in Table IV. The comparison is reasonably satisfactory. 


a3 
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Tasie IV 
Values are averaged to represent the whole period 1947-54 





R A F 28 12 6 
a or « (27-day) 0°59 0°53 0°72 0°59 0°67 0°70 
Mean life (in months) 1°7 1°4 2°7 1°7 2°2 2°5 
Half mean-life difference 0°0 —O'ls +050 roe) +O'25 +-o40 
Delay from R (in months) ole) —0O'4 +06 —o'2 orl Log 

















Points of interest in Table IV are as follows. There is an increase of mean-life 
and delay as the frequency is decreased. At 2800 Mc/s there is no significant 
delay as compared with sunspot activity. The longer life and increased delay 
of faculae in comparison with sunspots has long been known (16). ‘The difference 
between R and A is quite surprising. It is apparently due to an effect whereby 
the number of spot groups on the Sun tends to increase one or two rotations 
after, and in the neighbourhood of, an enormous sunspot. ‘This result appears 
clearly for the six longest sunspots ever observed (14), half of which come within 
the period studied. 


(ii) Slowly varying radiation.—The various yearly means and factors from 
which the slowly varying radiation has been determined are assembled in Table V. 


The units are Ziirich numbers for R, millionths of the visible hemisphere for A 


and F, and 10-* wm-* (c/s)-! for radiation; corresponding units are used for 
the slopes. The correction factors (b), (e), and (g) have been explained in 
Sections 3 and 4. A certain amount of smoothing of the data has been adopted ; 
in many cases the value quoted is the mean of the individual and the smoothed 
value. The slowly varying component of radiation in the last column is the 
product of the five columns preceding it. 

The interagreement of the values in the last column obtained almost 
independently from R, A, and F should give an indication of the reliability of 
the procedures used. Table Vi gives the weighted mean of the same values—the 
probable errors should be about 1 or 2 units in this table. The fluxes associated 
with various measures of activity (each representing approximately a normal 
sunspot maximum) are given in Table VII. The slowly varying flux is pro- 
portional to sunspot area but the relation with sunspot number and faculae is 
not quite linear. 

The flux is seen to be iispreeienle proportional to the radiation frequency 
in the range studied. Hgwever this relation does not persist for frequencies 
above 2800 Mc/s. At 3750 Mc/s Nagoya observations (7) show rather less flux 
for a specified sunspot area than at 2800 Mc/s, and at g400 Mc/s Sydney 
observations (7) show barely detectable slowly varying radiation—perhaps about 
1/10 the flux at 2800 Mc/s. These points have been clarified by a graph plotted 
out by S. F. Smerd and-kindly sent to the author. 


(iii) Quiet Sun radiajion. —The flux from the quiet Sun remains after the 
slowly varying component has been substracted from the total. The yearly 
values are shown in T able VI and the values associated with defined measures 


of activity are given in Table VII in terms of both flux and apparent solar 
temperature, 





Pere ee er Tt pswad, 











Certs tt 


toi ear ene 


ca es Nea Bie 





No. 2, 1957 


Radiation 
Activity 


28 from R 


28 from A 


28 from F 


12 from R 


12 from A 


12 from F 


6 from R 


6 from A 


6 from F 


Year 


53 
54 


of decimetre-wave radiation with solar activity 


ot bl] 


152 

136 

135 
8 


32 
14 


2640 
2010 
2160 
1300 
1170 

430 

150 


2890 
2330 
2600 
1750 
1380 
710 
330 
140 


162 
133 


69 
32 
14 


2320 
2070 
1300 
1170 
430 
150 
40 
2940 
2650 
1750 
1380 
710 
330 
140 


162 
133 


69 
32 
14 


2320 
2070 
1300 
1170 
43° 
15° 
40° 
2940 
2650 
1750 
1380 
710 
33° 
140 


Slope 


087 
0-78 
o°77 
0°67 
0°68 
050 
0°44 
o'41 
00336 
00350 
0°0332 
0°0314 
00312 
0'0317 
00302 
0°0281 
0055 
0052 
0°057 
0049 
0°046 
0°037 
0030 
0-025 
0:28 
o'2I 
0°20 


8888S SoSkeei 


600000000 00900000 


8 888888 


TABLE V 
(6) 


0°73 
o'74 
o°75 
0:80 
0°83 
°'90 
0°94 
0°98 
0°95 
0°96 
0°96 
0°97 
0°98 
0°99 
1'00 
1°00 
°o°'78 
0°78 
0°78 
o-78 
o°77 
0-80 
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o'90 
0°68 
°'70 
0°76 
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0:98 
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TasB.Le VI 


Annual mean values of total radiation, slowly varying radiation, and quiet radiation 


Unit= 10-2 wm-? (c/s)! 
































Year Total Slowly Varying Quiet 

28 12 6 28 12 6 28 12 6 

1947 220 96 64 94 39 27 126 57 37 
48 176 ene aie 75 ioe nee IOI ee aw 

49 178 76 55 78 26 15 100 50 40 

50 129 61 41 47 16 10 82 45 31 

51 121 40 27 40 12 8 81 28 19 

52 87 34 21°4 15 6 3°6 72 28 18 

53 75 31 20°1 6 2 7 69 29 18 

54 72 29°5 16°7 2 06 o'5 70 29 16 

Tas_e VII 
Radiation associated with specified measures of solar activity 
Radiation | R= A=F=0| R=100 | A=2000| F=2000 

Slowly varying flux 28 ° 55 72 58 
(Unit = 10~-*? wm~? (c/s)~!) 12 ° 20 26 19 
; 6 ° 13 16 13 
Quiet flux 28 70 90 100 gI 
(Unit= 10-2? wm~? (c/s)~") 12 29 4! 46 42 
6 17 27 32 28 
Quiet Sun apparent 28 43 55 61 55 
temperature (Unit= 1000°K) 12 96 137 153 140 
6 230 360 425 370 


























At all three frequencies there is an increase of quiet radiation with increased 
solar activity. At each frequency the flux varies approximately quadratically 
with the solar activity in accordance with the empirical relations (in units 
10-*? wm? (c/s)"1): 


Quiet flux 2870 + 20(R/100)?=70 + 7°5(A/1000)*, 
1229 + 12(R/100)?=29 + 4°3(A/1000)?, 
6=17 + 10(R/100)?=17 + 3°8(A/1000)?. 


The accuracy is not sufficient to establish the quadratic form of this variation 
but the change between sunspot maximum and minimum is much greater than 
any likely error. The relative change of flux is systematically greater for smaller 
frequencies. At sp-min the flux is approximately proportional to the frequency 
(as already found for the slowly varying component) but this does not apply at 
sp-max. 

In Fig. 6 the quiet flux values (expressed as apparent temperature of the 
Sun’s disk) for A=o and A=2000 are compared with a recent survey (18). 
A= 2000 can be regarded as representing the recent sp-max activity. The new 
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Fic. 6.—New values of apparent solar temperature compared with a recent survey (17). 


measurements call for no great change in the previously published curves but 
suggest that the two curves sp-max and sp-min instead of being parallel should 
be closer at the high frequency end and wider at low frequencies. 


6. Conclusion —The main conclusion coming from the analysis is that the 
variation of quiet Sun decimetre radiation is small but real. This conclusion 
differs from others recently published (4, 5, 15). 

It is important to enquire whether the analysis has omitted any factor that 
could influence this result and we are reminded that the complication (/) of 
Section 1 is not yet resolved. Are there any local sources uncorrelated with our 
R, A or F? We consider first whether uncorrelated courses could be expected 
physically, and secondly whether there is any evidence of them. 

The local sources appear to be hot regions, i.e. regions where the high coronal 
temperatures have penetrated into the higher density material (19, 20). There 
is every reason to think that faculae represent these hot regions in the high 
photosphere. We would expect the emission sources to be closely connected 
with faculae and therefore with the other phenomena that show an almost one 
to one correspondence with faculae such as calcium plages, regions of magnetic 
field, regions of 5303 emission, and sunspot groups. The only difference between 
these from the statistical or correlation point of view concerns the delay and mean 
life for which adequate allowance has been made. If the emission is a thermal 
effect we would not expect it from sources which did not show the similar thermal 
evidence. ‘Therefore we would not expect any local sources that are uncorrelated 
with sunspots or faculae. If the same analysis could be carried out using plage 
brightness or area of high 5303 coronal emission it should give the same result 
again. 

Turning to the evidence from observation we find that nearly all the observed 
local sources of emission can be identified (21, 22, 23, 24, 25) with sunspot 
groups, plages, or 5303 regions. However one of the local emissions at the 
eclipse of 1948 November 1 (21) was identified with a stable prominence. Stable 
prominences (seen as filaments on the disk) are practically uncorrelated with the 
solar activity data we have used. Numerical data on prominences from Catania, 
Tokyo, Madrid, and Kodaikanal in the period 1947-54 have been examined 





188 The variation of decimetre-wave radiation with solar activity Vol. 117 


with the view to seeking a prominence component in the slowly varying radiation. 
Unfortunately the prominence measurements do not appear to be consistent 
enough to use them for this purpose. There was only one clear feature that 
appeared on all the data and was not seriously confused with other activity. 
This was a maximum of prominences in 1951 September—October. From Fig. 3 
we see no corresponding increase in 28 or 12 but some evidence of an effect in 6 
(the frequency of the 1948 eclipse result). We might conclude that if there is a 
prominence effect it increases with decreasing frequency. 

These arguments are somewhat contradictory. It cannot be decided with 
certainty whether uncorrelated local sources (i.e. sources that are not associated 
with centres of sunspot activity) exist or whether they produce sufficient radiation 
to affect the issue. The author’s personal opinion is that the total emission from 


such sources is negligible and that the data of Tables VI and VII can be used 
directly for quiet Sun problems. 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 
1957 January. 


References 


(1) A. E. Covington and W. J. Medd, 7.R.A.S. Canada, 48, 136, 1954. 
(2) J. L. Pawsey and D. E. Yabsley, Aust. 7. Sci. Res. A, 2, 198, 1949. 
(3) W. N. Christiansen and J. V. Hindman, Nature, 167, 635, 1951 
(4) J. H. Piddington and R. D. Davies, M.N., 113, 582, 1953. 
(5) M. Waldmeier, Z. f. Ap., 36, 181, 1955. 
(6) A. E. Covington, Tables of daily intensity and outstanding disturbances of 10°7 cm solar 
radiation, 1946-49, N.R.C. Canada, Ottawa, 1950. 
(7) Quarterly Bulletin of Solar Activity. 
(8) F. J. Lehany and D. E. Yabsley, Aust. 7. Sci. Res. A, 2, 48, 1949. 
(9) F. J. Lehany and D. E. Yabsley, Aust. 7. Sci. Res. A, 3, 350, 1950 and private 
communication with J. L. Pawsey and S. F. Smerd. 
(10) M. Waldmeier, Ergebnisse und Probleme der Sonnen-forschung, 2nd edn., p. 145, Leipzig 
1955. 
(11) U.S. Naval Obs. Circulars. 
(12) Greenwich Photoheliographic Results. 
(13) Royal Greenwich Observatory. Summaries of spots and faculae in M.N. 
(14) Royal Greenwich Observatory, Sunpots and Geomagnetic Storm Data, 1955. 
(15) H. W. Dodson, McMath—Hulbert Obs. Repr., No. 41, 1955. 
(16) C. W. Allen, Terr. Mag., 51, 1, 1946 ; 53, 433, 1948. 
(17) G. Righini and G. Godoli, Oss. e Mem. Arcetri, No. 67, 1, 1952. 
(18) C. W. Allen, J.A.U. Symposium No. 4 on Radio Astronomy, p. 253, 1956. 
(19) J. H. Piddington and H. C. Minnett, Aust. ¥. Sci. Res. A, 4, 131, 1951. 
(20) M. Waldmeier (and H. Miller), Z. f. Ap., 27, 58, 1950 ; 40, 221, 1956. 
(21) W. N. Christiansen, D. E. Yabsley and B. Y. Mills, Aust. 7. Sci. Res. A, 2, 506, 1949 
(22) T. Hatanaka et al., Rep. Ionosph. Res. Jap., 9, 195, 1955. 
(23) M. Laffineur et al., Ann. d’ap., 17, 358, 1954. 
(24) A. E. Covington and N. W. Broten, Ap. 7., 119, 569, 1954. 
(25) H. W. Dodson, Ap. 7., 119, 564, 1954. 











sis st ic Si dial en ORS SR 


2 Raa se eal 











se iy it aN 


te RE ANE i Patt Bele Be 


i 
3 
: 
i 
£ 
; 











ea lat csc uly 


eelintinie a wise 
SD a ae nial te id 


_ RE Teint alba coat ictal ea 


Neb RY 5 ate 


Raden IRR 


Ni i as Sat 


AER tnt i a bay 


C1 NL ADEN ok oy 


eS onkepeents aay, 


THE 2°S—2°P, 3P, 35D, 43D EXCITATIONS OF HELIUM 
ATOMS BY ELECTRONS 


B. L. Moiseiwitsch 
(Communicated by D. R. Bates) 


(Received 1957 January 25) 


Summary 


Born’s approximation is used to calculate the cross sections for the 
excitation of the 2°P, 3°P, 3°D and 4"D states of helium from the metastable 
2°S state by electrons with energies ranging from the thresholds of the 
processes to about 100 eV. In the case of the 2°S—2*P excitation a close 
coupling formula due to Seaton is also used. 





1. Introduction.—Goldberg (1) has investigated the relative intensities of the 
helium lines for a nebula in radiative equilibrium, making no allowance for the 
effects of collisions with electrons. ‘The lack of accord between theory and 
observation which was obtained might well be due to the neglect of collisions (2). 
Unfortunately it has not yet been possible to take them into account, due to the 
lack of knowledge of the relevant cross sections. It seems desirable, therefore, 


to carry out calculations on the collision processes which are of importance in 


gaseous nebulae. ‘Those involving the 2°S metastable state of helium are of 


particular interest in this connection. 
In this paper, Born’s approximation has been used to obtain the cross sections 


for the 2°S—2°P, 3°P, 3°D, 4°D excitations of helium atoms by electrons. 
The energy range investigated in these calculations is from the thresholds of the 
processes to about 100eV. Born’s approximation tends to overestimate the 
cross sections, especially when the coupling between the levels involved is strong. 
As this occurs for the 2*S —2P excitation, a close coupling formula based on 
Bethe’s approximation which is due to Seaton (3) was also used in this particular 
case. 

2. Theory.—The Born approximation to the cross section for the excitation 
of the 1sm/*L state of helium from the metastable 1s2s*5S state is given by the 


expression 
i 8am? pKmax, . ac ST \iae ek 
O(2s*S —nl®L)= Fa5 ,. ~ | A(2s8S —nl8L)PKdK (1) 
where 
A (28*S —nlL)= | | {$(2s9S5 ry. e2)O(l Ls ty.) V(t P25 Fs) 
expiK.r, dr, dr, dr, (2) 
and 
K =k; —ky, kj = 27my;/h, kp=27mv_/h; (3) 


vy; and vy are the initial and final velocities of the incident electron; r, and r, are 
the position vectors of the bound electrons relative to the helium nucleus and r; 
that of the the free electron; the ¢’s are the wave functions of the helium atom 





190 B. L. Moiseiwrtsch Vol. 117 


in the states indicated and V is the interaction potential which is given by 





rf 
V(ryt23 ts)= — + = == (4) 
"13 23 3 
If we write the helium wave functions in the form 
d(nl®L; ry, 62) = V2 = {u( 1s; r,)u(nl; r2)—u(is; rp)u(nl; r)} (5) 
where the atomic orbitals u are normalized — mutually orthogonal, we obtain 
N (2838 —nlL) = = " |F (28 —nl)| (6) 
where 
JI (2s —nl)= [ u(as; r)u(nl; r)expit.rdr (7) 
and t=Ka,, (8) 


r being in units of dp, the radius of the first Bohr orbit of hydrogen. 
For the 28S —2°P, 3D and 43D excitations, the helium 2s atomic orbital was 
chosen to have the form 





ans op _ a +B) {1— f(a + B)rje*” (9) 
with “= 2-00, B=0-695 


while u(zp; r), u(3d; r) and u(4d; r) were assumed to be hydrogenic. ‘The 
2*S —3°%P cross section was found to be very sensitive to the parameters of both 
the initial and final orbitals. Consequently, for this particular case, use was 
made of the more accurate atomic orbitals obtained by Morse, Young and 
Haurwitz (4) and by Goldberg and Clogston (5); they are 


u(2s; r)= nigel M(re~“’ — Ae~*) (10) 
V 47 
u(3p;r)= . | 3 cos@Nr(re-” — Be-"*) (11) 
47 
‘with 
M=0°356, A=2:93, a@=1'57, p=0-610 
and 


N=o0-00990, B=14:2, b=0-490, A=0°318. 
In the case of the 2°S'— 2°P excitation of helium, the value of the line strength 
S given by the expression 

S=3| { u(2s; r)u(2p; r)rcos 6 dr)|? (12) 
is 19°3. This implies that the coupling between the initial and final states is 
strong and hence that Born’s approximation must greatly overestimate near the 
threshold. Seaton (3) has shown that in such cases a more accurate procedure 

can be based on Bethe’s approximation 

mi 87a,’ ki tke 

Overt ee “ake Slo 0$ «=F ki — kt’ (13) 


If the coupling is strong the first few partial cross sections Oj, will exceed 
the theoretical maximum 


Oster = Sy (2l+1) 





Pete eee 
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imposed by charge conservation. Suppose /, is the value of / for which 
Osetne Oe. (Sh) 
<@i. (i> ly). 
Putting O'=30,!,, for = and O'=OQ,,,. for >, Seaton obtains 


7 , lo 
Q _ ki? (2, + 1) + Osetie re 2 Bethe (14) 


This formula, as well as Born’s approximation, has been used in the present 
paper to obtain the cross section for the 2S —2°P excitation. 

3. Results.—Displayed in the table are the values of the cross sections for the 
2°S—25P, 3°P, 3°D and 4°D excitations of helium «atoms by electrons with 
energies ranging from the thresholds of the various processes (1:14eV, 
319 eV, 3:25 eV and 3-92 eV respectively) to about 1ooeV. 

TaBLe | 
Cross sections Q for the 2*S—2°P, 3°P, 3°D and 4°D excitations of helium atoms by electrons 


QO (2°S—2°P) 





Electron (in 77a,") QO (2°S— 3°P) O (2°S—3*D) QO (2°S—4°D) 
Energy , A . (in za,") (in 7a") (in 7a”) 
(in eV.) Born Equation (14) 
1‘14 ° ° 
222 173 
1°72 79 
2°18 294 94 
3°40 248 110 1:08 3°60 
4°89 202 114 1°29 6°94 1°69 
6-66 167 112 1*13 6°34 1°73 
8-70 138 105 1°07 5°41 1°51 
13°6 IOI gI 0°99 3°85 1°09 
19°6 76 "7 ool 2°81 080 
30°6 54 60 0°78 1°87 0°53 
54°4 34 41 0°59 1°08 O31 
85°0 24 30 0°46 0°70 0°20 
122 18 23 0°37 0°49 O14 


The thresholds for the 2°S—3°P, 3°D and 4°D excitations occur at 3:19 eV, 3°25 eV, 
3°92 eV, respectively. 


The effect of using Seaton’s modification of Bethe’s approximation instead 
of the Born approximation is to reduce the cross section for the 2*S—2°P 
excitation by nearly one-third close to the threshold and also drastically to alter 
the shape of the peak. At higher electron energies the two approximations 
agree quite well. 

The cross section for the 2?S—3%P excitation is very much smaller than 
that for the 2°S—23P excitation, their maximum values being about 1-37a,? 
and 1157a 9" respectively. The former cross section is abnormally small, mainly 
due to the fact that the line strength S for the 2°S—3°P transition is only 0-92. 
Consequently, the 23S —3*D and even the 2°S— 4D excitation cross sections 
are larger than the 2°S — 33P cross section near the threshold, having maximum 
values of about 77a)? and 1-87a,? respectively. 
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RADIAL VELOCITIES AND SPECTRAL TYPES 
IN THE GALACTIC CLUSTERS M25 AND NGC 6087 


M. W. Feast 
(Communicated by the Radcliffe Observer) 


(Received 1957 January 21) 


Summary 


67 spectra of 42 stars in M 25 and 38 spectra of 25 stars in NGC 6087 have 
been used for the determination of spectral types and radial velocities. The 
majority of the stars observed in each cluster are late Band early A main sequence 
stars. The radial velocity results support the contention that the classical 
cepheids U Sgr and S Nor are members of M 25 and NGC 6087 respectively. 





Introduction.—The southern galactic clusters M 25 (1.C. 4725) and NGC 6087 
are of interest since near each of their centres is a classical cepheid (U Sgr and 
S Nor, respectively). 


Cluster a(ig50) 8(1950) l b Type (1) Distance (1) 
M 25 18529 —19°10 3413 —6:0 2b 980 psc 
NGC 6087 1615 —57°47 295°3 —6°3 1-2b 870 psc 


It is clear that if the cepheids are actually cluster members we have an 
independent method of determining the zero point of the period—luminosity 
relation and the clusters may also assume some importance in theories of the 
evolution of galactic clusters. Dr J. B. Irwin and Dr A. J. Wesselink both drew 
the writer’s attention to the interesting problem of these clusters and, during his 
. recent visits to the Royal Observatory, Cape and the Radcliffe Observatory, 
Pretoria, Dr Irwin obtained three colour photoelectric observations of the 
cepheids and of a number of stars in the two clusters. ‘These observations are 
now being prepared for publication. It was clear that spectroscopic observation 
of the stars in the two clusters was desirable and it was arranged with Dr Irwin 
and Dr Wesselink that the writer should carry this out. -The purpose of the 
spectroscopic programme was twofold : 

(a) To compare the radial velocities of a number of stars in each cluster 
with the radial velocities of the two cepheids; this is important since, according 
to Irwin (2), the available proper motion data are too poor to decide whether or 
not the cepheids belong to the clusters. 

(5) ‘To classify the cluster stars on the MK system for combination with the 
photometric results in a discussion of the clusters. 

Observations.—The observations were made between 1955 July and October 
using the two-prism spectrograph at the Cassegrain focus of the 74-inch Radcliffe 
reflector. Either the f/3-7 camera (49 A/mm at Hy) or the f/2 camera (86 A/mm 
at Hy) were used. In the case of the f/2 camera the spectra were considerably 
widened for spectral classification. The present discussion is based on 67 spectra 
of 42 stars in M25 and 38 spectra of 25 stars in NGC 6087. 
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The majority of the stars observed in the two clusters are main sequence 
stars of early A or late B type. For these stars the wave-length system recom- 
mended by the Victoria workers has been employed with both cameras. Tests 
of this system on the f/3°7 camera have been reported previously (3). For the 
f/2 camera it is intended to publish results of tests on standard stars later. So 
far such tests have shown no significant systematic errors between the two 
cameras. The lines in the AV-BV stars are generally broad and few in number 
so that it would not be possible to determine accurate velocities for these stars 
without a much more extensive programme than that reported here. It has 
therefore been decided to give only the mean velocity of the main sequence 
stars for each cluster. The programme also included a few late-type stars. 
Spectra of these with the f/3-7 camera were measured, as is the normal practice 
here, using the Victoria wave-lengths (4). ‘Those on the f/2 camera were 
measured on a Hartmann spectrocomparator, either y Hya (gG 5) or HD 119971 
(K 5) being used as standards. 

The spectral types were determined by comparison with standard stars on 
the MK system (§) and by use of the Yerkes Atlas. 

Discussion.—The radial velocity results and their standard errors are given 
for the two clusters in Table I. 


TABLE | 
Radial velocities in M 25 

U Ser 4 km/sec 

54 spectra of 35 stars (BV- AY) 4+4 km/sec 

No. 150 (Gs III) 2+5 km/sec 

No. 193 (G2 III) 1+2 km/sec 

No. 26 (K2 V:) —35+3 km/sec (non-member) 
No. 70 (Ag V) —56+8 km/sec (non-member) 

Radial velocities in NGC 6087 
S Nor 


3 km/sec 


27 spectra of 18 stars (BV—AV) 2+ 3 km/sec 


No. x (Mrz III) 2+4 km/sec 
No. f (G4 IID) —27+2 km/sec (non-member) 


No. B (A4 V) —14+3 km/sec (non-member) 
No. t (Bog III:) —10+7 km/sec 


The letters and numbers used to designate the stars in Tables I and II are 
due to Irwin and will be identified by chart in his forthcoming discussion of the 
photometry of the clusters. 

The normal radial velocities of U Sgr and S Nor are taken from the recent 
investigation by Stibbs (6). According to Stibbs “ the estimated average probable 
error of a single normal velocity is of the order of 1 km/sec”’.. The 35 stars in 
Mzs and 18 in NGC6087 used in forming the mean “main sequence” 
velocities in Table I are indicated by asterisk in Table II. ‘There would appear 
to be little doubt that they are all cluster members. ‘The agreement of the mean 
““main sequence” velocities in each case with the velocities of the cepheids 
strongly supports the contention that the cepheids are indeed cluster members. 

Stars other than the BV—AV stars are listed separately in ‘able I. However, 

omitted from this table are the M giants 6087/e and 25/49 for which only one 
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f/2 spectrum each was available. The velocities obtained for these two stars 
are not of sufficient accuracy to decide on cluster membership. The A type 
supergiants 25/124 and 25/192 are also omitted from Table I. They are 
presumably non-members though their radial velocities are too poorly determined 
to draw definite conclusions from,them. The stars 6087/g (B8V or III) and 
6087/v (AoIII: (V?)) are the only other stars not listed in Table I. It seems 
likely that they are main sequence cluster members and that their velocities 
might be incorporated in the ‘‘main sequence” mean for NGC 6087 but in 
view of the uncertainty in their spectral types it has been thought best to omit 
them. 

Turning now to the stars listed separately in Table I, the two clusters may 
be discussed in turn. In M25 the G-type giants No. 150 and No. 193 appear 
from their radial velocities to be cluster members, a conclusion that is not at 
variance with a rough estimate of their magnitudes compared with those of the 
main sequence cluster members. On the other hand the stars No. 26 (K2 V:) 
and No. 70 (Ag V) whose velocities indicate that they are non-members appear 
in fact to be too bright for membership and to be foreground objects. In 
NGC 6087 star No. x (Mr III) appears to be a member from its radial velocity 
whilst No. f (G4 III) seems to be a non-member. This latter result is confirmed 
by three f/2 and one f/3-7 spectra. 6087/8 (A4 V) was not included in the main 
sequence mean as it is clearly a non-member. It lies on the outer edge of 
NGC 6087 and is the brightest star in the region of the cluster (brighter than 
SNor itself). ‘Thus since it is a main sequence object it must lie in the 
foreground of the cluster and it is reassuring that the radial velocity result confirms 
this. In the case of 6087/t the radial velocity results do not give any very clear 
indication of membership. The star is classified as Bg III: it is peculiar in 
that it possesses hydrogen lines having very sharp absorption cores and 
moderately extensive wings. 

The standard errors given for the BV-AV means in the two clusters 
(+4 km/sec for M25 and +3 km/sec for NGC 6087) include both observational 
errors and true velocity dispersion in the clusters. However since more than 
one spectrum of a good many of the stars is available it has been possible to 
estimate the standard errors of the means to be expected from observational 
errors alone. In both clusters this is found to be +3 km/sec so that the velocity 
dispersions are too small to be detected from the present results. In view of the 
broadness of the lines of the majority of main sequence cluster members, a 
very extensive programme would be necessary to obtain radial velocities of 
sufficient accuracy for studies of velocity dispersion 

Table II lists the spectral types in the two clusters. It also indicates the 
number of plates available with the two cameras. HD types are shown where 
these are available. In two cases, 6087/h and 6087/0, Hf emission is suspected. 

It is not profitable to enter into a discussion of the spectral types at this stage 
since magnitudes and colours of the stars will soon be published. However the 
spectral types alone do not indicate any abnormalities in the clusters. It is 
perhaps desirable to point out that, whilst the majority of stars are classified as 
late B-early A luminosity class V and have been referred to as main sequence 
stars, it is not beyond possibility that some of them may be luminosity class IV. 
Distinction between luminosity classes V and IV appears to be somewhat 
difficult in the case of broad-lined stars of early type. 
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C.P.D. No. 
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TABLE II (continued) 
Spectral types in NGC 6087 


Irwin Cordoba Spectra 
No. No. fiz f13°7 HD .MK Type 
a 19 I oa : A1 V 
b 23 I * Bo B8 V 
c 27 I ° te Bo V 
d 30 I - 2 B8 V 
e 28 I Ma M3III: member? 
f 26 3 I Ko G4lll non-member 
g 33 I Ao B8 V or III member? 
h 39 I ete si B8 V(e) 
i 40 I I *  B8 Bs V 
j 47 I I * Bs B6 V 
k 61 I . B8 V 
] 45 I - 2 Bs V 
m 55 I I ss Bs V 
n 51 3 be Ao V 
o 53 I wa ee Bs V(e) 
p 76 I I * Bs B7 V 
q 78 1 ig A2 V 
r 79 I * B8 V 
s 80 I I > 8B B6 V 
t 84 I I B8 Bg III: member? 
u 83 I * Ro Bo V 
v 89 I oa Bo Ao III : (V?) member? 
x 73 I I Ma MrlIill 
y 87 I . Ao V 
B 9 I Az A4V non-member 


Notes to Table IT: 
(1) The stars used'in forming the BV-AV means of Table I are denoted by an asterisk. 
(2) The stars in M 25 are identified by their numbers in the Cape Photographic Durch- 
musterung (Cape Annals Vol. 3). The stars in NGC 6087 are identified by their 
numbers in the Cordoba Photographs of Clusters (Cordoba Resultados Vol. 19). 
NGC 6087 = Cordoba A326 (Norma). 
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MOTIONS OF STARS OF TYPE A 
PERPENDICULAR TO THE GALACTIC PLANE 


R. v. d. R. Woolley 
(Received 1957 January 2) 


Summary 


Formal analogies are examined between density in directions perpen- 
dicular to the galactic plane and density in globular clusters. In the light 
of these, and otherwise, data on the density of A-stars as a function of z 
(including recent results published by van Rhijn) and the statistics of A-star 
radial velocities (from Wilson’s recent catalogue) are examined, with the help 
of Jeans’s theorem. It is concluded that the data demand a density of 0-2 solar 
masses per cubic parsec in the neighbourhood of the Sun. 





1. The relations between density, potential and stellar velocities of the 
material in the neighbourhood of the Sun offer to a limited extent some analogies 
with the corresponding quantities in a globular cluster. ‘Two methods of studying 
the former are open to us: we may either assume some morphology or model, 
and investigate the dynamical consequences of what we have assumed, or we may 
investigate the observed facts about stellar densities and velocities and find out 
as much as we can with the minimum appeal to theory. Since neither approach 
can as yet be entirely satisfactory, some sort of compromise may be attempted. 

Considering briefly the first or model approach, we may make as a first 
approximation the very simple assumption that the material is stratified in 
parallel planes—planes, that is to say, parallel to the plane of symmetry of the 
Galaxy. Then the densities and motions perpendicular to this plane exhibit 
a marked resemblance to the densities and motions in a globular cluster. In 
fact if all the mass is collected into stars, or if in some way the residual matter 
(dust and gas) distributes itself as do the stars of mean density, we have simply 
to deal with a one-dimensional case of the globular cluster. 


This may be quickly explored. Jeans’s relation* is deduced from Liouville’s 
theorem as follows. If 


f(x,y, 3, u, v, w) dx dy dz du dv dw 


is the number of stars in the volume dx dy dz having velocities between u and 
u + du etc., then in a steady state, i.e. provided that f is not a function of the time, 


Df|Dt =o. 


Jeans shows that if »(V,¢) dV is the number of stars per unit volume with speeds 
between V and V+dV at a place where the gravitational potential is ¢, the 


* Jeans, J. H., M..N., 76, 78, 1915. 
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vectors being distributed isotropically, and if 


WV, >) = V*f(V, 9), 
f(V,6)= F(V? — 29) (1.1) 


this being Jeans’s relation.* Its equivalent may be deduced by enumerating 
stars in various 2 levels in the case where ¢ is a function of a single coordinate z 
alone (Oort) or in various concentric shells where ¢ is spherically symmetrical 
(Woolley), but the particular investigations make no more and no less appeal 
to “equilibrium” than does the deduction from Liouville’s theorem: in all 
cases we simply assume that the state of the stellar motions is secularly stable. 
In the three-dimensional case we must also assume that the distribution of the 
vectors is isotropic, but if the analysis is restricted to a single coordinate this 
requirement vanishes.t However, if the velocities are Z restricted to the z 
coordinate, to obtain Jeans’s relation we must set 


(Z, $) = f(Z, 9). 
A particular solution of Jeans’s relation is 
F(Z? — 2) =const. x exp [ —j7(Z? — 2¢4)], 
where j is independent of Z and ¢ but may depend on the mass or class of the star. 
This is a requirement of full equilibrium, i.e. equilibrium stable against the 
effect of stellar collisions; but Oort employs a Gaussian distribution of stellar 
velocities (or a sum of such distributions) without supposing it to have been 
set up by stellar collisions. Instead, he infers it from radial velocity statistics. 

Jeans’s relation informs us that —2¢ must enter the distribution in the same 
way as Z*, If there are m classes of stars, of which the pth class has a Gaussian 
velocity distribution with a parameter j,,, then 

F(Z? — 2¢) = A, exp [—j,(Z? — 2¢)]. 
More generally, from a complete knowledge of the stellar velocity distribution 
we can construct the potential ¢, or its derivative the acceleration K(z). 

To proceed to the analogies with globular clusters, we examine theoretically 
the supposition that ¢ is a function of z alone, in three cases. 

Case I. Consider first the simple case where j has the same value for all stars. 
Let =o be the plane of symmetry of the Galaxy and let ¢, be the value of ¢ at 
z=0. Write 2j°7(4—¢))= —y%. Then 

u(Z, ¢) =const. x e*e-¥ 
and by integration with respect to Z, v(¢)=const. xe’. If p is the density and 
P/ Py = then 


then 


“eé 


n=e". 
We suppose that if there is any appreciable mass not concentrated into stars, 


* So called by A. S. Eddington (M.N., 76, 572, 1916) a few months after Jeans published the 
original paper. But most writers quote the relation in a different form, for example W. M. Smart, 
Stellar Dynamics, pp. 324-326, 1938. 

+ The empirical facts concerning stars in the neighbourhood of the Sun are of course that the 
velocities are not distributed isotropically but ellipsoidally, and the observed ellipsoid does not 
coincide with the ellipsoid deduced from Jeans’s relation for a rotating galaxy. (W. M. Smart, 
Stellar Dynamics, p. 395, equation (13); v. d. Pahlen, Einfiihrung in die Dynamik von Sternsystemen, p. 
126, footnote (18); S. Chandrasekhar, Principles of Stellar Dynamics, p. 125.) It is however quite 
natural to suppose that the state of the velocities in the galactic plane is not steady, yet the state of the 
velocities normal to this plane is steady, in the sense that the average numbers of stars crossing each 
plane parallel to the galactic plane are equal in the two senses upwards and downwards. 
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its density, as a function of 2, is similar to that of the stars. Poisson’s equation 
V°¢= —47Ip becomes 
d*4/dz* = — 471 poe” 
and if /=(87I'p,j?)-!? and x= 2/1 we have 
d*b/dx* =e” (1.2) 
subject to the boundary conditions ys;=o and =o at x=0. The solution is 


b= 2 Incosh (x/ 4/2) | (1.3) 


n=e-*=sech? (x/1/2) 
Equation (1.3) corresponds to the “‘equation of the isothermal gas sphere” in 
the case of spherical symmetry* but the solution does not give an infinite mass 
in a cylinder of finite cross-section: we have in fact 
@ 
| pdz=2y/2lpo, 


so that there is no need to invoke a “‘truncation”’ of the velocity distribution to 
get a finite mass. 

From equation (1.3) we can calculate the periodic time of a star to move up and 
down in z. Let the motion be such that when furthest from z=o the potential 
is %. At this point the velocity is zero and elsewhere 

jZ=(fo—4)*. 
The periodic time T is therefore 
Vo 
T=2y2.jl { (bo — )-2%(1 + e-¥)-12 di, 

0 
and since 2>1+e~-*>1 for all positive x, 

T = 4jliby'? x c (1.4) 
where 4/2>c>1, a formula which gives the order of magnitude of 7 without 
having to work out the integral. 

Case II. We next consider the case where the stars have a distribution in 
mass and where j? is proportional to this mass. This is the full equilibrium 
case. Select a suitable mean mass m and let p=m/m. Let v, be the number 


of stars per unit volume at z=o and let the number of stars per unit volume 
with reduced masses between p and n.+dy be 


veo) dy. 
In particular we consider the case & (u)=const. x p“e~’". For each mass the 
velocity distribution is proportional to exp (—,j?Z*) ex hypothesi. Set 


2j°(—$) = —. 
Then Jeans’s relation gives 
v,(Z, ) = const. x e""2" eo" 
and, by integration with respect to Z, 


v,(ib) = voEo(ue”. 


Since the density p is proportional to m pE (je du we have 
, 0 


nD 
p=const. x (b + s)-@+*) sitie-= dx 


“0 





* ef. G. 1. Camm, M.N., 110, 309, 1949. 
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and 
1=plpo=(1 + p/b)@*. (1.5) 
Using, as before, z=/x where now / = (8xI"p)j-2)-, Poisson’s equation gives 
ds b \at2 
a ~ (535) tin 


subject to {=o and b=0 at x=0. 

The value a=1 fits the distribution of stellar mass in the neighbourhood of 
the Sun reasonably well, and we may set b= 3 which defines the mass m. (Choice 
of m is a choice of /.)* ‘Then (1.6) becomes 


d*p/dx? =(1 + hy) 


of which the appropriate solution is 


1+4p=(1 +422)! | 


1=p/po=(1 +4x")*? 


The periodic time to move up and down is roughly 47j/ for moderate values 
of pp. 

Case III. Oort’s empirical representation of the velocity statistics may be 
examined. Suppose there are n classes of stars and that for the pth class the 
z component of velocity is distributed as exp(—j,?Z*). Suppose that the mass 
of each star in the pth class is m, and that the abundance of the pth class at z =o 


is &,(p) subject to 


(1.7) 


n 
> &o(P)=1. 
p= 
By Jeans’s theorem 
v,(Z, 6) =const. x exp [ —j,,°(Z? — 2¢)} 
and 


"=pife= 2 m,€o( Pp) exp [27,,°(¢ — dy)] > m,€o( p) 
po p= 


- 2 A, exp [2j,°(¢ — $o)] 


subject to } A,=1. Select a suitable value of j,”, namely j”, and set 
pi 


a,=j,2)%  %=2Xbo—d), [= (8aT pyj®) 22, x =s/l 
Then Poisson’s equation becomes 


d*b/dx? = y A,, exp (—4,) (1.8) 
p-l 

and we may impose the condition } A,a,=1, which merely defines p. In 
pl 


setting up this equation we suppose that ¢ is a function of z alone—that is, that 
the distribution of matter is uniform over planes normal to the z axis over distances 
large in comparison with the values of z with which we propose to deal. We 


* In setting up a model. In a real analysis of stars we know too little about the mass function 


to be sure of m, 
* 
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also assume, as usual, that the dust and gas are distributed as are the stars of a 
given class (or a set of classes). The solution in series of (1.8) is 


p= px2— tet +.. 
y=t— batt SA, = oe (2.9) 


The three solutions are very similar near x=0; we have in fact, from (1.3) 

p= hx?—-txt+..., n=1— }x?+ixt-..., 
from (1.7) : 

p= tx?—txt+..., n=1—4}x2+224- 
Numerical values of y% and d/dx from equations (1.3) and (1.7) are given in 
Table I. It will be seen that they are very similar as far as the table goes: either 
set represents the potential and its derivative as well as can be required, in view 
of the doubtfulness, in any application to the Galaxy, of the legitimacy of 
leaving the dust out of account and of assuming uniformity in planes normal 
to the z axis. In Table 1, z=distance from galactic plane =/x ; ¢= gravitational 
potential, d= —2j*d$; K(z)=acceleration= —dd/dz =(2)2/)1 dib/dx. In this 
1=(87I"p, py! 2 where py is the density at =o and j i is the mean parameter in the 
velocity distribution. 


TaBLe I 


Potential and acceleration for two solutions in which matter is distributed uniformly in 
infinite planes 


Reduced distance Case I Case II 
from galatic plane of All velocity distributions Equation (1°7) 
symmetry alike. Equation (1-3) 
x wb dis/dx yp dis/dx 
° 0°000 0-000 0*000 0°000 
o'2 0°020 0'200 0°020 0'200 
o"4 0079 0390 0079 0°390 
06 O°175 0°566 O'175 0°567 
o'8 0°304 0°724 0°304 0°726 
1‘O 0°463 o0'861 0464 0866 
1°5 0'962 I‘112 0'969 1°134 
2°0 1°557 1°256 1582 1°310 
2°5 2°207 1°334 2268 1°424 
3°0 2°885 1°374 3°000 I*500 
4°0 4°278 1°404 4°550 1*589 
oe) oa) 1°414 oa) 1°732 


The two solutions shown in Table I start together and diverge very little 
from each other even as far out as x=2. The acceleration dis/dx derived from 
(1.7) is pe against Oort’s values* of Ks) in Fig. 1. Since Oort adopts 
log p= —az*, where —a=g x 10-6 when 2 is in parsecs, he implies a value of /. 
We have, from (1.7), for small x, Inp= —}x?. This is equivalent to Oort’s value 
if /=155 parsecs, and this is adopted in drawing Fig. 1 . A further constant 
must be chosen and in fact dys/dx=1 was taken to give K’ (z)=3 x 10-* cm/sec’. 
Then if jis in km/secand /is in cm, (2j2/)-! = 3 x 1 10~® giving j= = 1-71 x 10%cm/sec 
or 17°1 km/sec. The corresponding value of py is 0-11 solar masses per cubic 
parsec, which is close to the value 0-09 given by Oort. 


* Oort, J. H., B.A.N., 6, 262, 1932'(Table 14). 
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2. If however one makes use of the published data in a way differing from 
Oort’s, one arrives at substantially different values of py. This is because the 
density that one finds is proportional to the square of a mean velocity that has 
to be adopted at some point. ‘ 
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Oort uses star counts by van Rhijn to find log, as a function of z for various 
absolute magnitude classes, for each of which he constructs a value of the mean 
square velocity as a function of z. ‘This is a complicated matter, and one may 
avoid complication by considering instead of all stars (or all stars of a given 
absolute magnitude), all stars of a particular class that can be isolated in some 
way. At present the best data seem to be those given by A stars. For these the 
density as a function of height z has been given by Lindblad.* It may also be 
estimated from counts in the Henry Draper Catalogue (a device used by Oort) 
and it is the subject of a recent study by van Rhijn.t As the results agree quite 
well, one may have some confidence in the observed distribution. The results 
are shown in Fig. 2. The values derived from the H.D. counts require explanation. 
Oort selects the area from 11" to 14"30™ in R.A. and from + 10° to + 50° Dec., 
with an extension to +5° Dec. between 11"30™ and 14", to represent the N 
galactic cap. Counts in this area in H.D. give 


Visual magnitude <5 5 to 599 6to6gg 7 to 799 8 to 89g 
Types Ao and A2 7 26 43 50 85 
Types A3 and As5 3 13 31 38 58 


‘These may be converted roughly to densities by using the simple assumption 
that all Ao and A2 stars have absolute magnitude +0°5 and all A3 and As stars 
have absolute magnitude +1-5. (Inspection of the count shows that this will 
work quite well, with approximately equal densities of stars in the two groups 
at each distance, since the number of A3 and As stars in the (m+ 1)th interval 
is always roughly equal to the number of Ao and Az stars in the mth.) Then the 
Ao and Az stars of apparent magnitude 5"™-o are distant 79:4 parsecs and those 
of 6™-o are distant 126 parsecs, so that the stars between 5™ and 6™ lie in a volume 
proportional to 126° minus 79°4°, etc. ‘The mean density of stars of apparent 
magnitude 5™-5, distant 100 pc from the Sun (and therefore 130 pc from the 
central plane of the Galaxy{) is taken to be proportional to the number counted 
(26 in this case) divided by the relative volume (126*— 79-4) and so for other 
entries. In this way we find, combining the two groups, 


Distance from the central 
plane in parsecs 93 130 190 280 430 

log relative number of stars 1°75 1°67 1°13 2°75 2°40 
The logarithms have an arbitrary added constant since the brightest stars are 
too few in number to give reliable densities for small values of z. 

If the distribution of velocities of the A stars is a simple Gaussian one, say 

n(Z)=const. x e“"4", 
then the distribution of density with potential is, by Jeans’s theorem, 
p=pe"” 
where os = 2j*(¢) — ¢),, and then 
log p — log py = — 0°43.43 pup, 

so that the distribution of density with galactic height z in Fig. 2 gives a direct 


* B. Lindblad, Ark. Mat. Ast. Fys., 19b, No. 15, 1927. 

t P. J. van Rhijn, Publ. Kapteyn Ast. Lab., Groningen, No. 57, 1955. ‘Table 34. 

} Throughout this paper it is assumed that the Sun is + 30 pe from the central plane of the Galaxy, 
as used by J. H. Oort, B.A.N., 236, p. 260, 1932. P. J. van Rhijn finds a lower value. : 
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determination of (a quantity proportional to) the potential % as a function of z 
This may be compared with the potential in ‘Table I, by adopting values of I 
and. In practice the quantity /u~? is fixed by the observations,* and /= 150pe, 
f/=2'5 give a good fit to the observations from z=0 to about s=250 pe, at 
which point the potential deduced from the observed densities begins to diverge 
from the theoretical potential of Table I (equation (1.7)). This is shown in 
Fig. 3. Ignoring this divergence for the moment we have /u~!?=95 parsecs. 
Since the average Z velocity is given by |Z|= (mj )? and since /? = (87I"p)j?)"! 
we have 

8I'py = w(|Z|)?0-* 
so that with /u-!?=95 pe we have the following table. 


Mean velocity Density p, 
|Z| Solar masses per 
km/sec cubic parsec 
6 o'12 
7 0°16 
8 O'2I 
9 0°26 


3. We now turn to the observed radial velocities of the A stars. Oort gives 
some results from data available when he wrote his paper which was published 
in 1932. ‘There is now available an extensive cataloguef of radial velocities 
published by R. E. Wilson in 1953. ‘This gives 124 stars of types Ao to A8 
(including two which are double stars with one F-type component) in Oort’s 
northern galactic cap area and 32 in the southern area. When these are corrected 
for the solar motion the velocities are grouped as shown in Table II. 

Tasce II 


Radial velocities of A stars in Oort’s galactic cap regions corrected for the solar motion 


Velocity No. of stars Gaussian 
km/sec N S Together distribution 
oto 4 40 10 50 54 
5to 9 51 9 60 458 
10 to 14 17 9 26 31 
15 to 19 6 I 7 15 
20 to 25 8 I 9 6 
30 we I I > 
34 I ae I . 
36 I sey I 24 
41 sia I I 
Average velocity | Z|, km/sec 79 9°6 8°3 
Median velocity M km/sec 6°7 7°8 6°8 


The Gaussian distribution shown is 156 stars with a median velocity of 6:8 km/sec, or 
j=0'070. 

In forming ‘Table II no correction was applied to allow for the fact that since 
the area extends to about 30° from the galactic pole, the radial velocities are not 
strictly speaking Z velocities. It is not possible to supply reliable Z components 
due to motion perpendicular to the line of sight, and it was simply supposed that 
the distribution is statistically similar to a Z distribution. 

* because, for relevant values of ys, 4 ~x!/* approximately. 
+R. E. Wilson, General Catalogue of Radial Velocities, Carnegie Institute of Washington 
Publication 601, 1953. 
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A survey was also made of radial velocities of stars brighter than 8™-o, of 


types Ao and Aon only, within 40° of the galactic poles. The results (for 
velocities corrected for the solar motion are as follows. 


Velocity No. of Ao Gaussian 
km/sec stars distribution 
oto 4 44 39 
5to 9 28 35 
10 to 14 21 22 
15 to 20 10 II 
20 to 24 4 4 
25 to 29 5 2 
30 and above 2 j 
Average velocity | Z| 8-9 km/sec. 
Median velocity 4/7 6-8 km/sec. 


The Gaussian distribution has a median velocity of 6-8 km/sec. 


These two surveys are not independent, as some of the data overlap, but the 
agreement is satisfactory. For a Gaussian distribution the average without 
regard to sign is 1:24 times the median, and 6-8 x 1:24= 8-4, in fair agreement 
for both surveys. ‘Taking these results as they stand, the average velocity |Z| 
is between 8 and g km/sec, and from the table given at the end of Section 2, 
py lies between 0-21 and 0-26 solar masses per cubic parsec—at least twice Oort’s 


value. There is always the possibility of a selective effect in any catalogue of 


radial velocities, some faint stars having been put on a programme because of 
large proper motion: especially in the Southern Hemisphere where the catalogue 
is very incomplete. For this reason the median should be preferred to the mean, 
as a measure of the velocity distribution. However, the median is reduced 
slightly, but perhaps significantly, if one rejects velocities whose quality is given 
as c or d by Wilson, when the median weneny, (of 54 Ao and Aon stars within 
40° of the pole) drops to 6-0 km/sec. 

4. In this section we discuss the effect bf a double Gaussian velocity 
distribution, if present, on the density, and examine the evidence for a distribution 
of this kind. Since the work of Section 1 shows that the potential % as a function 
of reduced galactic height x is not sensitive to a choice of model velocity 
distribution we supposed that for x <2 we may assume the values of ¢ in ‘Table II 
(Case II) as valid with sufficient accuracy. If now for a particular class of star 


n(Z) = const. x {(1 —0)j,e-%"2" + Oj,e-F*" (4.1) 
at x=o0. Then by Jeans’s theorem we must have elsewhere 
n(Z) =const. x {(1 —0)j,e°3"e-™ + Oj,e-F° Me" (4.2) 
where j,;7=,j" and j,?=p.j*._ Integrating with respect to the velocity Z we 
have, on division, 
1= plpp=e™*{1 + Ae” — 1)} (4-3) 
the relative density » referring to the particular class of star. 
Oort assumes a distribution of A-type velocities in which 6=0-043 and 
Jolj1= 20/84. For such a case, equation (4.3) gives values of » which are very 
similar to that given by a single Gaussian distribution, namely 
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for small values of 4, but shows important differences when ys becomes larger. 
lhe following short table shows how with a suitable adjustment of the constants 
the partial densities can be made practically identical out to z=200 parsecs. 


Taste II] 


Star densities arising from single and double Gaussian velocity distributions 


Single Gaussian Double Gaussian 

distribution distribution 

x ob z,pe —logy z,pe —logy 

o's 0°304 120 0°33 132 o°31 

1°O 0°464 150 0°50 140 0°47 

1°5 096g 225 1°05 210 o'9! 
2'0 1°582 300 272 280 1°2 

2°5 2°268 375 2°46 350 1°46 

3°0 3°000 450 3°26 420 1°52 


In this we have taken /= 150 pe and x =2°5 for the single distribution, and 
/,=140 parsecs and px, =2°5 for the double distribution. If the velocity distri- 
butions give the same average velocity at s=o0, then 


jrv= 1-138) 


and the central density (of the whole material) deduced from space and velocity 
distributions is reduced by about ten per cent when passing from the single 
distribution to the double distribution (i.e. in the ratio j?/*/j,7/,7). Fig. 3, in 
which both theoretical distributions are plotted against the observed densities 
of A stars, shows that the double velocity distribution gives a much better 
representation than does the single distribution, for z>200 parsecs. 

But there is no real support for the double Gaussian distribution in the 
statistics of A star radial velocities: in fact the contrary, since from equation (4.2) 
we can easily show that the average velocity should increase quite markedly 
with increasing ‘y%. Since this corresponds to increasing z or increasing distance 
from the Sun, the double Gaussian distribution requires an increase of average 
velocity with increasing stellar magnitude, and this is not observed. We have 
in fact 


|Z| = const. x il + o1g3e mv \ {0°957 fi 0-043 MW 


Now assume that all A stars have absolute magnitude +1: then the distance 
of those of apparent magnitude 6™-o is 100 parsecs, or z= 130 pe and we have: 


Apparent magnitude 2, parsecs fy Average velocity 
km/sec 
61-0 130 0°50 8:8 
J™-o 188 1°82 10°2 
Pa 230 2°34 15°2 


But the average velocity of 44 stars with m (vis)<5"'9 is |v] =7-8 km/sec, 
and the average velocity of 25 stars with m (vis)>7"0 is |v|=8-0 km/sec, 
and these figures are very strongly against the increase of velocity with distance 
required by the double Gaussian distribution. 
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5. If we reject the double Gaussian velocity distribution and accept the 
single Gaussian distribution as a satisfactory representation of the velocities, 
then Fig. 2 defines the behaviour of ¢ and it is quite clear that % as a function of = 
has a point of inflection. ‘This is impossible if the material is distributed 
uniformly in parallel planes, whether it consists of stars which can be dealt 
with by the methods of Section 1 or whether it consists partly of dust and gas. 
But of course there is nothing to say that the material of the Galaxy is distributed 
uniformly in parallel planes—quite the reverse: there are spiral arms. ‘This 
suggests that we examine the potential of matter distributed in (infinitely) long 
cylinders. We shall have to abandon Poisson’s equation and consider mere 
models, because the spiral arms are certainly not equilibrium figures unless 
galactic rotation is taken into account, and not even then, considering the 
difficulty of accounting (on equilibrium theory) for the observed axes of the 
Schwarzschild velocity ellipsoid. In other words we abandon the theoretical 
approach of Section 1 and seek a compromise between some vestige of theory 
(Jeans’s theorem) and the empirical facts. This consists of setting up as a model 
‘spiral arm”’ an infinitely long cylinder of square section with a uniform density. 

The potential of an infinitely long cylinder of uniform density p and square 
section of side 2a at a point distant z from the centre of the square, situated along 
a line passing through the centre normal to a face, is given by 

$y) —¢= 4! pa*[F(z/a+1)— F(z/a—1)—2F(1)] for z>a | , 
= 4) pa*[ F(1 + s/a) + F(1 —2/a)—2F(1)] for <a | (5-1) 
where ¢, is the potential at the centre and 
F(x) = 3(x?—1) cot-+ x + 7/4—}x[1 —In(1+.2°)). 
A point of inflection in ¢ occurs where s=a, and by inspection of Fig. 2 if this 
model is to represent the observed potential we must set a=150 parsecs (or 
possibly a little larger). Since F(2)=2-0g0, F(1)=0°632 and F(o)=o0, we have, 
when z=a4, 
b— dy = 4 pa? x 0°826. 
If as usual we set Y= 2j7(¢)—¢), we have from inspection of Fig. 2 log p= —0°5 
at S=150 pe or p=1-1513 at s/a=1. Then 
pop = 1°394[ F(2/a+ 1) — F(s/a—1)—2F(1)] s>a | 2 
= 1°394[F(1 + 2/a)+ F(1—2z/a)—2F(1)] z- <a{ / (5-2) 

Points from equation (5.2) are plotted as crosses in Fig. 3. The model 
represents the potential deduced from A-star densities extremely well, and 
appears to indicate that the effective width of the spiral arm in which the Sun 
is embedded is of order 2a or 300 parsecs. Since 

pb = 2uj*(by — b) = 8 uj? pa? x 0-826 at z=a, 
we have 
p=40°5(a%yj2)! = 40°s(a4,2) 
if j is in (km/sec)! and a in parsecs. This gives the following table : 


median velocity M J, (km/sec)~ p, solar masses per 
km/sec cubic parsec 
6°8 0°070 0°37 
6:0 0080 0°28 


a2 0°0gO0 0°22 








BE ees TARTAR RN tO roy 3 











As Actes 


No. 2, 1957 perpendicular to the galactic plane 209 


the densities being higher than before; since the attracting figure is now limited 
in extent, it requires a greater density to produce a given potential than does the 
unlimited figure. 


Z porsec 
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Fic. 3.—Theoretical distributions of A stars with x, 
@ Double Gauss -| Single Gauss 
< Square cylinder — Curve representing observations shown in 
Fig. 2. 


6. The distribution of radial velocities as actually observed differs from the 
true distribution on account of observational errors. In this section we make 
an attempt to determine an appropriate correction so as to infer the true 
distribution from the observed distribution. 

If a number of stars are members of an open cluster they may be supposed 
to have a common velocity with a very small dispersion about it. For example, 
the dispersion in velocity in the stars in the Pleiades is less than 1 km/sec.* 
Nevertheless the radial velocities of the A stars in the Pleiades in Wilson’s 
catalogue, differenced from the mean radial velocity of such stars, show a median 
of 2-6 km/sec and a root mean square of 3:0 km/sec: that these differences are 
almost entirely due to observational error and not due to a true velocity dispersion 
in the cluster is shown by the proper motions. 

Besides the Pleiades, there are three other clusters (Hyades, Coma Berenices 
and Praesepe) which give useful data about the dispersion of the catalogued 
A star radial velocities. Wilson assigns four qualities to the determinations, 
a, b, c and d, in descending merit, eighty per cent pf the stars being classed 
either b or c: but the clusters do not show a difference of quality between the 
b and c class velocities. 

* R. v. d. R. Woolley, M.N., 116, 296, 1956, quoting spectroscopic material given by Struve and 
proper motions given by Hertzsprung and by ‘Titus. 
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Collected results are as shown in Table IV. 


TaBLe IV 


Dispersion in catalogued radial velocities of A stars in open clusters 


Cluster Pleiades Praesepe Coma Hyades 
No. of A stars in Wilson’s 
catalogue 35 18 11 27 
median velocity difference, 
km/sec 2°6 2°8 2'0 1°9 
root mean square, km/sec 
velocity difference 3°0 3°9 3°0 3°3 


The median refers to the difference, without regard to sign, between the observed 
velocity of the star and the mean velocity of the A stars in the cluster, except in the case 
of the Hyades, where the observed velocity is differenced from the calculated v cos A 
velocity according to van Bueren. 


The evidence for a possible difference between the dispersions in class b 
and class c velocities is inconclusive, but as far as it goes it does not support any 
such difference. Only four stars in the Hyades and one in Coma are of quality c 
(and none in either case is of quality d). In Praesepe the c class velocities show 


less dispersion than do the b class velocities. In the Pleiades there are 13, 15 and _ 


7 velocities of classes b, c, and d, with practically identical dispersions. In all 
cases comparatively few stars of classes Ao and Ar occur. In the Hyades and 
Praesepe the A stars are late in type. In the Pleiades the dispersion increases 
from early A to late A type: the median is 2:0 km/sec for 23 stars of types Ao, 
Ar and A2, compared with 2-6 km/sec for 35 stars of all A types. 

From these results we conclude that a root mean square dispersion of + 3 km/sec 
on account of observational error can be ascribed to the catalogued A star velocities. 
With this figure we compute the following corrections : 


Observed median velocity, km/sec 6 6°5 7 
Observed value of j, (km/sec)~} 0°0795 0'0734 00681 
Corrected value of j, (km/sec)! 0°0748 00698 00692 
Correction to the deduced density: multiply by 0°89 0°90 0°92 


The deduced density is accordingly reduced by about ten per cent, by allowing 
for the observational dispersion in radial velocities, taken to be +3 km/sec. 
(The reduction would be approximately doubled if the dispersion were taken 
as +4 km/sec.) 

We now summarize the data concerning the median velocity, as follows: 

156 A stars of all types in Oort’s galactic cap regions 6-7 km/sec, 

115 Ao stars within 40° of the galactic poles 7:2 km/sec, 

55 Ao stars within 40° of the galactic poles a+b quality only 6-6 km/sec, 

If the adopted median velocity is taken as 6-7 km/sec, and it is assumed that the 
observational dispersion has a r.m.s. value of 3-0 km/sec, then 


j=0°0747 (km/sec) 
and with this value of j we find 


Py =0°18 solar masses per cubic parsec. 
This assumes that the velocity distribution is a single Gaussian distribution 
and that the material is stratified in parallel planes. 
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This density is high in comparison with values given by J. H. Oort* which 
range from o-11 to 0-08 solar masses per cubic parsec, though Lindbladf gave 
0°22 in these units and Jeans{ 0-14. ‘The value given here might be reduced if 
evidence could be adduced for the reality of a double Gaussian distribution in 
A star velocities, but would have to be increased if any allowance was necessary 
for departure from stratification of galactic material in parallel planes: and this 
departure must occur. The writer is inclined to the view that if evidence from 
stellar classes other than type A forces us to adopt a density smaller than, say, 
o-15 solar masses per cubic parsec it will be necessary to suppose that the A stars 
do not satisfy Jeans’s relation. This would be the case, for example, if they 
were an expanding system, but it does not seem to be an hypothesis that one 
would lightly adopt. 

Royal Greenwich Observatory 

Herstmonceux Castle, Sussex: 

1957 January 1. 


* J. H. Oort, B.A.N., 6, No. 238, 284, 1932. 
+ B. Lindblad, Upsala Medd., No. 11, p. 30, 1926. 
t J. H. Jeans, M.N., 82, 122, 1922. 





GENERAL RELATIVITY AND MACH’S PRINCIPLE 


W. Davidson 
(Communicated by W. H. McCrea) 


(Received 1957 January 25) 


Summary 


Sciama, among others, has taken the view that general relativity has failed 
to account satisfactorily for the inertial properties of matter. This paper 
shows that general relativity is entirely consistent in principle with Sciama’s 
ideas of inertia as an inductive effect predominantly of distant matter, and that 
therefore his remarks concerning general relativity are not justified. It is 
shown that general relativity provides a superior presentation of his idea of 
Mach’s principle and appears to be the general tensor theory he was looking 
for. 

Arguments are put forward to show that general relativity may fully incor- 
porate Mach’s principle contrary to Einstein’s own belief. This paper 
emphasizes the fitness of the steady state theory as a cosmological solution 
which permits this possibility. 





1. Introduction.—A tentative theory has been presented by Sciama (1), with 
Maxwell type equations, which is designed to provide a combination of Newton’s 
laws of motion and of gravitation with the inertial frames determined by Mach’s 
principle. In the introduction to his paper Sciama states that general relativity 
has failed to provide an adequate theory of inertia. He claims that his theory 
differs from general relativity principally in the following respects: 

(i) it enables the amount of matter in the universe to be estimated from a 
knowledge of the gravitational constant ; 

(ii) the principle of equivalence is a consequence of his theory, not an initial 
axiom; and 

(iii) it implies that gravitation must be attractive. 

The chief characteristic of Sciama’s theory is that ‘‘in the rest frame of any 
body the gravitational field of the universe as a whole cancels the gravitational 
field of local matter so that in this frame the body is ‘free’. ‘Thus in this theory 
inertial effects arise from the gravitational field of a moving universe.” For 
this purpose Sciama employs a scalar potential ® and a vector potential A to 
calculate gravitational effects, using Maxwell type field equations in flat space- 
time. 

It is the purpose of this paper to show that general relativity is fully consistent 
with this interpretation of Mach’s principle by Sciama, and to indicate that general 
relativity may fully incorporate Mach’s principle. 

2. Free particle in general relativity.—The motion of a free particle in general 
relativity, when the gravitational field is weak and when the reference frame is 
such that the spatial velocity of the particle is small compared with the velocity 
of light, can be described by a Maxwell-type pondermotive equation. This idea 
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is not new and has in fact, with limited application, been presented by Einstein 
(2). But since Einstein’s derivation of the result appears to contain errors of 
detail we give our own derivation here, before investigating its significance for 
Mach’s principle. 

Let Latin letters indicate space coordinates running over indices 1, 2, 3 while 
Greek letters cover the space-time indices 1, 2, 3, 4. The world line of a free 
particle in general relativity is a geodesic in the field having equations, in a 


standard form, 
d dx® Lag Ax* dx? 
(8 a) 7} ‘i xs (1) 
where 
ds* = g4,(dx")? + 2g, dx'dx” +g,,,dx? dx". (2) 
For «=? equations (1) may be written 


as \ed (de?) Bean (db \8d (do) te de Auge de 
dxt}] ds\°" ds = Ox! dx'} ds Sis Ty ds ox' dx! ; ‘ax’ dx’ dxt 


Write now x*= t, dx?/dt=v", and neglect squares and products of the spatial 
coordinate velocities 7”, getting 


ds d dt \_ ,&u _dsd . dt\ | Ba .. 
dt dt\ Sv ds” )~ 2 Ox’ — didt\ Fas) * Oxi °° (3) 


We shall now represent the metric (2) as that of a weak field in the form 
ds? = (1 + y44)dt? + 2yqydx? dt — (1 — y4,)(dx")? — (1 — y29)(dx*)? — (1 — y33)(dx)P? 
+ Ypydxrdx' (p49) (4) 
Here we take the velocity of light, c, as unity. We see that the y,, are the 
deviations of the g,, from the Galilean values in the so-called inertial frames. 
They are the y,,, of Einstein’s treatment of the problem except for sign due to his 
employment of imaginary x". 
We make the assumption that the squares and products of the y,,, and those 
of their derivatives can be neglected. In solving the field equations to this 
approximation Einstein showed (3) that the y,,, were the solutions of the equations 


a \3 a \2 a\2 /a\2 ml 
(ar) + (a) _— jase ere 


provided ge g2 
os- =i 5% v5) + . aeront 7 — 33% vp) = (6) 


to the order of the approximation. Here «=87G/c* where G is the Newtonian 
come of gravitation, and y%=8”y,, where 5% are the Galilean values of the 








ax"ox® 


g’®, Assuming the contribution of stress to the energy momentum tensor to be 


vanishingly small compared with the density and momentum components for 
the case considered by Einstein, equation (5) yields Einstein’s solution: 


[pe] ay 


[p] K 
yu=Ya=0=Yu= - {7 —dV, Yin= == | y 


4" (7) 


and 
Ypa =0, p -q 
In this solution p is the mass density, u” the space velocity, of the element of 
mass in volume dV at distance r from the point where the y,,, are evaluated, all 
quantities being measured by observers at rest in the reference frame. Square 





214 W. Davidson Vol. 11 


brackets indicate retarded values corresponding to the propagation of the field 
with the unit velocity. 

The integrals, supposed convergent, are over all matter producing the field. 
Such a solution of the wave equation of Lorentz (equation (5)) is well known to 
be valid only if the quantities solved for (the y,,,) tend to zero in a suitable way. 
Thus it is clear that the case considered by Einstein involves mass concentrations 
only in the neighbourhood of the space origin and a field metric which is Galilean 
at “infinity”. The integrals in (7) evaluated over such mass concentrations 
are therefore evidently convergent. ‘The solution has to: be consistent with 
conditions (6) which will be satisfied if the expressions 


0 
ax (V7 ie 58% y5) 


vanish, for all yx, to the first order in the y,,. Using (7) and the fundamental 
equations 7,” =o, it is easily seen that for integrals over a finite region of mass 
these expressions are indeed second order quantities. 
To this approximation therefore, retaining only first order terms, we can 
reduce equation (3) to 


d ( giv" y d , 
ai ( Si) = 2 8 - 5 (ga) + Sater, 
at V £44 Ox 


on using (2) to find the appropriate approximation for ii ineachterm. There 
is of course no summation over 7 on the left. Rearranging we can write: 











5(- se.) =-} = oe: = — 84) + 15a (g4;) — = (21,) bo» (8) 

so that 
re rade) = 28 ~ 5(-e) +H Gw)- Ow bor.) 
Write now bm wit pee Ar=— 46 (A) lav (10) 


so that by (7), 11 = Y22 = Y33 = Yaa = 29, Yap = — A”. (11) 


Equation (g) can therefore be written in vector form covering 7=1, 2, 3 
d 0A 
it ((1 —3¢)v¥} = —grad$— > +va curl A. (12) 


Equation (12) is the required Maxwell-type pondermotive equation of the 
field. ‘The assumptions made during its derivation are: 
see The particle velocity v in the reference frame is assumed small such that 

2 is negligible compared with v/c. 
“iy The deviations of the g,, from the Galilean values are small such that 
their squares and products and those of their derivatives can be neglected. 

(iii) The deviations y,,, vanish at “infinity” so that the quantities A, ¢ are 
defined in terms of convergent integrals. If in addition we now further assume 
that 

(iv) The source velocities of the field are also small in the reference frame so 
that the same remark as in (i) applies for them, then equation (12) reduces to 


a 
t a 


d 
<= ~ grad $ OA ev vcurlA, (13) 
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The equation obtained by Einstein was (our notation) , 
d oA 
ii [(1 —¢)v] = —gradd— my TVA curl A, 


Since he assumed condition (iv) as well as (i), (ii), (iii), his result is incorrect to 
the order he was considering, and misleading. In obtaining this result he put 


(ps, i] =1 (2 - Se) 


An A : 
ox? ox' 


thereby neglecting the term og;,,/0x* which, when p=i, contributes to our result 
in equation (g) as the term d(—y,;)/dt in the coefficient of v‘ in the left hand side. 
The neglect of this term is of course consistent with condition (iv), but on the 
other hand the retention of the term d(—¢)/dt in the coefficient of v, arising in 
our approximation from the term d (— }y,,)/dt in the coefficient, is not consistent 
with Einstein’s assumptions. 

3. Interpretation of the pondermotive equation.—As Einstein pointed out, 
equation (12) indicates that general relativity goes far towards incorporating 
Mach’s principle. It may be compared with the Newtonian equation, viz., 

: = —grad ¢. 
The additional terms are small in the quasi-Galilean frame considered by Einstein, 
and, as he said, beyond physical measurement. Nevertheless they show in the 
sense of Mach’s principle how concentrated matter affects the inertial mass of a 
freely moving particle, and the acceleration of its locally inertial rest frame relative 
to the given frame, in the following respects: 

(i) The inertial mass is apparently proportional to 1 — 3¢. 

(ii) The locally inertial rest frame of the particle is accelerated by means of : 

(a) gravitational attraction towards the local mass concentrations indicated 
by the term —grad¢; 

(5) an inductive effect of local accelerating matter in the same sense as the 
acceleration, indicated by the term —0A/ot; 

(c) an inductive effect of matter which is rotating relative to the compass of 
inertia (to use Gédel’s phrase) at “‘infinity’’, in the sense of the rotation, as indi- 
cated by the term vacurlA. This is of the same type as the “‘fictitious’’ Coriolis 
force familiar in Newtonian dynamics, when a reference frame is used which is 
rotating relative to the compass of inertia. Centrifugal force also arises in this 
case as a fictitious. gravitational force. 

It is clear therefore that general relativity certainly incorporates in detailed 
manner the aspects of Mach’s principle indicated above. For a satisfactory theory 
of Mach’s principle, however, Einstein realised the necessity of showing how inertia 
depended on the entire cosmic distribution of matter. Because he assumed that 
the metric was Galilean at “infinity” and therefore excluded any contribution to 
A and ¢ other than that of the matter concentrated near the space origin, he was 
unable to examine the cosmic influence on inertia. 

We shall here put forward an analysis to show that general relativity actually 
permits the same interpretation of inertia which has been presented by Sciama 
as the inductive effect of the whole universe. 
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4. Inductive effect of the universe in general relativity.—We shall investigate 
the extent to which we may generalize the circumstances when the motion of a 
free particle may be described by a Maxwell-type pondermotive equation. For 
this purpose we make the assumptions less restrictive than in the quasi-Galilean 
case as follows: 

(i) The particle velocity v in the reference frame is assumed small such 
that v?/c? is negligible compared with o/c. 

(ii) The velocities of the sources of the field, in the region of space-time 
coordinates with which we shall be concerned, are also small of the same order 
so that the same remark applies. 

(iii) The deviations of the g,, from the Galilean values are small in the above 
quoted range of space-time coordinates, such that their squares and products 
and those of their derivatives can be neglected. We do not however assume that 
these deviations vanish at “infinity”, nor that they even remain small outside the 
specified range. 

It is clear from equation (8) that the equation of motion of a free particle can 
in these circumstances be written 

v' a a e é 
or 1S 5a +{5,(e0)- Sle) ber (14) 





for i= 1, 2, 3. 

This equation is generally covariant, in the sense that, in all reference frames and 
regions of space-time which do not violate the assumptions above, it describes 
the space motion of the free particle in terms of the derivatives of the y,,, involved. 


We now generalize the quantities A, ¢ occurring in the quasi-Galilean analysis 
by defining 


(A, ®) = (—g4 3814)- (15) 
The three equations in (14) may then be written concisely 
WY = —grad@— + vacurlA: (16) 


The vector notation implies the vector character of the terms for purely spatial 
transformations. For space-time transformations however the quantities (A, ) 
do not transform as a 4-vector but as components of the tensor g,,._ This is 
because, unlike the corresponding electromagnetic pondermotive equation, the 
permitted transformations are not necessarily between inertial frames and 
therefore not in general linear. 

It is to be noted that here we have not as in the quasi-Galilean case identified 
A, ® with the deviations of the g,, involved, from their Galilean values but, 
consistent with our endeavour to account for the whole of inertia according to 
Mach’s principle, in terms of the total g,,._ The covariance of (16) is secured by 
the tensor character of the total g,, involved; the deviations do not transform as 
tensors for general transformations. Indeed according to the field equations it 
is the toral g,,, field that is related inseparably to the distribution of matter in the 
whole universe. 

Bearing in mind therefore the physical interpretation of the quantities A, ¢ 
in the quasi-Galilean case we should expect analogous interpretation of A, ® in 
(16) which would, if Mach’s principle is to be satisfied, take account of the 


distribution and motion of matter in the whole universe, relative to the particular — 


wt, 





owl 
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reference frame being used. It would be an immediate consequence of such an 
interpretation of the terms in (16) that the “fictitious” forces of Newtonian 
mechanics in accelerating or rotating reference frames would become directly 
attributable to the inductive effect of a moving universe. 

In particular, in a reference frame in which a freely moving particle was 
permanently at rest, equation (16) would reduce to 


A 


oA 
— grad — — ay =o (17) 


holding at the particle. This is the equation postulated by Sciama. To use 
Sciama’s expression the “‘gravoelectric” field of the whole universe would be 
zero at the particle and it would be gravitationally “free” in its own rest frame. 

For a reference frame at rest relative to the averaged motions of the rest of the 
matter in the universe (the “smoothed-out” universe) we should expect by Mach’s 
principle that, in the neighbourhood of the space origin, the derivatives of A, O 
on the right of (16) would vanish and therefore that the left hand side must vanish. 
The real existence of such frames which are locally inertial is the basis of 
Newtonian mechanics. ‘This aspect of Mach’s principle is built into general 
relativity theory since the field equations predict that such a reference frame will 
be Galilean near the space origin, because of the spherical symmetry about it. 
Thus in this neighbourhood the metric will approximate to 


ds? = dt? — dx? — dy? — dz. (18) 


It is emphasized that in this paper we attach importance for Mach’s principle 
to the total g,, involved in equation (16) and not just their derivatives. General 
reasons for this have already been given and further justification provided in 
Sections 5, 6. Accordingly it is important to obtain the total value of ®. It 
follows from equation (18) that the static potential ®), at the origin of such a 
frame, of the whole universe would be 


®, = dgqs(0) = } (or 4c? in general units). (19) 


The dimensions of ® and the significance we are trying to associate with it would 
require D, to be of order —GM/R where M is the effective gravitational mass of 
the universe and R its effective radius. Sciama’s approach is to define @, as 


™=RodV 

oe ° : ° ° 

-| —— where a is the gravitational mass density, and he gets 

r 
T=0 


G®, = —2. 


Both results are numerically of the same order. ‘The discrepancy in sign will 
occupy us later. Before investigating to what extent general relativity theory 
justifies this tentative physical interpretation of A, ®, we give some applications 
of our theory. 

5. Applications of the inductive theory in general relativity.—{i) Sciama 
considers the case of a free particle in rectilinear motion in the gravitational field 
of a mass M which is at rest relative to the smoothed-out universe. If we choose 
a reference frame at rest relative to the smoothed out universe with this mass M 
at the space origin, we can neglect in that neighbourhood the deviations from the 
Galilean values of the g,, as far as they arise from the universe as a whole, and 

15* 
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include only the deviations due to the mass M. ‘Thus to this approximation the 
metric will be 


tthe (: * xo dt? — (: + a) (dx? + dy? + da®). (20) 





r 


It is to be noted that, according to the ideas presented in this paper, the contri- 
bution to the g,, potentials from the universe as a whole is present in the Galilean 
terms of the g,,. 

Since the universe is at rest in this frame we have 


A=0 
while P=} (: ime =) . (21) 
J 





r 


Suppose the particle is moving freely towards the mass M along the x axis. 
If its space coordinates are (x,,0, 0) at coordinate time ¢ then its coordinate speed 
is dx,/dt= —v, where v>o. Make now the transformation to a suitable rest 
frame for the particle, by means of the relations 

x=X+x,y=Y,2=Z,t=T (22) 
yielding dx =dX —vdT, dy=dY,dz=dZ,dt=dT. We get therefore to sufficient 


order for the covariance of (16) 


ds =(1 wen sem \aT* + 20d XAT -(1 + —) (dX24+dY2+4dZ%). (23) 





Thus in the particle’s rest frame 


A=(-—v7,0,0) | 





=} (: -o- =) ' (24) 
2 - | 
Apply now equation (17) in the particle’s rest frame, yielding 
a 2GM a 
=- — — — — {—-v)=— 
ax{i(t v = )} ap v)=0 (25) 
leading to 
GM dv 
ont ar i + 77) =0 (26) 


on substituting the original coordinates. This is the Newtonian equation of 
motion of the particle and is also the equation which would follow from the general 
pondermotive equation (16), applied in the original frame, using (21). 

On examining (24) and (25) we see that the origin of the inertial term 


- (-2) 


in (25) lies in the relative motion of the universe, yielding A=(-—v, 0,0) in the 
particle’s rest frame, and thus creating an inductive field at the particle which 
balances the local gravitational attraction due to the mass M, thus connecting 
with Sciama’s ideas. 

We note also that the A, ® in (24) arise by transformation of the whole g,, and 
not just their deviations from the Galilean values, in accordance with our 


tentative interpretation of the Galilean values as the static potentials of the whole 
universe, 
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(ii) The other case considered by Sciama is that of a particle moving with 
uniform motion in a circle under the attraction of a mass M at the centre, this 
mass being again at rest relative to the smoothed-out universe. 

Transform therefore from the metric (20) to a suitable rest frame for the 
particle according to the relations 


aT 


x=XcoswT— Y sinw7 
y=YocoswT+XsinwT 





ae (27) 
t=T 
so that to sufficient order 
ds? =(1 -2GM/R —w* R®) dT? — 2w (— YdXdT +XdYd7) 
—(1+2GM/R) (dX? +dY?+dZ?) (28) 
with R?=X?+ Y?. 
Thus in this frame 
A=(-—w/Y,wX,0) 
o=1(1— 26M _ aan), (29) 
. R 
‘The equation (17) then yields 
- Se +w?R=0 (30) 


which is the Newtonian equation of motion, and also, putting R=r, what 
would be given by (16) in the original frame. 

Connecting with Sciama’s ideas we say that the gravitational attraction by M 
is balanced by the gravitational field induced by a rotating universe, whose 
rotational momentum is indicated by A in (29). 

(iii) As a final example we shall show how, by means of the covariance of (16), 
the Newtonian ‘“‘fictitious” forces may be attributed to the inductive effect of a 
moving universe in the most general Newtonian motion of the reference frame 
relative to a locally inertial frame. 

Consider a free particle at rest in a reference frame which is locally inertial, 
so that the metric is approximately as given by (18) in that region. Let r be the 
position vector of the particle in that frame. Then by (15), (16) we have 


r=constant. (31) 


‘Transform to a second frame whose space origin has variable velocity V and 
which has variable spin w relative to the first frame. If the position vector of 
the particle in this frame is R, then a well-known kinematic result of Newtonian 
motion gives 

r=V+R+waR (32) 


#=V+wraV+2waR+@aR+wa(waR)+R (33) 


differentiation being with respect to the common Newtonian time of either frame. 
Thus for the particle in the second frame 


* 


= —[V+waV+20aR+@aR+wa(waR)]. (34) 
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The transformation connecting the two frames is, by (32), in differential form 


ir =(V R)dT+dR 
aWVowaRitoi } Gs 
Hence ds? = dt? —dr* 

=[1 —V2—2V.(w a R)—(w a R)*]d7? — 2(V + w « R). dRdT — dR’ (36) 


so that in the second frame 


A=V+waR (37) 
® = }[1 —V?—2V.(waR)—(wa R)*] 37 
Now grad ® =waV+wa(waR) 


oA av ow 
oT 6OT)~=—shCOOT 
=V+waR 


AR 


while curl A= 2w. 
Hence by (16) R= —[V+waV+2waR+@aR+Wa(warR)] 
giving complete agreement with (34). 

Thus our theory gives an exact treatment of the fictitious forces as the 
inductive effect of a moving universe, 

6. Physical interpretation of A, ® in general relativity.—he analysis in this 
section is intended to be of a tentative nature, since complete rigour cannot be 
claimed for it. 

In Section 4, equation (19), we obtained the result 

Dy = 20° = 3844(0) 
for that value of the gravitational potential ® of the whole universe which enters 
into the pondermotive equation (16), when evaluated at the space origin of a 
reference frame locally inertial there. In order to interpret this result in terms of 
Mach’s principle we recall the expressions for ¢ in the quasi—Galilean case given 
by (7) and(11). Since the field equations are relations for the whole g,,, in terms 
of the matter in the whole universe, we make the tentative inference that in some 
way the Galilean terms themselves are related to world gravitation, so that inertia 
would arise in accordance with Mach’s principle. ‘To what extent does general 
relativity provide justification of this inference ? 

In all cosmological models of general relativity in which the average inertial 
density p does not vanish there is an effective radius R of the model which is the 
distance, measured in a suitable way, to the horizon of the model where the 
velocity of the matter relative to the space origin equals the velocity of light. 
For an observer at the origin matter which goes beyond this distance virtually 
ceases to exist because of the Doppler effect on its light and presumably on its 
gravitation. ‘l’hisis the case whether the model be of the homogeneous rotating 
type (Gédel’s models) or the isotropic expanding or contracting types. We 
shall discuss the latter as an example. ‘These have the general metric 


(dr? + r°d@ + r? sin? dd¢*) 





ds* = c?dt® — eM) 


(+7 4RS? ” 


where R,” may be positive, negative, or infinite, and the fundamental particles 
have constant 7, 0, ¢. 
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‘Lhe distance to the particle at (r, 6, 6), measured in the simultaneity of the 
fundamental observers at cosmological time t, from r=o is 
Tr 
Jato [" —__o . 
; a ol +r/4R,? 
Its radial velocity is therefore 


i=} gl (39) 
so that | /|=c when /= on 
lg| 
Thus R=+ z according as g20 (40) 
. cl 
and l=+ R: (41) 


The importance for Mach’s principle is that R is related to o, the cosmological 
density of gravitational mass. In general relativity theory, gravitational mass 
density is defined so as to lead to Gauss’ flux theorem for small regions of space 
(see, for example, Synge (4)), and for the isotropic cosmological models 
a=p+3p/c* where p is the pressure. The gravitational “force” on unit mass 
due to the field is in this case the proper acceleration relative to the space origin. 
With these definitions McCrea (5) has shown that the equations of general 
relativity for the isotropic models are consistent with the variation of the gravi- 
tational force according to the Newtonian inverse square law, using proper 
radial distance but Euclidean geometry, for spatial regions of any size. ‘The field 
equations applied to the metric (38) give, with A=o, in general units, 


87G Ce 3 





= — —¢ oO _ 5 252 
a ee 
302 3 (42) 
yy Det) 4. 2 G52 
87Gp = R° no 4 4° j 
so that 
8nGo= — 3(# + $2. 
Now for g>o0, g= = j= —2cR/R® by (40). Hence 
82Go= — oa (1 ~ Ric) 
whence, if g>0, GoR?= — ti ~Rjc) 
ae (43) 
and, if g<o, GoR*® = — a + R/c) 


which are the required relations between o and R, at time t. For g>o we see 
that o20 according as R>c. If R>c matter is entering the region bounded by 
the defined horizon; if R <c matter is passing beyond this horizon. Forg<o,a<0 
according as Rée. 
By (41) and (43) we get the Newtonian type equation 
|= —4$nGol (44) 


relating gravitational force and proper distance at time f. 
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Comparing equation (44) with the pondermotive equation (16) we see that, 
for a fundamental observer whose radial space coordinate is the proper distance / 
and whose time is the cosmological time t, A=o and —grad ®= — z mGol. 
Equation (44) however holds for all /<R and not just in the neighbourhood of 
the origin. Consider therefore the gravitational ‘“ work” done by the field when 
a particle of unit mass is moved from its actual position at time ¢ to the horizon 


and therefore beyond influence of the origin. This will be 
R 


hep t2G | oldl. (45) 


I 

This may be regarded as the analogue of the Newtonian potential at the distance /, 
in the gravitational field as witnessed by an observer at the origin. Both o and 
R will vary with / in this integral as the motion proceeds, according to (42), (43). 
However , may be evaluated as 


R. 3 
o,=| idl = 42-4 


= —(1-P/R?). (46) 
This is the potential at 1 at time ¢. 
Putting /=o we get @, =c?/2 (47) 
which therefore provides a physical identification, in a natural way, of the 
potential ®, arising in equation (19). 

The result given by (47) for ®, is got as the limit of ®, when / tends to zero 
irrespectively of the sign of o or g. For instance if o>o and g>o then R>c, 
so that for the field to carry the particle to the horizon of the space origin would 
mean going backwards in time. The discrepancy in sign between our ®, and 
Sciama’s, referred to at the end of Section 4, arises because of Sciama’s arbitrary 
definition of ® for an expanding universe. His definition appears to ignore the 
above considerations and in particular to presuppose the identity of the cosmo- 
logical gravitational mass density and the inertial density. 

For cogent reasons which have been put forward elsewhere (6) a stationary 
cosmological solution is to be preferred. ‘The only known stationary solution 
which does not contradict observational results (expansion, spatial isotropy) is 
the steady state theory proposed by Bondi and Gold (6). This has the metric 

ds? = c2dt® — e%t/B( dy? + 726? + 1? sin? Odd?) (48) 
where R is a constant which is the effective radius for the model. ‘The steady 
state model, unlike the general cosmological models of metric given by (38) for 
which equation (16) vanishes identically, allows a static metric to be used so that 
the motion of a particle, relative to the observer at the origin, is measured by the 
rate of change of the spatial coordinates. This is the De Sitter metric 

2 
ds* = c*(1 — [?/R®)dr? — ae — Pd@ — I? sin? 6d¢? (49) 
connected to (48) by a well known transformation. In this metric / is the distance 
from the origin, in the simultaneity of the fundamental observers, of our general 
analysis. The theory of Section 4 defines the ® involved in (16) as }g4y 
which in the case of the metric (49) gives 


O= ae ~ P/R?) (50) 
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agreeing with (46) and therefore having the physical interpretation associated 
with (46). 

It is to be noted that the steady state forms a natural cosmological background 
to mass concentrations. For instance the exact solution of the field equations 
for an isolated mass m superimposed on the steady state is 





‘ 2mG-sSCosi2 : dl oe ee ee 
at e(1 — ye — 7) | 2mG — /?d@? — *sin*6d¢?. (51) 
cl R 
For this metric 
2 G 
@=<(1-P/R)- = (52) 


to be interpreted physically as the work done by the field in removing unit mass 
from the point in question to the horizon of the model, regarding mG/R as 
negligible. 

A Newtonian type integral for ® in terms of the distribution of the mass does 
not follow simply in the case of the general models because of the stated dependence 
of o and R on cosmological epoch. However, for the steady state, 0 and R are 
constant and we may write for the potential of unit mass, at distance / from the 
mass odV constantly in the volume element dV, 


-R , 
a.-of Za 


J, E 
I I 


-R 
Thus ®,= — io | (; - R) 4nldl 





Equation (43) gives for the steady state GoR* = — 3c?/47 so that 

®, =c?/2 
in agreement with (47) and justifying our physical interpretation of ®, as the 
gravitational potential of all the matter in the universe apparent to an observer 
at the origin and having influence there. 

The quantity A of our theory defined in equation (15) of Section 4 is zero for 
the cosmological metric (38). On making a transformation such as that of the 
first example in Section 5, equation (22), a non-zero A arises by transformation of 
the g,,. If we accept the association of the Galilean values of the g,, with world 
gravitation, according to the tentative analysis presented above, the association 
of A with the relative momentum of the universe would also follow. While the 
Galilean g4,, viz. c?, is associated with ®) as 2, the spatial Galilean g,, are 
associated with A. Thus in the first example in Section 5 to get Ain the particle’s 
rest frame we have to multiply the Galilean g,, in the original frame, viz. 
—1=g,,(0), by v. According to the subsequent application of the ponder- 
motive equation —g,,(0) is proportional to the inertial mass of the particle, the 
whole equation indicating equality of gravitational and inertial mass in the case 
of a particle. Since —g,,(0)=g44(0)/c? = 2 /c? we see that inertial mass can 
therefore be associated with the influence of the whole universe, in accordance 
with Mach’s principle. 
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7. Comparison with Sciama’s theory.—In this final section we shall remark 
briefly on the three principal differences claimed by Sciama between his theory 
and general relativity, enumerated (i), (11), (ili) in the introduction to this paper. 

(i) It is evident from the analysis in this paper that a knowledge of G, 
occurring in the integral (45) leading to (47), together with R given as cT where 
T is the reciprocal of the Hubble parameter, leads to an estimate of the amount 
of matter in the universe in general relativity as well as in Sciama’s theory. 

(ii) Sciama states that in general relativity the principle of equivalence 
predicts that one gravitating mass in an otherwise empty universe produces the 
same inertial effects as in his theory, and since there is no universe in this case 
to give rise to the inductive field “‘it is difficult to see why the principle of equi- 
valence should be true’’. Such an argument however implies a solution of the 
field equations involving the use of coordinates for all points of space-time in a 
universe which, except for the isolated mass, is empty. Such coordinates are 
purely conceptual, defined without reference to matter and restoring to space an 
objective substance, independent of matter, which general relativity has sought 
to deny. The logical course for general relativity, according to the field 
equations, is to relate the Galilean g,, of special relativity to world gravitation in a 
full universe. That general relativity may be capable of doing so has been 
indicated in this paper, where, independently of the value of p as long as it does 
not vanish, } g4,(0) has been identified as ©) = c?/2, the potential of the universe. 
The case of the empty universe can only logically be approached as a limit where 
p->o and R-> o (equation (43)), so that inertia is always accounted for. 

(iii) Sciama’s determination of the sign of the field in his theory is of doubtful 
significance as it depends on his ® as defined turning out to be negative. Reasons 
for questioning this arbitrary definition of ® have been given in Section 6. 

In general relativity the coefficient of 7” in the field equations is chosen so 
as to make gravitation attractive on the small scale (the pressure being then 
relatively negligible). However, on the cosmological scale this leads to gravi- 
tational mass being interpreted as negative if the expansion is accelerating 
(equation (44)). Thus there would appear to be no intrinsic importance to be 
attached to the sign of the field in any theory that allows for factors at present 
unknown, that is whether the gravitational effect of a lump of matter is more 
primary than that of a cosmological region containing ‘‘ zero-point”’ stress. 

I am grateful to Professor W. H. McCrea for opportunities for discussion 
of this paper and for his criticisms which have stimulated several improvements. 
Acknowledgments are due also to the authorities of the Battersea College of 
Technology for a relaxation of lecturing duties to perform this research. 
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SURVEY OF BRIGHT GALAXIES SOUTH OF —35° DECLINATION*# 
G. de Vaucouleurs 


Summary 


A survey of the bright galaxies listed in the Shapley-Ames Catalogue 
(H.A., 88, 2, 1932) south of declination —35° was carried out between 
1952 and 1955 with the 30-inch Reynolds reflector of the Commonwealth 
Observatory, Mount Stromlo. About 250 one-hour exposures were 
obtained on Eastman 103a-o emulsion with the reflector diaphragmed to 20 
inches (f/6, scale 67”/mm). New or revised types and dimensions were deter- 
mined for 210 S—A objects (out of 230 listed) and in addition for 120 other 
NGC and IC objects and also for 130 uncatalogued objects, a total of 460 
galaxies generally larger than 1’ and/or brighter than the 14th magnitude. 

Morphological types were assigned in a revised classification scheme 
incorporating the improved Mt Wilson—Palomar system of Hubble and 
Sandage, as modified and supplemented by the present Mt Stromlo work. 
It includes four main classes: ellipticals (E), lenticulars (SO), spirals (S) 
and irregulars (I); two families of lenticulars and spirals : ordinary (SA) and 
barred (SB), with intermediate objects (SAB); two main waricties of lenti- 
culars and spirals: ringed (r) and spiral (s), with intermediate objects (rs); four 
or five sub-divisions or stages along the spiral sequences; early (a), intermediate 
(b), late (c, d) and magellanic (m), the latter forming the transition toward 
the magellanic irregulars I (m). Three stages, early, intermediate and late 
are also distinguished in the lenticular sequences (SO-, SO, SO*) between 
ellipticals and spirals. Outer ring-like structures (R) are frequently observed 
near the transition stage SO/a between lenticulars and spirals (Section 3 and 
Plate I). The apparent frequency of tvpes among the S—A objects in the zone 
is: E(23-4 per cent), SO (21-0 per cent), SA (24°4 per cent), SB (26-3 per cent), 
I (4-9 per cent), (Section 6). The possible evolutionary significance of the 
classification sequence is briefly discussed (Section 16). 

Dimensions obtained by visual measurement of the negatives and reduced 
to standard conditions are given for the bright, inner regions (D,;vd,) and for 
the faint, outer regions (Doxdy), (Section 4). The median outer diameters Dy 
for the S—A objects in the zone are: E (1'4), SO (2’0), S (3’2), (Section 7); 
there is no perceptible latitude effect (Section 10). However, apparent 
diameters depend on apparent flattening, i.e. tilt, and a correction, indepen- 
dent of type, must be applied (Section 9). The ratio of inner to outer 
dimensions (D,/D)), measuring the concentration, depends also on flattening 
(Section 8); the variations of corrected diameter and of concentration index 
vs. galaxy type are determined (Section 11). 

The axial ratio (d/D) of ellipticals is discussed (Section 12); there is no 
excess of either globular or elongated shapes compared with a random 
distribution of true oblateness and tilt. - On the contrary, lenticulars (Section 
13) show a large excess of strongly elongated shapes as expected; but spirals 
(Section 15) of corrected diameters larger than 3’ show a deficiency of elon- 
gated shapes which may indicate some restriction on the spatial orientation 
of the nearer galaxies. 


* The full text of this paper appears in Memoirs of the Commonwealth Observatory, Mount 
Stromlo, 3, No. 13. 
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Summary of Mount Stromlo Memoir Vol. 117 


Ring-structures in spirals are of special interest as possible geometric 
distance indicators. A preliminary discussion (Section 17) of 35 objects 
indicates that the mean diameter of the inner ring (r) is constant throughout 
the spiral sequences from SO/a to Scd; the range among objects of known 
velocities appears to be less than 2 to 1 with a mean value of about 2:5 kilo- 
parsecs and a variance of about + 15 per cent. Ring diameters appear 
potentially capable to serve as useful distance indicators out to about 100 
megaparsecs. ‘The fainter outer rings (R) are about 1°7 times larger on the 
average but are often too faint to be detected in small objects. 

Relative distances based on apparent diameters are determined for seven 
clouds or regions outlined by the apparent distribution of the S—A objects 
in the zone (Section 5 and Plate II). An excellent agreement is observed 
between the various determinations based on ellipticals, lenticulars, spirals 
and rings. Absolute distances based on ring diameters are derived (Section 
18) which vary from about 7 megaparsecs for a group of spirals in Grus 
(region VII) to 10 megaparsecs for the Pavo—Indus cloud (region VI). 
A dozen foreground objects of outstanding apparent diameters are at distances 
of from 1 to 3 megaparsecs. ' 

Relative diameters at unit distance are determined as a function of type 
along the classification sequence (Section 19). The range of mean diameters 
from late spirals to ellipticals is only 2 to 1 in contrast with early Mt Wilson 
results, but in good agreement with Harvard data on the Virgo cluster. The 
frequency functions of diameters give for most types some indication of a 
bimodal distribution with a small group of supergiant and a main group 
of giant galaxies; dwarf systems are practically excluded by the magnitude- 
surface brightness selection of the S—A Catalogue. 

Pairs and triplets are classified according to the visibility or otherwise 
of tidal distortion in the components (Section 20). For 30 pairs showing no 
clear sign of interaction the mean projected separation is 18 kiloparsecs 
(range 8 to 45 kpc); for 8 pairs showing distortion effects it is 6-5 kiloparsecs 
(range 3 to8 kpc). The mean critical separation within which tidal distortion 
becomes notable is of the order of 9 kiloparsecs (actual). Tzo colliding 
systems (NGC 1487, NGC 3256) are described; neither of them is a radio 
source. Some dense groups and clusters were noted in Fornax, Antlia, 
Centaurus, Telescopium and Pavo (Section 21). The densest, around 
NGC 4696, appears to be associated with a radio source. 

Comparisons with earlier surveys including objects south of —35° indicate 
notable under-exposure for the older Helwan. and Harvard (H.A., 88, 2) 
series; there is good agreement, however, with the more recent Harvard 
series (H.A., 88, 4, 1934), especially for spirals; the mean scale factor is 
HA 88,4/Stromlo (D,) = 1-2 and the average deviation + 8 percent. Reduction 
factors to the scale of the Heidelberg survey of northern nebulae (Reinmuth’s 
Catalogue) are obtained via H.A., 88, 4. 

Four lists give general information and details of the Reynolds reflector 
data on (I) S-A objects south of —35° (main catalogue), (II) additional 
NGC and IC objects, (III) large uncatalogued objects, (IV) pairs of galaxies, 
in the same zone. 

Eight plates give (1) a description of the classification system, (II) a map 
of the apparent distribution of the S—A objects south of -35°, (III to VI) 
Reynolds reflector photographs of 12 large spirals, and (VII, VIII) photo- 
graphs of 32 other bright galaxies illustrating the classification system and 
some interesting objects. 
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Summary of Mount Stromlo Memoir Vol. 117 


Ring-structures in spirals are of special interest as possible geometric 
distance indicators. A preliminary discussion (Section 17) of 35 objects 
indicates that the mean diameter of the inner ring (7) is constant throughout 
the spiral sequences from SO/a to Scd; the range among objects of known 
velocities appears to be less than 2 to 1 with a mean value of about 2°5 kilo- 
parsecs and a variance of about + 15 per cent. Ring diameters appear 
potentially capable to serve as useful distance indicators out to about 100 
megaparsecs. ‘The fainter outer rings (R) are about 1°7 times larger on the 
average but are often too faint to be detected in small objects. 

Relative distances based on apparent diameters are determined for seven 
clouds or regions outlined by the apparent distribution of the S—A objects 
in the zone (Section 5 and Plate II). An excellent agreement is observed 
between the various determinations based on ellipticals, lenticulars, spirals 
and rings. Absolute distances based on ring diameters are derived (Section 
18) which vary from about 7 megaparsecs for a group of spirals in Grus 
(region VII) to 10 megaparsecs for the Pavo—Indus cloud (region VI). 
A dozen foreground objects of outstanding apparent diameters are at distances 
of from 1 to 3 megaparsecs. 

Relative diameters at unit distance are determined as a function of type 
along the classification sequence (Section 19). The range of mean diameters 
from late spirals to ellipticals is only 2 to 1 in contrast with early Mt Wilson 
results, but in good agreement with Harvard data on the Virgo cluster. ‘The 
frequency functions of diameters give for most types some indication of a 
bimodal distribution with a small group of supergiant and a main group 
of giant galaxies; dwarf systems are practically excluded by the magnitude- 
surface brightness selection of the S-A Catalogue. 

Pairs and triplets are classified according to the visibility or otherwise 
of tidal distortion in the components (Section 20). For 30 pairs showing no 
clear sign of interaction the mean projected separation is 18 kiloparsecs 
(range 8 to 45 kpc); for 8 pairs showing distortion effects it is 6-5 kiloparsecs 
(range 3 to8 kpc). The mean critical separation within which tidal distortion 
becomes notable is of the order of 9 kiloparsecs (actual). Tzvo colliding 
systems (NGC 1487, NGC 3256) are described; neither of them is a radio 
source. Some dense groups and clusters were noted in Fornax, Antlia, 
Centaurus, Telescopium and Pavo (Section 21). The densest, around 
NGC 4696, appears to be associated with a radio source. 

Comparisons with earlier surveys including objects south of -35 indicate 
notable under-exposure for the older Helwan and Harvard (H.A., 88, 2) 
series; there is good agreement, however, with the more recent Harvard 
series (H.A., 88, 4, 1934), especially for spirals; the mean scale factor is 
HA 88,4/Stromlo (Dy) = 1-2 and the average deviation + 8 per cent. Reduction 
factors to the scale of the Heidelberg survey of northern nebulae (Reinmuth’s 
Catalogue) are obtained via H.A., 88, 4. 

Four lists give general information and details of the Reynolds reflector 
data on (I) S-A objects south of —35° (main catalogue), (II) additional 
NGC and IC objects, (III) large uncatalogued objects, (IV) pairs of galaxies, 
in the same zone. 

Eight plates give (1) a description of the classification system, (II) a map 
of the apparent distribution of the S—A objects south of —35°, (III to VI) 
Reynolds reflector photographs of 12 large spirals, and (VII, VIII) photo- 
graphs of 32 other bright galaxies illustrating the classification system and 
some interesting objects. 
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